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Chapter 1 



On the Dynamics of Large 
Particle Systems in the 
Mean Field Limit 



The general purpose of all the mean field equations considered below is to de- 
scribe the dynamics of a very large number of identical particles, assuming that 
the interaction between particles is known exactly. For instance, at a temper- 
ature of 273 K and a pressure of 1.01 • 10'' Pa, the number of molecules of any 
ideal gas to be found in a volume of 2.24 • 10~^ m^ is the Avogadro number, i.e. 
6.02 • 10^^. This is typically what we have in mind when we think of a number 
of particles so large that the dynamics of each individual particle cannot be 
exactly determined. Thus, in order to be of any practical interest, these mean 
field models should be set on (relatively) low dimensional spaces. In any case, 
they should involve fewer degrees of freedom (of the order of a few units) than 
the total number of degrees of freedom of the whole particle system. 

To be more precise: in classical mechanics, the number of degrees of free- 
dom of a single point particle moving without constraint in the d-dimcnsional 
Euclidean space R'' is d. The single-particle phase space is the set R'' x R"^ 
of pairs of all possible positions and momenta of an unconstrained single point 
particle in R''. 

For a system of N identical point particles moving without constraint in the 
rf-dimensional Euclidean space R'', the number of degrees of freedom is therefore 
dN. The space (R'' x R'')^ of 2A^-tuples of all possible positions and momenta 
of the N point particles is the N -particle phase space. 

Thus, the laws of motion of classical mechanics (Newton's second law) writ- 
ten for each molecule of a monatomic gas enclosed in a container form a system 
of diflferential equations set on the A''-particle phase space, where is the total 
number of gas molecules to be found in the container. With N of the order of 
the Avogadro number, this is of little practical interest. In the kinetic theory 
of gases, the evolution of the same gas is described by the Boltzmann equation. 
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an integro-differential equation set on the single-particle phase space. Although 
the kinetic theory of gases will not be discussed in these notes, all the mean 
field limits considered below will involve the same reduction of the A'"-particle 
phase space to the single-particle phase space. 

Here is a (by no means exhaustive) list of typical examples of mean field 
equations: 

(a) the particles are the ions and the electrons in a plasma; the interaction is 
the Coulomb electrostatic force; in that case, the mean field equation is the 
Vlasov-Poisson system of the kinetic theory of plasmas; 

(b) the particles are the nuclei and the electrons in a molecule; the interaction 
is again the Coulomb electrostatic force; the corresponding mean field model in 
this context is the Hartree equation or the system of Hartree-Fock equations in 
atomic physics; 

(c) the particles are vortices in an incompressible fluid in space dimension 2; the 
interaction is given by the Helmholtz potential; the corresponding mean field 
model is the vorticity formulation of the Euler equations of incompressible fluid 
mechanics 

All these models are obtained as approximations of the system of equations 
of motion for each particle in the limit when the number of particles involved 
tends to infinity. Rigorous justifications of these approximations are based on 
various mathematical formalisms that are introduced and discussed below. 

Excellent references on the issues discussed above are [HHl HOI HI] 

1.1 Examples of mean field models in classical 
mechanics 

In this section, we consider the examples mentioned above in which the motion of 
each particle is described in the formalism of classical mechanics, i.e. examples 
(a) and (c). All these examples have a common structure, to be discussed later, 
which involves the Liouville equation recalled below. 

1.1.1 The Liouville equation 

The Liouville equation governs the evolution of the distribution function for a 
system of particles subject to an external force field. 

The notion of distribution function is fundamental in the kinetic theories of 
gases and plasmas, and more generally in statistical mechanics. It was intro- 
duced by Maxwell in one of his most famous article^ 

The distribution function of a system of identical point particles is / = 
f{t, X, v), that is the number density of particles that are located at the position 

^Illustrations of the Dynamical Theory of Gases, Philosophical Magazine (1860); reprinted 
in "The Scientific Papers of James Clerk Maxwell", edited by W.D. Niven, Cambridge Uni- 
versity Press, 1890; pp. 377-409. 
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X and have instantaneous velocity v at time t. In other words, the number of 
particles to be found at time t in an infinitesimal volume dxdv of the single- 
particle phase space centered at (x, v) is ~ /(t, x, v)dxdv. 

Assume that a particle located at the position x with instantaneous velocity 
V at time t is subject to some external (or imposed) force field F{t,x,v). As 
a consequence of Newton's second law of classical mechanics, the distribution 
function / of the system of particles under consideration satisfies the Liouville 
equation 

dtf + v- + ^ dW,{F{t, X, v)f) = , 

where m > is the particle mass. The Liouville equation is a partial differential 
equation (PDE) of order 1, whose solution can be expressed by the method of 
characteristics. 

To the PDE of order 1 above, one associates the system of ordinary differ- 
ential equations (ODE) 

{x = v, 

\v^^F{t,X,V), 



with the usual notation 



These ODEs are referred to as the "equations of characteristics" for the Liouville 
equation. Denote by i i— > {X{t,toTX,v),V{t,tQ,XTv)) the solution of this ODE 
system such that 

X{t(3,t(^,x,v) = X , V{to,to,x,v) ^v; 

the map {t, to, x, v) i— )■ {X{t, to, x, v), V(t, to, x, v)) will be henceforth referred to 
as "the characteristic flow" associated to the Liouville equation above. 

One immediately recognizes in the system of ODEs above the equations of 
motion of classical mechanics for a particle of mass rn subject to the force field 
F = F{t, X, v) (the first equality being the definition of the velocity, while the 
second is Newton's second law). 

Assume that the force field F is such that the characteristic fiow {X, V) is 
globally defined (i.e. defined on R x R x R"^ x R'*). This characteristic flow is 
used as follows to express the solution of the Cauchy problem 

r dtf + v - V^f + ^AW,{F{t,x,v)f) ^{) , x,^;eR'^xR^ teR, 

For each test functiord cj) £ C^(IV^ x R''), and for each i € R, one has 

/(t, X, w)0(x, v)dxdv f "(x, v)(t>{{X, V){t, 0, x, v))dxdv , 



■^For each topological space X and each finite dimensional vector space E on R, we des- 
ignate by C{X,E) the set of continuous functions defined on X with values in E, and by 
Cc{X, E) the set of functions belonging to C(X, E) whose support is compact in X. For each 
n,k > 1, we denote by C^(R", E) the set of functions of class C'° defined on R" with values 
in E whose support is compact in R". We also denote C{X) := C(X, R), Cc(X) := Cc(X, R) 
and C^*^(R") := C^*={R",R). 
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and this completely determines the distribution function f(t,x,v). 

Exercise: Deduce from the equality above an explicit formula giving f{t,x,v) 

in terms of /*" and of the characteristic flow {X, V) . 

For a concise discussion of the method of characteristics (and a solution of 
the exercise above), see chapter 1, section 1, in ,17 . 

1.1.2 The Vlasov-Poisson system 

Our first example of a mean field kinetic model is the Vlasov-Poisson system 
used in plasma physics. 

Consider a system of identical point particles with mass m and charge q. 
The electrostatic (Coulomb) force exerted on any such particle located at the 
position X € R'^ by another particle located at the position y x is 



where eo is the dielectric permittivity of vacuum. This force is indeed repulsive, 
as it is of the form X{x — y) with A > 0. 

More generally, the electrostatic force exerted on a particle located at the 
position X by a cloud of particles with number density p{t, y) (which means that, 
at time approximately p{t,y)dy particles are to be found in any infinitesimal 
volume element dy centered at y) is 



X — y 



47reo |a; - y\ 



3 ' 




Recall that 



G{x) 



1 



\{0} 



47r|a;| ' 



is the solution of 




'0 



in I?'(R3) , 

as |x| — >■ cx) , 



where Sq designates the Dirac measure at the origin. 
Thus the electrostatic force F is given by 



F{t,x)=qE{t,x), 



where E = E{t, x) is the electric field, i.e. 



E{t,x) 



with electrostatic potential cf) given by 



£0 
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In particular 



-A(b{t, •) = A(G * Pit, •)) = --(AG) * Pit, •) 
Co eo 

= —So* pit,-) = —pit,-). 
eo eo 

The Vlasov-Poisson system describes the motion of a system of identical 
charged point particles accelerated by the electrostatic force created by their 
own distribution of charges — referred to as the "self-consistent electrostatic 
force" . It consists of the Liouville equation governing the evolution of the distri- 
bution function, coupled to the Poisson equation satisfied by the self-consistent 
electrostatic potential as follows: 

( idt + v- V.)/(t, X, v) - ^V^(t>it, x) - V„/(<, x,v)^Q, 



A^(/)(t,a;) = j-qpfit,x) , 



Pfit,x) = / fit,x,v)dv. 

JR3 



Exercise: Let / = fit,x,v) and (j> = 4>it,x) be a solution of the Vlasov- 
Poisson system. To avoid technicalities, we assume that / and belong to 
C°°(R+ X R3 X R3) and C'°°(R+ x R^) respectively, and that, for each t > 
the functions ix,v) i— fit,x,v) and x i-?- 4>it,x) belong to the Schwartz clasf|j 
5(R^ x R^) and 5(R^) respectively. 
1) Prove that 

dtpfit, x) + divx vfit,x,v)dv = 0, 
dt / mvfit,x,v)dv + divx / mv®'^ fit,x,v)dv 

JR3 Jr3 

-fgp/(i,x)V^0(t, a;) = 0, 
\m\v'\^ fit,x,v)dv + dwx / v\m\v'\^ fit,x,v)dv 



R3 



R3 



-Vx4>it,x) 



R3 



qvfit, X, v)dv = . 



These three equalities are respectively the local conservation laws of mass, mo- 
mentum and energy. 



^For each n > 1, the Schwartz class iS(R") is the set of real-valued C°° functions defined 
on R," all of whose partial derivatives are rapidly decreasing at infinity: 



5(R") := {/ e C°°(R") s.t. \xrd°'f{x) -s> as |x| -s> oo for all m > 1 and a G N"} . 
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2) Prove the global conservation of mass (or equivalently of the total number of 

particles) 

-j: f{t,x,v)dxdv = 0. 

3) Prove that, for each ij) e C^(R^), one has 

AV'VV' = div((VV')®^) - IVIV^I^ 

and conclude that 

^ r r 

mvfit, X, v)dxdv 



(This is the global conservation of momentum). 
4) Prove the global conservation of energy: 

4(// ^m\vff{t,x,v)dxdv + ^eo[ \V^(j){t,x)fdx] = . 

\iiR3xR3 iR3 / 

(Hint: the term 

/ (p{t, x)divx\ / qvf{t,x,v)dv\dx 

can be simplified by using the result in question 1).) 

The Vlasov-Poisson system above is written in the case of a single species 
of identical charged particles. In reality, a plasma is a globally neutral system, 
and therefore involves many different species of particles. Denoting by fk = 
fk{t,x,v) for k = 1,...,M the distribution function of the kth species in a 
system of M different kinds of particles, the Vlasov-Poisson describes the motion 
of this system of particles submitted to the self-consistent electrostatic force 
resulting from the distribution of charges of all the particles in that system: 

( {dt + v \7^)fk{t, x, v) - ^ Vx(/>(t, x) ■ Vyfkit, x,v) = 0, fc = 1, . . . , M , 



M 



- A^(l){t, a;) = ^ ^ qkPk{t, x) , 



fe=i 



pk{t,x) = / fk{t, 
JR3 



x, v)dv . 



Here ruk and qu designate respectively the mass and the charge of particles 
of the A:th species. In practice, considering systems of particles of different 
species in the Vlasov-Poisson system does not involve additional mathematical 
difficulties. Therefore, we shall consider only the (unphysical) case of a single 
species of particles in our mathematical discussion of this system. 
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There is a huge hterature on the Vlasov-Poisson system, which is of con- 
siderable importance in plasma physics. The global existence and uniqueness 
of classical solutions of the Cauchy problem for the Vlasov-Poisson system was 
obtained by Ukai-Okabe [94] in the 2-dimensional case, and independently by 
Lions-Perthame [681 ^-nd PfafFelmoser [T^ in the 3-dimensional case. For a pre- 
sentation of the mathematical analysis of this system, the interested reader is 
referred to [1T1|T7||H1]. 

1.1.3 The Euler equation for two-dimensional incompress- 
ible fluids 

The Euler equation for an incompressible fluid with constant density (set to 1 
in the sequel without loss of generality) governs the evolution of the unknown 
velocity field u = u{t, x) G and of the unknown pressure field p = x) £ R 
in the fluid, assumed to fill the Euclidean plane R^. It takes the form 

dtu{t, x) + {u ■ 'Vx)u{t, x) + VxP{t, a;) , div^^ u{t, x) = . 

The pressure field p can be viewed as the Lagrange multiplier associated to the 
constraint div^: u — 0. 

There is another formulation of the Euler equation in which the pressure 
field is eliminated. To the velocity field u = u{t,x) G R^, one associates its 
scalar vorticity field u = bj{t,x) defined as follows 

ui{t,x) := dx^U2it,x) - dx2Ui{t,x) . 

The vorticity field satisfies 



One can check that the Euler equation is (formally) equivalent to the system 



There is yet another formulation of this system, where the velocity field u is 
represented in terms of a stream function. Indeed, since div^; u = 0, there exists 
a stream function (j) = 0(i, x) such that 



dfUj + diva;(uaj) = . 



dtU! + divx{uuj) — , 
< divx u — , 

UJ = dlYx{Ju) , 



where J designates the rotation of an angle — ^: 




u{t,x) = J\7^(j)[t,x) . 



The vorticity field lu is related to the stream function as follows: 
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so that the vorticity formulation of the Euler equation becomes 

dtid + div2:(a; JVj:^) = , 

- Ax(f> = LO ■ 

In this last system, the unknown is the pair (w, 0); once the stream function <f> 
is known, the velocity field u is obtained as the spatial gradient of cf) rotated of 
an angle — ^. 

There is an obvious analogy between this formulation of the Euler equation 
and the Vlasov-Poisson system: the vorticity field w is the analogue of the 
distribution function / in the Vlasov-Poisson system, the stream function is the 
analogue of the electrostatic potential in the Vlasov-Poisson system, and the 
velocity field u is the analogue of the electric field in the Vlasov-Poisson system. 
There is however a subtle difference between both systems: the equations of 
characteristics associated to the vorticity formulation of the Euler equations are 

JX = -v.A{t,x), 

while those corresponding with the Vlasov-Poisson system are 

X = -^V^^{t,X), 
m 

(after eliminating V — X). The first ODE system is of order 1, while the second 
is of order 2, because the gradient of the solution of the Poisson equation is a 
velocity field in the vorticity formulation of the Euler equation, while it is an 
acceleration field in the case of the Vlasov-Poisson system. 

An excellent reference on the Euler equations of incompressible fiuids and 
on issues related to the mean field limit is [71] . 

1.1.4 The Vlasov-Maxwell system 

Observe that magnetic effects are not taken into account in the Vlasov-Poisson 
system. Yet charged particles in motion generate both an electric and a magnetic 
field; on the other hand, the electric and magnetic fields both act on charged 
particles via the Lorentz force. Studying magnetized plasmas (such as tokamak 
plasmas, or the solar wind, for instance) requires using the Vlasov-Maxwell, 
instead of the Vlasov-Poisson system. 

The unknown in the Vlasov-Maxwell system is the triple {f,E,B), where 
E = E{t,x) e R'^ and B = B{t,x) G R'^ are respectively the electric and the 
magnetic fields, while / = f{t,x,^) is the particle distribution function. Notice 
that the velocity variable v is replaced with the momentum ^ in the distribution 
function. In other words, f(t, x, ^) is the phase space density of particles located 
at the position x at time t, with momentum ^. The momentum and the velocity 
are related as follows. The relativistic energy of a particle of mass m with 
momentum ^ is 

e(C) := V"i2c4 + c2|^|2, 
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where c is the speed of hght (notice that e(0) = mc^, the energy at rest of a 
particle with mass to). Then 

viO = Ve(0 = , ""'^ . 

^ Vto2c4+c2|^|2 

The (relativistic) Vlasov-Maxwell system takes the form 

' dtf + v{0 ■ V,/ + q{E + viO X B) • V^/ - , 
div^ B = 0, dtB + rot:r E = Q, 

diva; E = ^qpf , - rot^c B = -^ioqjf , 

The first equation in this system is the relativistic Liouville equation governing 
/; the term 

q{E{t,x)+v{Ox B{t,x)) 

is the Lorentz force field at time t and position x exerted by the electromagnetic 
field (E, B) on a particle with charge q and momentum ^. The second and 
third equations are respectively the equation expressing the absence of magnetic 
monopoles and the Faraday equation, while the fourth and fifth equations are 
respectively the Gauss equation (as in electrostatics) and the Maxwell- Ampere 
equation. Since the source terms in the Maxwell system are the charge density 
qpf and the current density qjf, the Lorentz force in the Liouville equation is 
the self-consistent Lorentz force, i.e. the electromagnetic force generated by the 
motion of the charged particles accelerated by this force itself. 

For the same reason as in the case of the Vlasov-Poisson system, the case of 
a single species of charged particles is somewhat unrealistic. Physically relevant 
models in plasma physics involve different species of particles so as to maintain 
global neutrality of the particle system. 

Exercise: Following the discussion in the previous section, write the relativistic 
Vlasov-Maxwell system for a system of M species of particles with masses to^ 
and charges qk, for k = 1, . . . , M. Write the local conservation laws of mass, mo- 
mentum and energy for the resulting system, following the analogous discussion 
above in the case of the Vlasov-Poisson system. 

For more information on the Vlasov-Maxwell system, see [HI [T71 [53] . 
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1.2 A general formalism for mean field limits in 
classical mechanics 

We first introduce a formalism for mean field limits in classical mechanics that 

encompasses all the examples discussed above. 

Consider a system of N particles, whose state at time t is defined by phase 
space coordinates zi{t), . . . , ZN{t) € R**. For instance, Zj is the position Xj of 
the jth vortex center in the case of the two dimensional Eulcr equations for 
incompressible fiuids, and the phase space dimension is d = 2. In the case of 
the Vlasov-Poisson system, the phase space is x ^ R^, so that d = 6, 
and Zj = {xj,Vj), where xj and vj are respectively the position and the velocity 
of the jth particle. 

The interaction between the ith and the jth particle is given by K{zi,Zj), 
where 

K : R*^ X R'* -i> R'' 

is a map whose properties will be discussed below. 

The evolution oi zi{t), . . . , ZN{t) G R'' is governed by the system of ODEs 

dz- ^ 

^(*) = E % (0) , iJ = l,...,N. 

3 = 1 

Problem: to describe the behavior of Zi{t), . . . , 5jv(i) € R"^ in the large A'' limit 

and in some appropriate time scale. 

First we need to rescale the time variable, and introduce a new time variable 
t so that, in new time scale, the action on any one of the N particles due to the 
A'' — 1 other particles is of order 1 as AT — )• +oo. In other words, the new time 
variable t is chosen so that 

dz- 

= 0(1) for each i = 1, . . . , iV as iV oo . 

dt 

The action on the ith particle of the — 1 other particles is 

N 

^K{zi,Zj), 

3 = 1 

and it obviously contains A^ — 1 terms of order 1 (assuming each term K{zi, zj) 
to be of order 1, for instance). Set i = t/N, then 

dz 1 ^ 

j=i 
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From now on, wc drop hats on all variables and consider as our starting 
point the rescaled problem 

1 ^ 

3 = 1 
3Vi 

At this point, wc introduce an important assumption on the interaction 
kernel: the action of the jth particle on the ith particle must exactly balance 
the action of the ith particle on the jth particle. When the interaction is a 
force, this is precisely Newton's third law of mechanics. Thus we assume that 
the interaction kernel satisfies 

K{z,z') = -K{z',z) , z,z'gK'^. 

We have assumed here that the interaction kernel K is defined on the whole 
j^d ^ j^d gpace; in particular, the condition above implies that K vanishes 
identically on the diagonal, i.e. 

K{z,z) = 0, zeR'^. 

Hence the restriction j ^ i can be removed in the summation that appears 
on the right hand side of the ODEs governing the A^-particle dynamics: since 
K{zi{t), Zi{t)) = for alH = 1, . . . , A'', one has 

1 ^ 



At this point, we can explain the key idea in the mean field limit: if the 
points Zj{t) for j = 1, . . . , A^ are "distributed at time t under the probability 
measure f{t,dz)" in the large N limit, then, one expects that 

1 ^ r 

^.y^K{zi{t),Zj{t))^ / K{zi{t),z')f{t,dz') asN^+oo. 

This suggests replacing the A^-particle system of differential equations with the 
single differential equation 

z{t)= I K{z{t),z')f{t,dz'). 

Here f{t,dz) is unknown, as is z{t), so that it seems that this single differential 
equation is insufficient to determine both these unknowns. 

But one recognizes in the equality above the equation of characteristics for 
the mean field PDE 

dtf + div,{flCf)=Q, 
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where the notation JC designates the integral operator defined by the formula 

ICf{t,z) := / K{z,z')f{t,dz'). 

Now, this is a single PDE (in fact an integro-differential equation) for the single 
unknown /. 

A priori / is a time dependent Borel probability measure on R'', so that 
the mean field PDE is to be understood in the sense of distributions on R''. In 
other words, 

^[ <j>iz)f{t,dz)= [ ICf{t,z)-Vc^{z)f{t,dz) 
at jjid j^d 

for each test functiorH (j) G CK'R'^). 

A very important mathematical object in the mathematical theory of the 
mean field limit is the empirical measure, which is defined below. 

Definition 1.2.1 To each N -tuple Zn = {zi,...,zn) E (R"*)^ ~ R'^^, one 
associates its empirical measure 

1 ^ 

J = l 

The empirical measure of a iV-tuple G (R'')^ is a Borel probability 
measure on R''. As we shall see in the next section, the iV-tuple 

t^ZN{t) = (zi(i),...,zw(t)) 

is a solution of the TV-particle ODE system 

1 ^ 

if and only if the empirical measure ^ZN(t) ^ solution of the mean field PDE 
3tMz«(t) + div^(^z«(t)^/^Ziv(t)) = ■ 

^For each topological space X and each finite dimensional vector space E on R, we denote 
by Cj, (X, i?) the set of continuous functions defined on X with values in E that are bounded 
on X. For each n,k > 1, we denote by C^(R", E) the set of functions of class C* defined on 
R." with values in E all of whose partial derivatives are bounded on R": for each norm \ ■ \ e 
on E, one has 

C^iR",E) := {/ e C''(R",E) s.t. sup \d"f(x)\E < oo for each a e N"} . 



We also denote Ct{X) := Ct{X, R) and C*^(R") := C*={R", R). 
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We conclude this section with a few exercises where the reader can verify that 
the formalism introduced here encompasses the two main examples of mean- 
field theories presented above, i.e. the two dimensional Euler equation and the 
Vlasov-Poisson system. 
Exercise: 

1) Compute A In |a;| in the sense of distributions on (answer: 2tt5o). 

2) Define 

^ J{x — x') 
\x — x'l^ 

where J designates the rotation of an angle — ^ : 

'0 1' 



" \-l 

For each w = uj{t,x) belonging to C^(R+ x R-^) such that supp(w(t, •)) is com- 
pact for each t > 0, prove that the vector field u defined by 

u{t,x) := / K{x,x')uj{x')dx' 



is of class Cl on R+ x R'^ and satisfies 

divj; u{t, x) — , diva;( Ju)(i, x) = a;(i, x) . 

3) Conclude that the two dimensional Euler equation for incompressible fluids 
can be put in the formalism described in the present section, except for the fact 
that the interaction kernel K is singular on the diagonal of R^ x R^ . 
Exercise: Let (/, (p) be a solution of the Vlasov-Poisson system satisfying the 
same assumptions as in the exercise of section [1.1. 21 Assume further that 

f{Q,x,v)dxdv — l, and // vf{0,x,v)dxdv = 0. 

1) Prove that 

f{t,x,v)dxdv — l and // vf{t,x,v)dxdv — for alH > . 

R3xR3 JjR3xR3 

2) Set z — (x, v) and 

if (z, z') = K{x, V, x' , v') := 

Prove that 



K{x, V, x', v')fit, x', v')dx'dv' = {v, -^V,</.(t, x)) , 

R3xR3 



where 



-A^(f>{t,x) = — f{t,x,v)dv. 
eo 7r3 



3) Conclude that the Vlasov-Poisson system can be put in the formalism de- 
scribed in the present section, except for the fact that the interaction kernel K 
is singular on the set {{x,v, x' ,v') G (R'^)'* s.t. x — x'}. 
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1.3 The mean field characteristic flow 

Henceforth we assume that the interaction kernel K : IV^ x R"* R'' satisfies 
the following assumptions. 
First K is skew-symmetric: 

(HKl) K{z, z') = -K{z', z) for all z, z' eR'^ . 

Besides, K G C^{IL'^ x R''; R''), with bounded partial derivatives of order 1. 
In other words, there exists a constant L > such that 

{HK2) snp \VzK{z,z')\< L, and sup \V ^, K{z, z')\ < L . 

Applying the mean value theorem shows that assumption (HK2) implies that 
K is Lipschitz continuous in z uniformly in z' (and conversely): 

sup \K{zi,z')~K{z2,z')\ < L\zi - Z2I , 
z'e'R'' 

sup \K{z,Zl) ~ K{z,Z2)\ < L\zi- Z2\ ■ 

zeR^ 

Assumption (HK2) also implies that K grows at most linearly at infinity: 

\K{z,z')\<L{\z\ + \z'\), z,z'eR'*. 

Notice also that the integral operator JC can be extended to the set of Borel 
probability measuretH on R'' with finite moment of order 1, i.e. 

ViiR'^) := |p e V(R'^) s.t. J \z\p{dz) < ooj , 

in the obvious manner, i.e. 

ICp{z) := / K{z,z')p{dz'). 

The extended operator K, so defined maps 7'i(R'^) into the class Lip(R''; R'') of 
Lipschitz continuous vector fields on R''. 

With the assumptions above, one easily arrives at the existence and unique- 
ness theory for the A^-body ODE system. 

Theorem 1.3.1 Assume that the interaction kernel K £ C^(R'^ x R'',R'') 
satisfies assumptions (HK1-HK2). Then 

a) for each N > 1 and each N -tuple = (2^1", • . • , z^f}), the Cauchy problem 
for the N -particle ODE system 

N 



[ Z^{0) - Z 



^Henceforth, the set of Borel probabiUty measures on R'' will be denoted by 'P{R'*). 
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has a unique solution of class on R 

t K> Zjv(t) = {zi{t), . . . , ZN{t)) =: TtZ'j^ ■ 

b) the empirical measure f{t,dz) :— fJ-T^z^" weak solution of the Cauchy 
problem for the mean field PDE 

r at/ + div,(//c/) = o, 

Statement a) follows from the Cauchy-Lipschitz theorem. Statement b) fol- 
lows from the method of characteristics for the transport equation. For the sake 
of being complete, we sketch the main steps in the proof of statement b), and 
leave the details as an exercise to be treated by the reader. 
Exercise: Let h = b{t,y) € C([0,r];R'') be such that Dyb e C([0,t];R'^) and 

{H) \b{t,y)\<K{l + \y\) 

for all t e [0,t] and y E R'', where k is a positive constant. 

1) Prove that, for each t G [0, r], the Cauchy problem for the ODE 

iY{s) = b{s,Yis)), 
\Y{t)=y, 

has a unique solution s t-^ Y(s, t, y). What is the maximal domain of definition 
of this solution? What is the regularity of the map Y viewed as a function of 
the 3 variables s, t, yl 

2) What is the role of assumption (H)? 

3) Prove that, for each ii,t2,^3 G [0,t] and y E R"^, one has 

Y{t3,t2,Y{t2,ti,y)) ^Yit3,ti,y) . 

4) Compute 

dtYis,t,y) + b{t,y)-VyY{s,t,y). 

5) Let e C^(R'^). Prove that the Cauchy problem for the transport equation 

r dtfit,y) + bit,y)-Vyfit,y) = 0, 

has a unique solution / e (^^([0, r] x R''), and that this solution is given by the 
formula 

/(i,2/) = r"(r(0,t,2/)). 

6) Let fi™ be a Borel probability measure on R''. Prove that the push-forward 
measured 

^i{t) :=y(t,0,-)#M'" 



^ Given two measurable spaces (X, ^) and (Y,B), a measurable map <!? : {X,jV) — {Y,B) 
and a measure m on {X,A.), the push-forward of m under 3> is the measure on (Y,B) defined 
by the formula 

$#m(S) = m,(<I>-^(B)) , for all B 6 6. 
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is a weak solution of 

{dtH + divydib) = . 
Hint: for <p G C^(R''), compute 



4/ 0(y(t,o,2/)V-(dj/), 



dt 



7) Prove that the unique weak solutioi^^ fi € C([0, r], w ~V{'R'^)) of the Cauchy 
problem considered in 6) is the push-forward measure defined by the formula 

^i{t) :=y(t,o,.)#M^" 

for each t e [0,t]. (Hint: for e C^(R'^), compute 

^(r(O,t,-)#A^(i),0) 

in the sense of distributions on (0,t).) 

For a solution of this exercise, see chapter 1, section 1 of |T7] . 

Our next step is to formulate and solve a new problem that will contain both 
the A^-particle ODE system in the mean-field scaling and the mean-field PDE. 

Theorem 1.3.2 Assume that the interaction kernel K g C^{'R'^ x R'',R'*) 
satisfies assumptions (HK1-HK2). For each S R'' and each Borel probability 
measure fi™ € ViCRJ^), there exists a unique solution denoted by 

R9 t Z(t,C",M'") e R'' 

of class of the problem 

dtZ{t,r,t^n = {iCfimzit,r,f^n), 

^i{t) = z{t,■,^l^-)#^i^-, 
z(o,C",m'") = C"- 

Notice that the ODE governing the evolution of t M- C*", /i*") is set in 
the single-particle phase space R"*, and not in the A^-particle phase space, as is 
the case of the ODE system studied in Theorem ll.3.11 

Obviously, the ODE appearing in Theorem ll.3.2l is precisely the equation of 
characteristics for the mean field PDE. Henceforth, we refer to this ODE as the 
equations of "mean field characteristics", and to its solution Z as the "mean 
field characteristic flow" . 

How the mean field characteristic flow Z and the flow Tf associated to the 
A^-particle ODE system are related is explained in the next proposition. 



^We designate by ui — P{R'') the set 'P(R'^) equipped with the weak topology of probability 
measures, i.e. the topology defined by the family of semi-distances 



runs through Ci,(R'*). 



/ 4>{z)^{dz) — / <j>{z)v{dz) 
JR'' Jr.'' 
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Proposition 1.3.3 Assume that the interaction kernel K G C^(R'' x R'',R'^) 
satisfies assumptions (HK1-HK2). For each = {z™ , . . . , z™), the solution 

TtZ^;} ^{z^{t),...,ZN{t)) 

oj the N-body problem and the mean field characteristic flow Z{t, C™, /i™) satisfy 

z,(t)^Z{t,zr\fizi^), t = l,...,N, 

for all ten. 

Proof of Proposition ll.3.3i Define 

Q{t):=Zit,zr,t,z^.), i = l,...,N. 

Theil 



1 



for all i € R. Therefore, Q satisfies 

C.(t) = ilc^^mc,it)) - - ^ Kicm, W) , * - 1, . . . , ^ , 

for all t e R. Moreover 

Therefore, by uniqueness of the solution of the A'^-particle equation (Theorem 
ll.3.ip . one has 

C^{t)^Z,{t), 

for alH = 1, . . . , iV and all i e R. ■ 



The proof of Theorem 11.3.21 is a simple variant of the proof of the Cauchy- 
Lipschitz theorem. 



*The reader should be aware of the following subtle point. In classical references on distri- 
bution theory, such as I53| , the Dirac mass is viewed as a distribution, therefore as an object 
that generalizes the notion of function. There is a notion of pull-back of a distribution under a 
C°° diffeomorphism such that the pull-back of the Dirac mass at yo with a C°° diffeomorphism 
X : -5> satisfying xi^o) = VO is 

This notion is not to be confused with the push-forward under x of the Dirac mass at S^g 
viewed as a probability measure, which, according to the definition in the previous footnote 
is 

X#^xo = &yo ■ 

In particular 

unless X h^-s Jacobian determinant 1 at xq. 
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Proof of Theorem 11.3.21 Let /i*" e 'Pi(R'^), and denote 



Ci := / \z\ti"\dz)^ 



Let 



X := e C(R'^;R'^) s.t. sup 
which is a Banach space for the norm 

\\v\\x ■■= sup 



zGR" 1 + 1^1 

1 + ■ 



< oo 



By assumption (HK2) on the interaction kernel K, for each v,w ^ X, one 



K{v{z),v{z'))^i'"{dz')- / K{w{z),w{z'))fi'''{dz') 
<L f {\viz) - wiz)\ + \v{z') - u)(z')|)Ai*"(rf^') 



has 



Rrf 

< - w\\x{l + \z\) + L\\v- w\\x / (1 + \z'\)fi"'{dz') 

= L||w-u)||x(l + k| + l + Ci) 
<L\\v-w\\x{2 + Ci)il + \z\). 

Define a sequence (2^,i)„>o by induction, as follows: 

Z„+i(<,C) =C+ r / A'(^n(t,C),^n(i,C')K"(rfC')rfs, n>0, 

Jo JRii 

One checks by induction with the inequality above that, for each n G N, 

((2 + Ci)L|i|)" 



\Zn+lit,-)- Zn{t,-)\\x < 



-\\Z,{t,-)-Zoit,-)\\ 



X 



Since 



l^i(i,C)-CI= / / K{C,C)^^'"{dC)ds 
Jo Jtl'' 

Jq Jr'' 

1*1 

L{\C\+C^)ds<L{l + C,){l + \C\)\t\, 





one has 



\\Zn+l{t,-)-Z„{t,-)\\x < 



((2 + Ci)LK|)"+i 
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Thus, for each t > 0, 

Zn(t, •) — > Z{t, •) in X uniformly on [— r, r] 
where Z g C(R;X) satisfies 



^(^,0 = c + 



JR' 



i^(Z(s,C),^(.s,C'))M™«)rfs 



for alH G R and aU C e R''. 

If Z and Z e C(R; X) satisfy the integral equation above, then 

Z{t, C) - Zit, {K{Z{s, C), Z{s, C')) - K{Z{s, C), c')))m'"«) : 

so that, for all t e R, one has 

This implies that 



\\Z{s,-)-Z{s,-)\\xds 



\\Z{t,-)-Z{t,-)\\x=0 

by Gronwall's inequality, so that Z = Z . Hence the integral equation has only 
one solution Z g C(R;X). 

Since Z e C{R+;X), K e C^{Il'^ x R'*,R'') satisfies (HK2) and fi'" G 
Vii'R'^), the function 

f K{Z{sX),Z{sX'))r{dO 
J If 

is continuous on R. 

Using the integral equation shows that the function t i— > Z{t, C) is of class 
C^ on R and satisfies 

dtZ{tx)^ I K{z{tx),z{tx'W{dC,'), 
z(o,c) = C- 

Substituting z' ~ Z{t, (') in the integral above, one has 



K{Z{t, C), Z{t, C WidO = / K{Z{t, C), z')Z{t, W^idz') 



SO that the element Z of C(R;X) so constructed is the unique solution of the 
mean field characteristic equation. ■ 

References for this and the previous section are [201 ES] • 
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1.4 Dobrushin's stability estimate and the mean 
field limit 

1.4.1 The Monge-Kantorovich distance 

For each r > 1, we denote by Vri'R-'^) the set of Borel probabihty measures on 
R"* with a finite moment of order r, i.e. satisfying 

/ \z\^ii{dz) < oo . 

Given /i, G 'Pr(R-'^), we define n(/x, z^) to be the set of Borel probability 
measures tt on R'^ x R'' with first and second marginals /i and respectively. 
Equivalcntly, for each tt e V(R'^ x R'^), 

TT e n(^, zy) ^ // {(f){x)+^{y))'K{dxdy) ^ ct){x) fi{dx) + ^{y)v{dy) 

for each (j},ip e C(R'') such that 0(z) = 0(|z|'') and tp{z) = Od^T) as \z\ -> oo. 

Probability measures belonging to n(/i, i') are sometimes referred to as "cou- 
plings of n and v" . 

Exercise: Check that, if ji and i' e 7'r(R'') for some r > 0, then one has 
n(^,iy) C VrCR"^ X R''). 

With these elements of notation, we now introduce the notion of Monge- 
Kantorovich distance. 

Definition 1.4.1 For each r > 1 and each ii^v ^ 'Pr(R'^), the Monge-Kantoro- 
vich distance dist MK.rifJ-ji^) between /i and v is defined by the formula 

distMK,r(Mif^) = inf I // \x - y\'''n{dxdy) 

These distances also go by the name of "Kantorovich-Rubinstein distances" 
or "Wasserstein distances" — although the minimization problem in the right 
hand side of the formula defining dist MK,r had been considered for the first time 
by Mongelfl and systematically studied by Kantorovich. 

We shall use the Monge-Kantorovich distances only as a convenient tool for 
studying the stability of the mean field characteristic flow. Therefore, we shall 
not attempt to present the mathematical theory of these distances and refer 
instead to the C. Villani's books [95l |96] for a very detailed discussion of this 
topic. 

However, it is useful to know the following property that is special to the 
case r = 1. 

'■^Mongc's original problem was to minimize over the class of all Borel measurable trans- 
portation maps T : R'' — > R'* such that T#/x = u the transportation cost 

/ \x - T{x)\fi(dx) . 
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Proposition 1.4.2 The Monge-Kantorovich distance with exponent 1 is also 
given by the formula 



dist 



MK,1 



sup 

■/,eLip(R,<i) jRd 
Lip(*)<l 



(j){z)n{dz) 



(j){z)iy{dz) 



with the notation 

\<P{x) - (l>iy)\ 
Lip((/)) := sup ■ ■ 

x^yeR'' F y\ 

for the Lipschitz constant of 4>. 



The proof of this proposition is based on a duahty argument in optimization: 
see for instance Theorems 1.14 and 7.3 (i) in ^5] . 



1.4.2 Dobrushin's estimate 

As explained in Proposition 11.3.31 the mean field characteristic flow contains 
all the relevant information about both the mean field PDE and the iV-particle 
ODE system. 

Dobrushin's approach to the mean field limit is based on the idea of proving 
the stability of the mean field characteristic flow Z{t, C™, /i™) in both the initial 
position in phase space and the initial distribution /i™. As we shall see, the 
Monge-Kantorovich distance is the best adapted mathematical tool to measure 
this stability. 

Dobrushin's idea ultimately rests on the following key computation. Let 
Cr, C2" e R'', and let ^if, M2" S 'Pi(R'^). Then 

z(i,ci,A^r)~^(t,C2,A^r) = ci-C2 

+ / / K{Z{sXi.tJ^'),z')p.^{s,dz')ds 

Jo J-R'i 

- [ f K{Z{sX2,fiT),z')m{s,dz')ds. 
Jo Jw 

Since Hj{t) — Z{t, •, ^™)#^™ for j = 1, 2, each inner integral on the right hand 
side of the equality above can be expressed as follows: 

= / K{z{s,Q,^l^),z{s,c^,^^M'(dO 

Jr-^ 
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for j = 1,2. Therefore, for each coupUng tt*" e ■Pi(/x^"', /Xj"), one has 



so that 

z{t,c,,^iT)-z{t,C2,^lT) = Cl-C2 

{K{Z{sXiyr),Z{sX[.I^T)) 

R<^ 



This last equality is the key observation in Dobrushin's argument, which ex- 
plains the role of couplings of /U™ and jJ^ in this problem, and therefore why 
it is natural to use the Monge-Kantorovich distance. 

After this, the end of the argument is plain sailing. By assumption (HK2) 
on the interaction kernel K, for all a, a', 6, h' G R'*, one has 

\K{a, a') - K{b, b')\ < \K{a, a') - K{b, a')\ + \K{b, a') - K{b, b')\ 
<L\a-b\+L\a' -b'\. 

Therefore 

\z{t,Cuf^t)-z{t,C2,iJ-T)\ 

< ICi - C2I + i / \Z{sXi,t^T) - Z{sX2.ii'i')\ds 







+l/ // |z(s,ci,Mr)-^(s,c2,Mr)k"«dC2Ms. 

Jo JJR'^xR'' 

It is convenient at this point to introduce the notation 

D\i:\{s) := [[ \Z{s,(:[,nT) - Z{s,Cif^T)\n{d(:[dQ 
for each tt e ■Pi(R'* x R**). Thus, the previous inequality becomes 

iz(t,ci,Mr)-^(t,C2,/^r)i<iCi-C2i 

+ l[ \Z{s,(:i,fiT)-Z{s,C2,fiT)\ds + L f D[Tr'"]{s)ds . 
Jo Jo 
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Integrating both sides of the inequahty above with respect to 7r™(dCic!C2) 
leads to 

< // ICi - C2|7r"(dCirfC2) + 2L f D[^'^]{s)ds 

J JU'^xIl'' Jo 

^ D[Tr'"]{0) + 2L f D[Tr"']{s)ds . 
Jo 

By Gronwall's inequahty, we conclude that, for all t E R, one has 

D[n'"]{t) < D[7r™](0)e2^l*l . 

Now we can state Dobrushin's stability theorem. 

Theorem 1.4.3 (Dobrushin) Assume that K e C^(R'' x R'',R'') satisfies 
(HK1-HK2). Let nf,nf £ ViiBJ^). For all t € R, let 

r fii{t)^z{t,-,fi'r)#fi\", 

\^i2{t)^z{t,■,^iT)#^lT, 

where Z is the mean field characteristic flow defined in Theorem \1.3.2l 
Then, for all t G R, one has 

distMK.i{l^i{t),pi2{t)) < e'^l*l distMKAl^T,l^T) ■ 

Proof. We have seen that, for aU nY\ E Vi{'R'^) and all tt™ e n(/xl", 
one has 

L>K"](t) < ^[Tr™] (0)6^^1*1 

for all t G R. 

Since Z{t, •, = ^^ (t) for j = 1, 2, the map 

<ft: (Ci,c2)^(^(t,ci,A^r),^(i,C2,A*r)) 

satisfies 

$t#7r" = ^(t)en(A*i(0,M2(i)) 

for all t £ R, since tt™ G n(^5'\/x^"). 
Thus 

distAfK,i(AiiW>M2W) = inf // ICi - C2|7r(dCic?C2) 



< 



inf . //' |Z(t,Ci,A^r)-^(i,C2,M2")k"(rfCidC2 



inf < e^-^l*! inf L'[7r'"](0) 

7r*"en(/i'j",/j^") 7r'"en(Ai5",^4") 

= 6^^1*1 distMK,l(A'r,A^2") 

which concludes the proof. ■ 

The discussion in this section is inspired from [29]; see also [71]. The in- 
terested reader is also referred to the very interesting paper [53] where Monge- 
Kantorovich distances with exponents different from 1 are used in the same 
context — see also US]- 
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1.4.3 The mean field limit 

The mean field limit of the TV-particle system is a consequence of Dobrushin's 
stability theorem, as explained below. 

Theorem 1.4.4 Assume that the interaction kernel K £ C^(R'^ x R*^) and 
satisfies assumptions (HK1-HK2) . Let f™ he a probability density on TV^ such 
that 

[ < C50. 

Jn-' 

Then the Cauchy problem for the mean field PDF 

dtfit, z) + div, {f{t, z)JCf{t, z)) = , z e R^ t e R , 

r I pin 

j\t=0 ~ J 

has a unique weak solution f £ C(R; L^(R'')). 

For each N > 1, let 5{N) = (z™^, . . . , z]^ j^) £ (R"*)^ be such that 

1 ^ 
=1 

satisfies 

distA/K,! (mscw) , r") -> as N -^oo. 

Let t i-> Tt3(iV) = (zi,7v(i), • • • ,ZAr^Ar(t)) £ (R'^)^ be the solution of the 
N -particle ODF system with initial data "^{N), i.e. 

( I N 

The^ 

l^Tt3{N)^f{t, as N^oo 

in the weak topology of probability measures, with convergence rate 

distMK,iW,3W,/(i,-)^'') < e'^l*l distMKAP3iN)jn ^ 
as N ^ oo for each t £ R. 

Proof. By Theorem II. 3.2l and questions 6 and 7 in the exercise on the method 
of characteristics before Theorem 11.3.21 one has 

f{t, O^'^ = Z{t, •, /™^'')#/"^'^ 

for all t £ H. This implies in particular the uniqueness of the solution of the 
Cauchy problem in C(R; Li(R'^)) for the mean field PDE. 

^"Thc notation Jf'^ designates the Lebesgue measure on R'*. 
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By Proposition [TTX31 

I^T,3{N) = ^(i, •,M3(w))#/^3(Af) 

for all i e R. 

By Dobrushin's stability estimate, 

distMK,i(A'T,3W,/(t,-)=S^') < e2^l*ldistMK,i(M3(iV),r") 
for all t G R, and since we have chosen 3(^) so that 

distA./K,i(M3(A'),r") ^ 
as — ?► oo, we conclude that 

distMK,i(/iTt3(JV),/(i,-)^'') ^ 

as — > oo for each t e R. 

As for weak convergence, pick (f> G Lip(R'^); then 



iz)tJ-n3(N){dz) - (j){z)f{t,z)dz 



< 



{(j){x) - (j){y))TT{dxdy) 
\(t){x) - (j){y)\'K{dxdy) 



R'^xR'' 



< Lip(0) / / |a; - y\-K{dxdy) 

for each tt e U{fiT,5(N), f{t, ■)^'^)- Thus 



R'i 



{z)nTt5{N)idz) - (t>{z)f{t,z)dz 



R<i 



< Lip(<^) 



7ren(/jj 



inf 



\x — y\-K{dxdy) 



Lip((/)) A:iiiiMK,i{p-Tt^(N) J{t,' 







for each t G R as — ?► oo. (Notice that the inequality above is an obvious 
consequence of the definition of distMif,i, so that the equality in Proposition 
11.4.21 is not needed here.) 

This is true in particular for each (p £ C^{Il'^), and since C^(R'^) is dense in 
Cc(R''), we conclude that 



R<i 



(k{z)iJ.Tt^(^N){dz) ^ I <j){z)f{t,z)dz 



as — > oo for each cf) £ Cc(R''). Since 



R'' 



tJ-Tt3iN){dz) ^ J{t,z)dz = l 



R'' 
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for all i S R, we conclude that the convergence above holds for each (j> € C(,(R'^), 
which means that 

MT,3(A') ^ fit, O^' 

as iV cxD in the weak topology of probability measures, by applying Theorem 
6.8 in chapter II of [70], sometimes referred to as the "portmanteau theorem". 
■ 

The theorem above is the main result on the mean field limit in [TSJ [201 • 
1.4.4 On the choice of the initial data 

In practice, using Theorem ll.4.4l as a rigorous justification of the mean field limit 
requires being able to generate A^-tuples of the form 3(^) = {zT'ni • ■ • i ^^n n) ^ 
(R'')^ such that 

1 ^ 

satisfies 

distA/K,i(M3(W)> r") ^0 as iV ^ oo . 
Assume that /™ is a probability density on R'' such that 

\z\^f{z)dz < oo . 

Let $1 := (R'')'^ , the set of sequences of points in R^ indexed by N*. Let 
T be the cr-algebra on ft generated by cylinders, i.e. by sets of the form 

Jl Bn with Bn Borel set in R"* 

n>l 

and Bn — IV^ for all but finitely many n . 

Finally, we endow the measurable space {^,J-) with the probability measure 
p ._ (^pn^<^(X) ^ defined on the set of cylinders of il by the formula 

\n>l j n>\ 

(Notice that p"{B„) — 1 for all but finitely many n, since Bn — tV^ except for 
finitely many n.) 

Theorem 1.4.5 For each z™ = (4")fc>i e ^, let Z'^' = {zl'^, ...,z'^). Then 

distM/f.i(A^z-,/"^')^0 
as TV oo for P-a.e. z™ G fi. 
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Proof. For (j> G Cc(R'') or (l){z) = \z\, consider the sequence of random variables 
on (51, J^) defined by 

y„(z) = (j>{Zn) , 

where 

z := (zi,...,2:„,...) G n. 
The random variables Yn are identically distributed, since 



P(r„ > a) = / l^(,)>,r"(z)dz 



is independent of n. 

The random variables Yn are also independent, since for all > 1 and all 
ffi, • ■ • G CbCR), one has 



N ^ N 



EP(gi(yi) . . .gjv(r^)) = n / .g,(</>(z))r"(z)dz = [] ^^i9k{Yk)) 

k^i-^^" k=i 



Finally, the random variables Yn have finite variance since 



\z\'P'^{z)dz <OQ. 



By the strong law of large numbers (see Theorem 3.27 in [21]), one has 

N ^ , N 

4>{z)r{z)dz 



k=l I k=l 



for P-a.e. z. 

Since Cc(R'') is separable, one can assume that the P-negligible set is the 
same for all (p £ Cc(R''), and take its union with the one corresponding to 
(f>{z) = \z\. This means precisely that 



1 ^ 



k=l 

weakly in 'Pi{'R.'^) for P-a.e. z e il. One concludes the proof with the lemma 
below. ■ 

Lemma 1.4.6 The Monge-Kantorovich distance distjv/if.i metricizes the topol- 
ogy of weak convergence on Vi^R"^). In other words, given a sequence (/x„)„>i 
of elements of Vi(R'^) and /i £ 7'i(R'^), the two following statements are equiv- 
alent: 

(1) AislMK.iil^m l^) —i' as n oo; 

(2) fin fJ- weakly in V^R"^) as n ^ oo and 

sup / \z\l^z\>RfJ'n{dz) ^ OS i? — ^ oo . 
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For a proof of Lemma 11.4.61 see [5S] . 

Exercise: The reader is invited to verify the fact that one can choose the P- 
neghgible set that appears in the proof of Theorem 11.4.51 to be the same for all 
(j) G Cc(R'^) and for (j){z) = \z\. Here is an outline of the argument, 
a) Let i? > 0; let En be the set of real- valued continuous functions defined on 
[~R,R]'^ that vanish identically on d[—R,R]'^, equipped with the sup-norm 

\\<j)\\ := sup \(j}{x)\ . 

xe[~B.,R]'' 

Prove that En is a separable Banach space. 
Denote by Af^ be the set of z e such that 



1 ^ 



k=l 



does not converge to 



^{z)riz)dz 



as N oo. Let R > and let (0„)„>i be a dense sequence of elements of E^, 
extended by to R''. Define 



b) Prove that 



cp{z)r{z)dz 



as — > oo for all (f) G Eb, and all z <^ Mr. (Hint: pick G Er and e > 0, and 
choose m :— m{(j), e) such that \\(p — <j>m\\ < £• With the decomposition 



1 ^ 



fc=i 



<j>iz)nz)dz 



1 ^ 

fc=i 



1 ^ \ f 



(t>Uz)r{z)dz 



,{z)~(t>(z))r{z)dz, 



prove that 



0(z)r"(z)rfz 



< 3e 



for all z ^ Mr provided that iV > A^o = ^o(e, 0)-) 
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c) Complete the proof of Theorem 11.4.51 

Thus, using Theorem II .4.41 to prove the mean field limit requires choosing 

5{N)^{ziy,...,zJ-^^)e{IlT 

for each iV > 1 so that 

Theorem 11.4.51 provides us with a strategy for making this choice, which is 
to draw an infinite sequence at random and independently with distribution 
pnc£d^ and to set •= -z™. This strategy avoids the unpleasant task of 
having to change the first terms in 3(-^) as — t- oo. 

Since Dobrushin's estimate bounds distjvfif,! (/(i, O-^''; MTt3(Ar)) in terms of 
distMi(:,i(/*"-S?'*, M3(Ar)); having an exphcit bound on distMK,i(/*"-5?''^, M3(Ar)) 
would provide us with a quantitative error estimate for the mean field limit. 
Such a bound will be given below — see Theorem 11.6.41 

More details on the topics discussed in the present section are to be found 
in [20], as well as a precise statement concerning the behavior of fluctuations 
around the mean field limit — in some sense, the asymptotic behavior at next 
order after the mean field limit (see Theorem 3.5 in ^20j). 

1.5 The BBGKY hierarchy and the mean field 
Umit 

In the previous derivation of the mean field limit of the iV-particle system with 
interaction kernel K satisfying assumptions (HK1)-(HK2), we benefited from 
a happy circumstance, i.e. the fact that the empirical measure built on any 
solution of the A^-particle ODE system is an exact solution of the mean field 
PDE. This is why the mean field limit was reduced to the stability of the solution 
of the mean field PDE in terms of its initial data, which follows from Dobrushin's 
stability estimate. 

However, there are other situations in statistical mechanics where the em- 
pirical measure built on solutions of the TV-particle ODE system may not be an 
exact solution of the target equation — the best known example of this being 
the Boltzmann equation of the kinetic theory of gases. There are various exam- 
ples of such situations where the mean field limit can nevertheless be justified 
rigorously — see for instance 73 , and [92 in the case of random dynamics. 

There are also situations where there is no clear notion of empirical measure 
— think for instance to the A^-body problem in quantum mechanics: in that 
case, it is impossible to exactly localize any one of the N particles in phase 
space, according to the Heisenberg uncertainty principle. 

In the present section, we present another approach to the mean field limit 
of A^-particle systems, that is, in some sense, more systematic than the method 
based on empirical measure and that can be applied to a greater variety of 
situations (including quantum models, as we shall see later). 
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1.5.1 iV-pcirticle distributions 

The state at time t of a system of N identical particles located at the posi- 
tions zi{t), . . . ,ZN{t)) in the single-particle phase space R'' was described in 
the previous section by means of the empirical measure 



where Zpf{t) :— {zi(t), . . . , ZN{t)). The empirical is a probability measure 
in the single-particle phase space R'^ as mentioned above. This measure is 
parametrized by the element Z^it) of the A^-particle phase space (R'')^. 

Another way of describing the state of the same system of N particles at 
time t is to use its AT-particle distribution function, that is 

-Fjv(t, z\,..., zn) ■ 

More generally, one could think of the AT-particle distribution as being a prob- 
ability measure on the A''-body phase space (R"^)^ 



The meaning of this A''-particle distribution is as follows. Let Aj c R be Borel 

measurable sets for j — 1, . . . ,N; then, the joint probability at time t to have 
particle 1 in. Ai, particle 2 in A2. . . and particle A'' in ^jv is 



Now, we are interested in situations where all the particles in the A'-particle 
systcim considered are identical. (For example, all electrons in the universe are 
identical; ions of any given species in a plasma are identical too.) Therefore, 
for any permutation a € &n, the joint probability of having particle 1 in Ai, 
particle 2 in ^2 • • • and particle A' in An is equal to the joint probability of having 
particle 1 in A^— particle 2 in A^— ip). . . and particle N in A^-i(^n)- This 
is indeed obvious since it is impossible to distinguish particle 1 from particle 
cr~^(l), particle 2 from particle a~^{2). . . and particle N from particle a~^{N). 




1=1 



FN{t, dzi . . . dzjsr) . 




Thus 




Equivalently 




1AixA2X...xAjv {zi,Z2, . . . , ZN)F{t, dzidz2 



dZN) 




dZN) ■ 
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For a G 6j\j, define 

Scr ■ {Zl,. . ■ , Zm) ^ (Zo-(l), • ■ • , 2;cr(JV)) ■ 

The equality above is recast as 

1Ai xA2X...xA„ {zi,Z2, ZN)F{t, dzidz2 . . . dzpf) 
1AixA2X...xA„ ° Sa{zi, ZN)F{t, dzidz2 ■ . ■ dzN) 
for all Ai, ^2, . . . , Ajv Borel subsets of R'^. This is equivalent to the equality 

for all cr S ©AT. 

Obviously, when is a probability density instead of a probability measure, 
the condition 

S„#F[t, ■)^'^^ = F{t, ■)^'^^ 
for all cr e ©AT is equivalent to the condition 

F(t, 2^(1), . . . , 2:<T(Ar)) = F(t, . . . , Zat)) = F{t, Zl,..., zn) 

for all cr £ ©AT and for all zi, . . . ,zn G R''- In other words, the function 
{zi,...,zn) H' F(t,zi,...,ZN) 

is symmetric. 

1.5.2 Marginal distributions of symmetric A^-particle dis- 
tributions 

There is however one serious difficulty in considering A^-particle distributions 
in the context of the mean field limit. Indeed, this limit assumes that N — > oo, 
so that one would have to deal with "functions of infinitely many variables" in 
this limit, which does not make much sense at first sighil^^l. 

A traditional way of circumventing this difficulty is by considering the string 
of marginal distributions of the iV-particle distribution. Before giving precise 
definitions, let us explain the idea in simple geometrical terms. 

Consider a sphere centered at the origin in the 3 dimensional Euclidean 
space. The only missing information in order to completely define this sphere 
is its diameter. In other words, the sphere is completely determined as soon 
as one knows its orthogonal projection on any axis passing through the origin. 
However, if one does not know a priori that the object is a sphere, its orthogonal 



^^This last statement is not completely correct, as P.-L. Lions recently proposed a well 
defined mathematical object that would play the role of a "symmetric function of infinitely 
many variables that is slowly varying in each variable" : see 1671 and section 11.7.31 below. 
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projections on each axis passing through the origin is not enough in order to 
reconstruct completely the object, since they will not distinguish between a ball 
centered at the origin and its boundary that is a sphere of equal radius. 

The situation that we consider here is slightly more complicated, since the 
group of symmetries is not the orthogonal group, but the group generated by 
the reflections exchanging two coordinate axis in the N dimensional Euclidean 
space. Knowing the orthogonal projection of a set that is invariant under the 
action of this group on the first coordinate axis is again not suflScient as it 
will not distinguish between a sphere of radius r centered at the origin and 
the (hyper)cube of side 2r centered at the origin with edges parallel to the 
coordinate axis. Knowing the orthogonal projection on any one of the planes 
defined by two coordinate axis removes this ambiguity. 

Considering marginals of an iV-particle distribution is the analogous oper- 
ation on probability measures. Denote by Psym{{^'^)^) the set of symmetric 
probability measures on the A/'-particle phase space, i.e. 



for all fj e Sat and all zi, . . . ,zn ^ R . 

Definition 1.5.1 For each N S N*, each Pn G Vsym{i^'^)^) o.nd each k e 
{!,..., N}, the k-particle marginal of Pn is the element of'Psym{(R'^)'') defined 
by the formula 



for each test function (j) G Cb{{'R'^)'') . We shall systematically use the convention 



PN{dzi . . . dZN) = Fn{zi, . . . , ZN)dZ\ . . . dZN 

where F^ is a symmetric probability density on (R'')^, then, for each k = 
1,. . . ,N, one has 




Sct{zi,...,Zn) = {z^{l),---,Za{N)) 



I 4>{z-i,...,Zk)PN:k{dz-i...dZk)= I 4>{zu---,Zk)PN{dz-i...dzN) 



Pnm = whenever j > N . 



If 



PN:k{dZl . . . dZk) 



FN:k{zl, 



Zk)dzi . . . dzk 



with 




Obviously, F^-.k is also a symmetric probability density on (R'*)'^. 
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The following elementary exercise confirms the analogy between the orthog- 
onal projections of a subset of the Euclidean space on the subspaces generated 
by the coordinate axis and the marginal distributions associated to a symmetric 
probability on the A''-particle phase space. 

Exercise: Consider for each N e N* and each k = 1,. . . ,N the orthogonal 
projection 

: (R'^)^ 3 (21, ...,zm)^ (zi, ...,Zk)e (R'^)'= . 
Check that, for each i-V & 'Psymii^'^)^), one has 

PN:k = P^#-PiV , 

and that 

{PN:k):j = Pn-.j for all j, k such that 1 < j < k < N . 

An important example of symmetric A/'-particle distributions is the case of 
factorized distributions. Given a probability density / on R'', we consider for 
each N e N* the iV-particle probability density Fm defined by the formula 

N 

Fn{zi, . . . ,zn) --^Yi ■^(^■J') ■ 
This -/V-particle probability density is denoted as follows: 

Obviously = is a symmetric iV-particle probability distribution, and 
its marginals are also factorized distributions, since 

Fn = f^"" ^ FM:k = 

forallfc = l,...,A^. 

There is a very nice characterization of factorized distributions in terms of 
entropy. We shall not use it in the sequel. Nevertheless, it is important to know 
it, and we leave it as an exercise. 

Exercise: For each probability density / on R** and each N e N*, define 

EN{f) — {Fn symmetric probability density on (R**)^ s.t. F^-.i = /} • 
We want to prove that Fn = J®^ realizes 

inf / Fn\^Fn{zi, . . . ,ZN)dz\ . . .dzN ■ 

1) Let a = (ai, . . . , ad) G (R!J_)'' satisfy ai + . . . + = 1, and consider 
M(a) := {A = e Md(R;) s.t. (1, . . . , 1) • A = a} . 
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Find the critical points of the function 

d 

H : M{a)3 A^ ^ A^j In Aij g R . 
2) Prove that, for each x,y > 



(l){x,y) := 




-x + y> 0, 



with equahty if and only ii x — y. 

3) Express in terms of (j){Aij,aiaj) the quantity 

d 

[Aij In Aij - aiUj ln{aiaj)) . 

4) Find 

inf H(A) . 

AeM(a) 

5) Using the intuition provided by questions 1-4, solve the minimization problem 

inf / Fjy In Fn{zi, zj\i)dzi ... dzN ■ 

-FiveEjv(/) J(Rd)N 

References for this and the previous section are chapter 3 in [17] and chapter 
3 in [26]. 

1.5.3 The A^-particle Liouville equation 

We have explained above how the state of a system of N identical particles is 
described by a symmetric probability measure on the A^-particle phase space. 
Our next task is to define the evolution of such a probability measure, knowing 
that the positions of the particles in phase space are governed by the system of 
A^-particle ODEs 

r 1 ^ 

. z.(o) = zr . 

As explained in Theorem ll.3.1[ whenever K S C^(R'^xR'^, R"*) satisfies assump- 
tions (HK1-HK2), the system of ODEs above generates a flow on the A^-particle 
phase space (R'')^ denoted by Tt and defined by the formula 

Tt{zr,...,z^j^) {z,{t),...,ZN{t)) 

for all t eR. 
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Given a A''-particle symmetric probability measure F™ S Vsymii'^'^)^), we 

set 

FN{t):^Tt#FI^, ten. 

This formula defines FN{t) as the unique weak solution in C (R; w — V {{'R.'^)^ )) 
of the Cauchy problem for the iV-particle Liouville equation 



1 ^ 



Exercise: Check this, by using the method of characteristics as in the exercise 
following the statement of Theorem 11.3.11 (For the solution, see chapter 1 in 

m)- 

As explained above, it is important that the probability measure in the iV- 
particle phase space describing the state of a system of identical particles should 
be symmetric. Whether this symmetry property is propagated by the flow of the 
iV-particle Liouville equation is therefore a very natural question. The answer 
to that question is given by the following proposition. 

Proposition 1.5.2 Assume that the interaction kernel K e C^(R'* x R'',R'^) 
satisfies the assumptions (HK1)-(HK2) so that the N -particle ODE system de- 
fines a unique flow Tt on (R'')^ as proved in Theorem \1.3.1\ For each a G ©at 
we denote as above by Sa the transformation on (R"*)^ defined by 

Saizi, Zn) = (2:cr(i), . . . , Za(Af)) 

for all Zi, . . . , zn G R''. 

1) For all (7 G &n and all < G R, one has 

TtSa — Sa-Tt , 

viz. 

t I— >■ • • • , •2(T(Ar)(i)) is the solution of the N -particle ODE system 

with initial condition (•^"(i), • • ■ , z^a{N)) i 

2) For each F™ G P((R'*)^) and all t G R, the probability measure 

F^it) rt#F™ 
is symmetric if F^ is symmetric, i.e. 

S^#FN{t) ^ FN(t) for allt&H, 

if 

C 41 rpin Tpin 
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Proof. An elementary computation shows that 
satisfies the same A^-particle ODE system as 

t K> {Zl{t), . . .,ZN{t)) . 

By uniqueness of the solution of the Cauchy problem for this ODE system 
(Theorem ll.3.1|) . this is therefore the unique solution of that problem with 
initial data {z^Jf-^y ■ • ■ , z^JIn-))- In other words, 

TtSa{zl", . . . , z]^) = . . . , = SaTt{zl", . . . , Z™) 

for all (z™, . . . , z™) G (R^)^ and all t G R, which proves statement 1). 
As for statement 2), observe that 

S,#FN{t) = S,#{Tt4F'^) = {S„Tt)#F'^" 

= {TtS.)#F]? = TtMS.#F'r?) = Tt#F'^ = FN{t) 
for alH e R and all a G ©^v, which is precisely the desired relation. ■ 

Finally, we discuss the growth of F^it, zi, . . . , zjv) as |zi| + . . . + |zAr| — > oo. 

Lemma 1.5.3 Under the assumptions (HK1)-(HK2) on the interaction kernel 
K , one has 

||T,(zi",...,zi,")||i<e2^l*l||(zr,...,ziv")lli 

with the notation 

||(zi, . . . ,zjv)||i = |zi| + . . . + |zjv| . 
In particular, if F^ G 7'i((R'*)^), then Tt^F'j^'' G 7'i((R'')^), and one has 

\\{zi, . . . , ZN)\\iTt#F'^^{dzi . . .dzN) 
< e^^l*! / ||(zi,...,ZAr)||iF™(dzi...dz^) 

J(R<i)« 

for all < G R. 
Proof. Observe that 

du.. . 1 ^ ...... ^ . u.. 



,|(zi(t), . . . , zjv(t))||i = ^ E i^(z,(i), z,(i)) 

<-Y^L{\zm + \zM) 

= Y.L(\zdt)\ + -\\{z,{t),...,zNm\ 



= 2L\\{zi{t),...,ZNm\i, 
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and conclude by the Gronwall inequality. 
Moreover, if F^" G Vi{(R'^)^), then 

/ \\{zu...,ZN)\\lTt#F},^{dZi...dZN) 

= / \\Tt{z,,...,ZN)\\iF}i\dzi...dzN) 

< 6^^1*1 / \\{zi,...,ZN)\\lF}^{dZi...dZN)<^, 

SO that rt#F™ e Pi((R'^)^). ■ 
1.5.4 The BBGKY hierarchy 

The curious designation for this procedure finds its origin in the names of N.N. 
Bogoliubov, M. Born, H.S. Green, J.G. Kirkwood and J. Yvon, who introduced 
it in various contexts. 

Before presenting the BBGKY hierarchy in detail, we first discuss the main 
reason for considering it in the first place. 

As explained above, the A''-particle distribution Fn is defined on the A''- 
particle phase space (R'')''^, whose dimension increases as iV — > oo. In other 
words, the number of variables in Fjv goes to infinity with A'^, so that the 
exact nature of the limiting object associated with F^ is not entirely obvious. 
Therefore, we seek to describe the behavior of in the large N limit by 
considering instead its first marginal F^-.i in that limit. By doing so, we avoid 
the problem of having the number of variables going to infinity with A', at the 
expense of losing some amount of information by reducing Fn to Fpja. 

To fulfill this program, it would be desirable to know the evolution of F^-.i, 
typically by means of a PDE in the 1-particle phase space to be satisfied by 

Fn:1- 

Unfortunately, because of the interaction modeled by the kernel K, it is 
impossible to find a closed equation for i*V:i. Instead, the equation for F]\[:i 
deduced from the A^-particle Liouville equation involves Fn:2- Again it is impos- 
sible to find a closed equation governing the evolution of Fjv:2 when Af > 2, as 
the pairwise interaction between particles will involve F;v:3 in the equation for 
Fn:2 deduced from the TV-particle Liouville equation. By the same token, one 
can check that it is impossible to derive from the A^-particle Liouville equation 
a closed system of equations for finitely many — i.e. m < N — of the marginals 

FN:k- 

Even though the system of PDEs governing the evolution of the marginal 
distributions F]\f:k cannot be put in closed form, these equations are nevertheless 
interesting, and we explain how to derive them below. 

We start from the A'-particle Liouville equation satisfied by Fn: 

1 " 

dtFN + j^Yl div^i(^('2i' ^j)^^) = . 
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That t I— >■ Fjv(t) is a weak solution of this equation means that, for each $ = 
$(zi, . . . , zjv) in Ci((R'^)^), one has 

dt ^N{zi,..-,ZN)FN{t,dZi...dZN) 

J(R<i)» 

1 ^ /• 

in the sense of distributions on R. 

Wc first seek an equation for the first marginal i^Ar^. Let G CliTif): 
writing the weak formulation of the iV-particle Liouville equation for the test 
function ^{zi, . . . , zjq) = 0(zi), we see that 



/ (t){zi)FN{t,dzi . . .dzN) 
1 ^ f 

= T?y^. K{zi,Zj)-'\/4>{zi)FN{t,dzi...dzN)- 



(Notice that the term associated to j = 1 in the summation on the right hand 
side of the equality above vanishes identically since K vanishes on the diagonal 
by assumption (HKl).) 
In the term 

/ K{zi, Zj) ■ \/(l){zi)FN{t, dz-L... dzN) 

we exchange the variables Z2 and Zj. Denoting by (J2j & &n the transposition 
exchanging 2 and j, one has 

/ K{zi, Zj) ■ \7(l){zi)FN{t, dzi... dzN) 

= / K{zi,Z2) ■V(l){zx)Sa^ FN{t,dZx . . .dZN) 

J(R<i)'v 

= / K{zi,Z2)-V<i){zi)FN{t,dzx...dzN)- 

by symmetry of FN{t). 
Therefore 



dt / (l){zi)FN{t, dzi... dzN) 

J(R<i)« 

N -I f 

= — — / K{zi,Z2)-'V(l){zi)FN{dzi...dzN). 
The integral on the left hand side of the equality above is recast as follows: 



/ 4>{zi)FN{t,dzi . . .dzN) = / 
J(R<i)« Jr' 



(zi)FiV:i(i, dzi) . 

/(R<i)« ^R<* 
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By the same token, the integral on the right hand side becomes 

/ K{zi,Z2) ■'V(t){zi)FN{t,dzi . . .dzN) 

= / K{zi,Z2) ■'V(l){zi)FN:2{t,dZidZ2) . 

Thus, for each e C^(R''), one has 

f N -1 f 

dt 4>{zi)F]si;i{t,dzi) = — I K{zi,Z2) ■V4>{zi)F]si;2(t,dzidz2) , 

jRd -'^ J(R<')2 

which is the weak formulation of 

- 1 f 

dtFN:l-\ r^div^i / K{zi,Z2)FN:2{-,dZ2) = . 

Equivalently 

N -1 

dtFNA + -Jj— div^, [K{zi, Z2)Fn:2]:1 = , 

where K{zx,Z2)Fm:2 designates the Radon measure defined on (R'')^ as the 
linear functional on bounded continuous functions given by the formula 

{K{zi,Z2)FN:2,i') ■■= j ^l){zl,Z2)K{zi,Z2)FN:2{dZidZ2) , 

J(R<i)2 

while [K{zi,Z2)Fiq.2]:i designates the Radon measure defined on R'' by the 
formula 

{[K{zi,Z2)FN:2\:l-,(t>) ■= / (t){zi)K {zi, Z2)Fn ■.2{dzidZ2) ■ 

J(R<*)2 

(Equivalently, 

[K{zuZ2)Fn:2U ■= VI#[K{zuZ2)Fn:2] 

where we recall that 'P\ is the orthogonal projection defined by the formula 

P^(zi,Z2) := Zl.) 

In any case, as anticipated, the equation for the first marginal distribution 
Fjv:i involves the second marginal distribution Fjv:2- 

Wc next proceed to derive the equations satisfied by the sequence of marginal 
distributions F^-.j for j = 2, . . . , N; this derivation will proceed as in the case 
3 = 1, except for one additional term. 

For 1 < J < A^, wc write the weak formulation of the A^-particle Liouville 
equation with test function ^{zi, . . . , Zn) = </'(^i, • • • , Zj), where (j) S C!b{i^'^y). 
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Thus 



dt (t){zi,...,Zj)FN{t,dzi...dzN) 

J(R<i)« 





Notice that the range of the index I is hmited to {1, . . . , 7} since the test function 
(j) does not depend on the variables Zj+i, . . . , zm- The range of the index k re- 
mains {1, . . . , iV}, and we have decomposed it into {1, . . . , j} and {j + 1, . . . , N}. 
This decomposition is quite natural, as the sum involving k,l G {1, . . . , j} ac- 
counts for the pairwise interactions between the j particles whose state is de- 
scribed by FM:j^ while the sum involving I S {1, . . . , j} and k e {j + 1, . . . , TV} 
accounts for the pairwise interactions between each one of the j particles whose 
state is described by Fjv:j and the N — j other particles in the system. 
As in the case j = 1, 

/ (l){zi,...,Zj)FN{t,dzi...dZN)= / (j){zi,...,Zj)FN:j{t,dzi...dZj), 

and by the same token, if 1 < k,l < j 



lil<l<j<k<N, denote by crj^i^k G ©at the transposition exchanging 
j + 1 and k. Then 




,Zj)FN:j{t,dZi...dZj) . 



,Zj)FN{t,dzi . ..dzN) 




.dZN) 




,Zj)FN{t,dzi 



FN{t,dzi 



.dZN) 



.dZN) 



by symmetry of Fiq{t). Then 




,Zj)FN{t,dzi 



dZN) 




K{zi, Zj+i) ■ \/zi<t>{zi,. ■ ■ , Zj)FN:j+l{t, dZi... dZj+i) . 
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Finally, we obtain the equality 

ri r 

(Zi,..., Zj)FM:j{t, dzi... dZj) 



d 

di 



1 ^ f 

H r^y" / K{zi,Zj+i) - V ^^(j){zx,. . . ,Zj)FM:j+i{t,dzi . . .dzj+r) 



to be verified for each </> e Cl{{BJ^y). 

This is the weak formulation of the equation 



dtFN:j + ^^-^ V div^, I K{ZI, Zj+i)FN:j+l{-, dZj+l) 

^ 1=1 



1 ' 

- ^ div,,(K(z;,Z,)FjV:,)=0. 

Equivalently 



;=i 



1 

+ iV ^ div.,(i^(^i,^fe)f'jv:i) = 0, 

A;,;=l 

where /^(zi. Zj+i)i<Ar:j+i designates the Radon measure defined on (R'')-'+'^ as 
the linear functional on bounded continuous functions given by the formula 

lyK{zuZj+x)FN:j+X,'>l)) 

■■= / Xp{zi,...,Zj+i)K{zi,Zj+i)FN:j+l{dZi...dZj+i), 
J(R<i)j + i 

while [K{zi, Zj+i)FN:j+i]:j designates the Radon measure defined on (R'')-' by 
the formula 

{[K{zi,Zj+l)FN:j+l]:j,(j)} 

:= / (j){zi,...,Zj)K{zi,Zj+i)FN:j+i{dzi...dzj+i) . 
Equivalently, 

[K{zi,Zj+-i)FN:J+l]:j := P^j^i#[K {zi, Zj+i)FN:j+l] , 
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where we recall that 'Pj^i is the orthogonal projection defined by the formula 
P]+iizi,...,Zj+i) {zi,...,Zj). 

The equation obtained in the case j = N is nothing but the iV-particle 
Liouville equation itself since FjV:Ar = Fjv and -FV:j = for all j > N: thus 



dtFN:N + ^ div^^{K{zi,Zk)FN:N) =0. 
We summarize the above lengthy computations in the following theorem, 

yr, 



where the notation Vi^symii'R-'^)^) designates 



Theorem 1.5.4 Assume that the interaction kernel K £ C^(R'^ x R'',R'') 
satisfies (HK1-HK2). Let F}^ e Vi,sym{{R'^)^), and let Fnit) = TtH^F^^ for 
all t G R; where Tt is the flow defined on (R'')^ hy the N -particle ODE system 
as in Theorem \1.3.1\ The sequence of marginal distributions Fj^-j of F^ with 
j — 1, . . . ^ N is a weak solution of the string of equations 

TV - 1 

dtFN-.i H div^^[K{zi, Z2)Fn.,2]:1 = 0, 

dtFN-.j + ^^^^ div^, [K{zi,Zj+i)FN:j+i]:j 
1=1 



1 ^ 

- ^ dw,,{K{zi,Zk)FN..,)=0, 



j=2,...,7V-l 



k,l = l 

1 ^ 

dtFN:N + X! 'i^Vzi{Kizi,Zk)FN:N) = 0, 

fe,;=i 

and satisfies the initial conditions 

This string of equations bears the name of BBGKY hierarchy for the TV- 
particle system whose dynamics is defined by the ODE system 



1 ^ 



i=i 

In spite of all the (somewhat) technical computations involved in the deriva- 
tion of the BBGKY hierarchy, the careful reader will notice that 

a) the BBGKY hierarchy is a consequence of the A^-particle Liouville equation, 
but 
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b) it contains the A'^-particle Liouville equation (which is the last equation in 
the hierarchy). 

This observation might cast some doubts as to the interest of considering 
the BBGKY hierarchy instead of the Liouville equation itself, since both contain 
exactly the same amount of information. 

References for this section are chapters 3-4 in [26] , [88l [90] , chapter 3 in [T^ 
and [ID]. 



1.5.5 The mean field hierarchy and factorized distribu- 
tions 

In the present section, our discussion of the mean field limit becomes purely 
formal. 

Our aim is to pass to the limit in each equation in the BBGKY hierarchy 
as — >■ CO, keeping j > 1 fixed. Assume that F^.j — ^ Fj as oo (in some 

sense to be made precise) for all j > 1. Then, in the limit as N oo 



N-j 
N 

while 

so that 



Kizi,Zj+i)FN.j+i{dzj+i) ^ / K{zi,Zj+i)Fj+i{dzj+i) , 



1 

N 



— K{zi,Zk)FN:j 



dtF, + V div„ / K{zuZj+i)Fj+i{dz.j+i) = , J > 1 . 
1=1 -^R-" 

This hierarchy of equations is henceforth referred to as the mean field hierar- 
chy, or the Vlasov hierarchy. The similarities between this mean field hierarchy 
and the BBGKY hierarchy are striking. Yet there is an important difference: the 
mean field hierarchy is an infinite hierarchy of equations — unlike the BBGKY 
hierarchy, which contains only N equations, where N is the total number of 
particles. The physical meaning of this infinite hierarchy of equations will be 
explained in section H. 7. 2 1 

However, this infinite hierarchy is directly related to the mean field equation 
by the following observation. 

Proposition 1.5.5 Assume that the interaction kernel K e C^(R'' x R'',R'^) 
satisfies assumptions (HK1)-HK2). Let he a smooth (at least C^) probability 
density on R"^ such that 



\z\f"{z)dz <oo. 
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Assume that the Cauchy problem for the mean field equation 



dtfit, z) + div, (^/(t, z) j^^ K{z, z')fit, z')dz'j = , 

lt=0 ~ •' ' 

has a classical (at least of class C^) solution f = f{t, z). Set fj{t, •) = f{t, ■)^^ , 
i.e. 

3 



fj{t,zi,...,Zj) = Yl f{t,Zk) 



fe=l 



for each f G R and each zi, . . . ,Zj S R . 

Then the sequence {fj)j>i is a solution of the infinite mean field hierarchy 

dtfj{zi,...,Zj)+y2'^^^^i / Zj+i)fj+i{zi, ■ ■ ■ , Zj+i)dzj+i =0 

for all j > I. 

Proof. Since / is of class at least, one has 

j j 

dtfjit, zi,..., Zj) = X! n ^Odtfit, Zk) 

k=l ' = 1 

= - E n ^^^-^ f I ^(^fc' ^')dz'] 

= -^div,, [flf{t,zi) j K{zk,z')f{t,z')dz'\ 

= -X]*^^^^'= (^J^^Ki^k,z')fj+i{t,zi,. . . ,Zj,z')dz'^ , 

which is precisely the j-th equation in the mean field hierarchy. ■ 

This crucial observation suggests the following strategy to prove the mean 
field limit by the method of hierarchies. 

Choose factorized initial data for A/'-particle Liouville equation: given /™ a 
probability density on R"^ such that 



/ \z\f{t,z)d. 



Z < 00 . 



define 

Fj^ = for each A^- > 1 , 
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i.e. 

N 

F}^{t,z,,...,zr,) = l[r{t,Zk). 

k=l 

Let Fn be the solution of the Cauchy problem for the iV-particle Liouville 
equation 

f 1 ^ 

I dtFN + ^ I] div,,{FNK{zk,zi)) = 0, 

\ k,l = l 

{FNi^^^irr'' . 

with that initial data. 

Assume that one can prove 

a) that F^-j Fj (in some sense to be made precise) for each j > 1, where Fj 
is a solution of the infinite hierarchy, and 

b) that the Cauchy problem for the infinite hierarchy has a unique solution. 

Let / be a solution of the mean field PDE with initial data /*" . 

Since the sequence fj := f^^ (for j > 1) is a solution of the infinite mean field 
hierarchy with initial data (y™)'^J by Proposition 11.5.51 statement b) implies 
that it is the solution of the infinite hierarchy for that initial data. Therefore 

FN:j Fj = J®^ as iV oo for all j>l. 

In particular, for j — \, one finds that the solution F^^ of the Liouville 
equation satisfies 

Fn:1 ^ / as N oo . 

In other words, the first marginal of the solution of the A^-particle Liouville 
equation with factorized initial data converges to the solution of the mean field 
PDE in the large N fimit. 

This is precisely the strategy outlined by Cercignani [25j for justifying rig- 
orously the Boltzmann equation in the case of the hard sphere gas. 

Notice that, in this approach, one needs to know that the Cauchy problem 
for the mean field PDE is well-posed (i.e. that it has a unique solution for 
each initial data in some appropriate functional space). On the contrary, in the 
previous approach based on the notion of empirical measure, the existence of a 
solution of the Cauchy problem for the mean field PDE is a consequence of the 
existence for finitely many particles and of the mean field limit itself, and the 
uniqueness of that solution is a consequence of Dobrushin's estimate. 

The reader might be under the impression that proving the uniqueness of 
the solution of the Cauchy problem for the infinite mean field hierarchy is a 
matter of pure routine, since the mean field hierarchy is a linear problem, at 
variance with the mean field PDE, which is nonlinear. This is obviously wrong, 
since the uniqueness of the solution of the infinite mean field hierarchy implies 
the uniqueness of the solution of the mean field PDE. In fact, the uniqueness 
property for the infinite mean field hierarchy is a very strong property and 
proving it is by no means obvious. See |89j and section 11.7.21 for a precise 
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discussion of this point — as well as of the physical meaning of the infinite 
mean field hierarchy. 

However, the strategy described in the present section (starting from the 
iV-particle Liouville equation, deriving the BBGKY hierarchy, passing to the 
limit to arrive at an infinite hierarchy of equations and concluding with the 
uniqueness of the solution of the infinite hierarchy with given initial data) has 
been used successfully in a greater variety of problems than the mean field 
limit considered in this course. For instance, the only rigorous derivation of the 
Boltzmann equation of the kinetic theory of gases known to this date (proposed 
by Lanford [65] ) follows exactly these step^lj — see also chapter 3 of [17] and 
chapters 2-4 of [26] for an account of this fundamental result. The recent paper 
[57] extends Lanford's result to short range potentials other than hard spheres, 
and gives a more detailed presentation of the Boltzmann-Grad limit than all 
previous references, even in the hard sphere case. 

1.6 Chaotic sequences, empirical measures and 
BBGKY hierarchies 

Our discussion of BBGKY hierarchies shows the importance of the following 
property of symmetric A''-particle probability measures Fjy: 

Fn-.j — > weakly as — > oo for all j > 1 fixed. 

Of course, if is a probability density on R"*, one has 



which is in general a non-factorized symmetric probability density on (R*^) 



for all J > 1 as iV — oo. 

Thus, the property above can be verified, in the limit as oo, by se- 

quences of A'^-particle probability measures that are not factorized exactly for 
each finite > 1. 

^^For the case of the Boltzmann-Grad limit for a system of N hard spheres, the infinite 
hierarchy cannot be derived rigorously from the Liouville equation by passing to the limit in 
the sense of distributions in each equation of the (finite) BBGKY hierarchy: see the discussion 
on pp. 74-75 in 1261 . The infinite Boltzmann hierarchy is derived by a different, more subtle 
procedure that is the core of the Lanford proof — see section 4.4 in |26| . 



But if ip is another probability density on R**, defining 




k=l 



then 
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Definition 1.6.1 Let p be a probability measure on IV^. A sequence Pjv of 
symmetric N -particle probability measures on (R"*)^ for all N > 1 is said to be 
chaotic, and more precisely p- chaotic, if 

Pn:j P®' weakly m P((R'*)^') 
as iV —> oo, for all J > 1 fixed. 

The notion of chaotic sequences appeared in the context of the derivation 
of kinetic equations from particle dynamics, for the first time in |58j . Perhaps 
the reason for this terminology is that this property corresponds to asymptotic 
independence of the A^-particles in the large N limit. 

1.6.1 Chaotic sequences and empirical measures 

We begin our discussion of chaotic sequences with a characterization of chaotic 
sequences in terms of empirical measures. 

Theorem 1.6.2 Let p € P(R'*); and let Pn e Vsym{{^'^)^) for each N >l. 
Then the two following properties are equivalent: 

(a) for each J > 1 

Pn:j ^ P®' 

weakly in V{{T{.'^y) as N oo; 

(b) for each (j> £ Cb(R'') and each e > 0, 

PN{{ZNe{R.X s.t. >e})^0 
as N ^ oo, where we recall that 

1 ^ 

k=l 

with Zn = (-zi, • • • 7 zn)- 
In very informal terms, 

weakly in 7'((R'')^) as — > oo for all j > 1 if and only if 

Pn Sp 

"weakly in V{V{Il'^)y^ as — > oo. In this statement, each iV-tuples Zjy is 
identified with the corresponding empirical measure fiz^ ^'^d Pn is viewed as 
a probability measure on P(R'') that is concentrated on the set of iV-particle 
empirical measures. This identification goes back to Griinbaum [47^ and is 
discussed in section [1.7.31 below. 
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Remark: Observe that property (a) is equivalent to 
(a') Pn:1 P and Pn:2 -*p<E)p weakly in 'P(R'^) and V{{Ii'^ f ) as N ^ oo. 
In fact, as we shall see, the proof given below establishes that (a')=>(b)=> (a). 

Proof. First we prove that property (a) implies property (b). 
Applying Bienayme-Chebyshev's inequality shows that 

Pn{{Zn e (R'^)^ s.t. \{iizn -P,^)\ > e}) < ^E^- K/xz, - p, . 
Then we compute 



Observe that, by symmetry of Pn, one has 

E^"(</.(^,)</.(zfe)) 



\e^-(.^(^i)2) ifj = fc. 



Thus 



hvE'^(^^) =^ Ee^"(<^(^^)'^(^^)) 



^{p^\r') = {p,<i>y 



as A'' ^ oo by (a), while 



E 



/ 1 ^ \ 1 ^ 



1 ^ 



^ • 1 



as N ^ oo, again by (a). 
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Therefore, property (a) implies that 

E^- Kmz. - P, </')P ^ {P, cl>f + (p, <l>f - 2{p, cpf = 
as A'' ^ 00, and therefore 

Pn{{Zn e (R'^)'^ s.t. Kmz, -p,^)\ > e}) < ^ -0(1) 

by Bienayme-Chebyshev's inequahty, which is precisely property (b). 

Next we prove that, conversely, property (b) implies property (a). 

Step 1: let us prove first that property (b) implies that Pn-.i P weakly as 
N ^ CO. 

Let (j) e C6(R'*). Denote 

U%{cj,) {Zn e {'R.T s.t. -p,4>)\>^}, 
V^{ct>) :=(RT\C/^(0)- 

Then 

|E^"(^z„,0) - = \^EP-{^lz. -P,4>)\ 

< E^« ( I (mz« - P, </>) I Ic/^ (^) ) + E^« ( I (mz« - P, </>> I ly^ (^) ) . 

Obviously 

E^"(Kmz. -p,<^)|lf/^(<^)) < 2\\<I>\\l^Pm{W^{^)) , 

and 

By (b), there exists N{e,(j)) such that 

N>N{e,cl>)^Pj,{UI,{(l>))<e, 

so that 

|E^-(/xz„,.^) - (p,<^)| < (2||.^|Uoo +l)e. 
On the other hand 

E^«((MZ.,<^))=E^« (^^E'^(^^)j =E^«(</.(^i)) = (^'iV:l,<A), 

so that 

AT > N{e,<i>) ^ \{Pn:1,^) - (P,<^)| < (2||.^||loo + l)e. 
which concludes step 1. 

Step 2: next we prove that property (b) implies that Pn-.j P^-' weakly as 
iV DO, for all j > 1. 



50 



CHAPTER 1. MEAN FIELD LIMIT 



Let (j) G CbiW^) \ {0}. Let Sj^ := {1, . . . , iV}{i- 'J> (the set of maps from 
{1, . . . , j} to {1, . . . , A^}) and let A'j^ be the set of one-to-one elements of Sj^. 
Then 

= E E^"(0(z,(i))...^(z,(,))) 
+ 7^ E E^"W^«(i))--.'/'(^.0)))- 

For s G one has 

E^-(0(Z,(1)) . . .0(Z,(,))) = E^"(0(zi) . . .^(z,)) = (P^:„0®^) 

by symmetry of Pjv, while, for all s & 

Now, for all j fixed 

r = Ar(Ar - l) . . . (AT - j -f 1) - A^^' as ^ oo , 
1 #{£'n\K) = -N{N-l)...{N-j + l) = o{N^) , 



so that 



Introduce 



NO 



• </.) = {Z^ G (R'')^ s.t. KmI^ - p^^,cj>^^)\ > e} , 

, Y^{j,^) = {R''f\X%{j,^). 



Observe that 



fe=l 

so that 



Therefore, property (b) implies that 
as AT ^ oo for all j > 1, all e > and all (/> G Cb{R'^) \ {0}. 
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In particular, for each j > 1, there exists Nj{e, (f) > j such that 
iV>iV,(e,0)^PAr(X^(j») <e. 

Thus 

<2||0||i^P^(X^(j,0))+e 

so that 

W > iV,(e,<^) ^ |E^«(a*|^ I < (2110111^ + 1)6. 

Therefore, for each e > and each j > 1, and for aU N > Nj{e, (p), one has 



< {2U\\i^ + l)e + 2^^ E^UAUi 



and since - = o{N^) as iV ^ oo. 



hm (p«^-PjV:,,0^^) <(2|10||i., + l)e. 
Since this holds for each e > 0, wc conclude that 

{Pn:j , 0®^) ^ (p®^' , (j}'^^) as AT ->c5o . 

This property holds for each (j> e Cfc(R'') \ {0} so that, by a classical density 
argument, we conclude that the sequence Pn is p-chaotic. ■ 

Exercise: The purpose of this exercise is to complete the "classical density ar- 
gument" used at the end of the proof of Theorem ll.6.21 For each $ e C((R'')-'), 
denote by 5*^$ the element of C((R'')^) defined by the formula 

Finally, for each (f>i, . . . , £ C(R'*), we denote hy 4>i ® . . . (pj the function 
01 (g) . . . (X) 0j : {zi, . . . ,Zj) h-> 01 (zi) . . . (I)j{zj) , 



and by 0*^^ the j-fold tensor product of with itself, 
a) Prove that 

1 

S',(0i ® . . . $5 0,) = -— ^(^101 + • . • + ^,0,)' 
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b) Prove that, for each Q S ^^^^^((R'')-') and each $ e Cf,((R'')^ ), one has 

c) Leilil <i> e Co((R'')-') and e > 0. Prove that there exists ai, . . . , a„ G R and 
a jn-tuple {(f>i.m)i<i<j.i<m<n of elements of Co(R'') such that 



< e 



where 



1*11= sup |*(zi,...,Zj)|. 

zi,...,2,-eR'' 



(Hint: apply the Stone- Weierstrass theorem.) 

d) Let Q e Vsymmy) and {Qn) N>i be a sequence of elements of sym 
such that 

(gAr>^^'>->(Q,<^®-'') 

sls N —>■ oo for each € Co(R'^). Prove that 

(Qw,<i>)^(Q,<i>) 

as ^ oo for each $ e Co{{R'^y). 

e) Conclude that 

(Qiv,*) -> (Q,$> 

as A'' — )• 00 for each $ e Cb((R'')^). (Hint: apply Theorem 6.8 in chapter H of 

m-) 

References for this section are chapter 1, section 2 in [92], and chapter 4, 
section 6 in 1261. 



1.6.2 Prom Dobrushin's theorem to the BBGKY hierar- 
chy 

In the proof of the mean-field limit based on the notion of empirical measure, 
one proves the stability of the limiting, mean field PDE in the weak topology 
of probability measures on the single-particle phase space. The method based 
on the BBGKY hierarchy involves the A^-particle phase space, and the need for 
considering marginals of the A^-particle distribution in the limit as iV — >■ oo. 
Perhaps the best way to understand the relation between these two approaches 
of the same problem is to think of the nature of the A^-particle empirical measure. 
Indeed 

Aiz„(dz) 

^''For each locally compact topological space X and each finite dimensional vector space 
E on R, we denote by Co{X, E) the set of continuous functions on X with values in E that 
converge to at infinity. We set C(i{X) := Co(X, E). 
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is a symmetric function of the N variables Z]\[ = [zi, . . . , zat), defined on the 
N -particle phase space (R'')^, and with values in the set oi probability measures 
on the single-particle phase space R'' — in the notation above, z is the variable 
in the single-particle phase space. 

Already the characterization of chaotic sequences in terms of empirical mea- 
sures obtained in the previous section clarifies the respective roles of the single 
particle phase space and of the A^-particle phase space in this limit. The A^- 
particle symmetric distribution function can be viewed as a probability on 
7'(R'^) concentrated on the set of empirical measures — in other words, the 
A^-tuple Zjv is, up to permutations of its N components, identified with the 
empirical measure ■ (This idea can be found in |47j and will be discussed 
in more detail in section [1.7.31 below. ) 

Although this is a static picture, it provides the right point of view in order 
to unify the two approaches of the mean field limit presented above, i.e. the 
approach based on empirical measures and Dobrushin's estimate, and the one 
based on the BBGKY hierarchy. More precisely, we shall prove that the sequence 
F]\[{t) of solutions of the A^-particle Liouville equation with factorized initial 
data is /(t).if''-chaotic for each t g R, where / is the solution of the mean field 
PDE. 

The proof of this result goes as follows. 

Since is factorized, it is of course chaotic (factorized probability measures 
being the first example of chaotic measures). But in fact, F}^ satisfies a stronger 
estimate than property (b) in Theorem ll.6.21 in terms of the Monge-Kantorovich 
distance distMX,2 (see Theorem 11.6.41 below). This estimate is propagated for 
all t S R by Dobrushin's inequality, involving the weaker distance distMif,i- 
The resulting bound implies that property (b) in Theorem 11.6 . 21 holds for fAr(t) 
for all t with p f{t)^'^. Applying Theorem ll.6.21 we conclude that Fff{t) is 
/(^^''-chaotic. 

Our main result in this section is summarized in the following statement. 

Theorem 1.6.3 Assume that the interaction kernel K G C^(R'' x R'',R'' sat- 
isfies assumptions (HK1)-(HK2) . Let /*" be a probability density on R'^ such 
that 

f \z\'^+^ f{z)dz <oo. 

Let Fj^ = ( jm_2"^)®^, and let Fn (t) = Tt #F}:^ be the solution of the N -particle 
Liouville equation 

( I N 

I dtFN + i;7j2 div.,(^^(^», Zj)Fn) = , 
Then, for each j > I 

FN-At)^{f(t,-)^''r^ 
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weakly in ■p((R'')^) as N ^ oo, where the probability density f{t, •) is the solu- 
tion of the mean field PDE 

dtfit, z) + div, (^/(i, z) J^^ K{z, z')f{t, z')dz'^ = , 

, -I \t=Q ~ J 

The rather stringent moment condition on the initial single-particle density 
comes from the following important result from statistics, that can be viewed 
as a quantitative variant of the law of large numbers. 

Theorem 1.6.4 (Horowitz-Karandikar) For all p G 'P(R'^) such that 

a:=(p,|z|'^+5) 

one has 

distMK,2(MZ«,p)V^(d^A^) < 



where C{a, d) is a positive constant that depends only on a and the space di- 
mension d. 

References for this result are the original article [57], and the monograph 



Taking this estimate for granted, we give the proof of Theorem ll.6.31 
Proof of Theorem 11.6.31 By the Cauchy-Schwarz inequality, for all /x, G 
7^2 (R'^) and all tt G n(/i, i/), one has 

1/2 

\x — y\TT{dxdy) < [ I I \x — y\'^TT{dxdy) 

Taking the infimum of both sides of the inequality above as tt runs through 
n(/i, v) shows that 

distMif,i(/i, i^) < diatMK.2{f^, v) ■ 

Hence 



/ 

Since Fjv(t) = Tt#F'^ , one has 



FNit){{ZN e {RT s.t. Kfiz. - > e}) 

e (R-^)^ s.t. I (mt,z"^ - f{t).^',4>)\ > e}) . 
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{(t>{x) - (j){y))n{dxdy) 



< Lip(0) / / 



X — y\TT{dxdy) . 



Taking the inf of both sides of this inequahty as tt runs through the set of 
couphngs U(firp zir^,f{t).S!f'^) shows that 

By Proposition II .4.21 and Dobrushin's inequahty (Theorem II. 4. 3p . 

\{^^nz-^ - < Lip(0)distMKaKzj,",/W=S^") 

< Lip(</>)e2^l*l distMif,i(Mzj?: Z"-^'') • 

Therefore 

FNimZ^ e {KT s.t. - /W^^</>)| > e}) 

< F'nUZ'n e (R'^)'^ s.t. distMK,i(Mzj,",r"^') > e-2i|*le/Lip(0)}) . 

Applying the Horowitz-Karandikar theorem recahed above and the Bien- 
ayme-Chebyshev inequahty shows that 

F^^{{Z^;} e (R"^)^ s.t. distMKalA'z^",/™^') > e-'^l*le/Lip(0)}) 

e^^l*! Lip((/.)2 C(a,d)' 

- £2 iV2/('i+4) ■ 

Hence 

as iV ^ oo for each e L°° n Lip(R''). By density of L°° n Lip(R'^) in Cbi^"^) 
and the imphcation (b) (a) in Theorem ll.6.21 we conclude that 

Fj,..,{t)^{f{t)J^T' 

weakly in 'P((R'')^) as TV ^- cx) for aU j > 1. 
In particular 

FN:l{t) ^ /W^'' 

weakly in P(R'') as — > oo, where f{t, •) is the solution of the mean field PDE, 
and this concludes the proof. ■ 
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1.7 Further results on mean field limits in clas- 
sical mechanics 

1.7.1 Propagation of chaos and quantitative estimates 

Following carefully our proof in the previous section, the convergence result 
stated in Theorem ll .6.31 could obviously have been formulated as a quantitative 
estimate of the form 

distMKAPN:j{t), {f{t)J^T') < , 

where ej{N) — for each j > 1 in the limit as iV ^ oo. Such an estimate would 
obviously imply Theorem 11.6.31 since the Monge-Kantorovich distance distA/i^^i 
metricizes the weak topology on 7'i(R") for all n > 1 — see Theorem 6.9 in 
[55] . although what is needed here is a consequence of Proposition ll.4.21 

In fact, there is another approach of Theorem 11.6.31 which is perhaps worth 
mentioning, since it provides additional information on the relation between 
the approach with the empirical measure and the approach based on BBGKY 
hierarchies. 

Specifically, one can prove that, if e Vi^symii'R-'^)^) and if t F^it) is 
the solution of the Cauchy problem for the iV-particle Liouville equation with 
initial data F^, assuming that the interaction kernel K satisfies the assumptions 
(HK1-HK2), then 



where Tt is the flow defined by the A^-particle ODE system, while RN,m{t) is a 
positive Radon measure on (R'i)™ whose total mass satisfies 

, „ , N , mim — 1) 

(i?jV,™W,l> -1-77^ TT^< 



{N - m)!7V'" - 2N 

This explicit formula, which can be found in [44 expresses the m-th marginal 
of the iV-particle distribution i^Ar(i) in terms of the empirical measure of the 
iV-particle system at time t, i.e. fJ-TtZ]^ j np to an error that vanishes as iV — >■ oo. 

With this formula, one easily arrives at the following quantitative estimate 
for the propagation of chaos in the mean field problem. Assume that the initial 
data is factorized, i.e. 

pin (J*" ^d-\^®N 



and that 



< oo. 



Then 



C{a,d) ' 
Ari/(d+4) 
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for alH G R and all iV > m > 1. In particular, for m = 1, one has 

distMKAPN:i{t)J{t)^'') < C7(a,d)e2^l*l/iVi/(<i+4) 

This estimate can be found in |44) : it is established independently and in a 
slightly different manner in |73) . 

In fact, the idea of obtaining quantitative estimates for the propagation 
of chaos in various situations of non equilibrium classical statistical mechanics 
has been systematically pursued in the reference [73], which bears on a more 
general class of models than the one considered so far in the present notes. 
In particular, the discussion in |73^ applies to situations where the empirical 
measure of the A'^-particle system is not an exact weak solution of the limiting 
mean field equation. 



1.7.2 Infinite iiierarchies and statistical solutions 

While the physical content of the BBGKY hierarchy is transparent (it is the 
string of equations satisfied by the marginals of the A^-particle distribution) , the 
physical meaning of the infinite mean field hierarchy is somewhat less obvious. 
We discuss this point in the present section. 

Consider the Cauchy problem for the mean field PDE: 

at/ + div,(//C/) = 0, 

J \t=0 ~ J ' 

where we recall that 

/C/(t,z) := / K{z,z')f{t,dz'), 

and where the interaction kernel K £ C^(R'^ x R'',R'^) is assumed to satisfy 
conditions (HK1-HK2). 

Denote by Gt : 7'i(R'*) i-> 7'i(R'^) the nonlinear 1-parameter group defined 

by 

Gtr 

where / is the unique solution of the Cauchy problem above. 

The notion of statistical solutions of the mean field PDE is defined by anal- 
ogy with the case of an ODE. Consider the Cauchy problem for an ODE with 
unknown 1 1-^ x{t) £ R" in the form 

x{t) = v{x{t)) , 
x(0) = xo ■ 

Assuming that v G Lip(R",R"), the Cauchy-Lipschitz theorem provides the 
existence of a global solution flow X : Rx R" — >• R", so that the map 



1 1-^ X{t,xo) 
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is the solution of the ODE above satisfying X{0,xo) = xq. The flow X corre- 
sponds with a completely deterministic notion of solution: knowing the initial 
condition xq exactly determines the solution 1 1— > X{t, xq) for all times. Suppose 
that, instead of knowing exactly the initial condition xq, one is given a proba- 
bility distribution /ip on the set R" of initial data xq. In other words, /io can 
be regarded as a "statistical initial condition" for the Cauchy problem above. 
With the flow X and the statistical initial condition /iq, one defines 

fi{t):^X{t,-)#fio, ten. 

Applying the method of characteristics (see the exercise following Theorem 
ll.3.ip shows that 

dt^i{t) + diY{p{t)v) = , 
fi{0) = /io . 

In other words, we recover the well-known fact that solutions of the transport 
equation (a 1st order PDE) can be viewed as statistical solutions of the ODE 
defined by the characteristic field of the transport operator. The idea of consid- 
ering the time-dependent probability /Lt(t) (or distribution function, when ij,{t) 
is absolutely continuous with respect to the Lebesgue measure on the single- 
particle phase space) instead of the deterministic solution X{t,xo) for each 
initial data xq lies at the core of the kinetic theory of gases. 

Let us return to the problem of defining a notion of statistical solution of 
the mean field PDE recalled above. First we need a probability measure vq on 
the space of initial conditions /*", in this case on 7^i(R'^): this probability i/q 
will be the statistical initial condition for the mean field PDE. By analogy with 
the case of the simple ODE presented above, we define 

„{t) := Gt^i^o , ten. 

The next question is to find the analogue of the transport equation satisfied 
by /i(0- 

First we need to have a better grasp on ^{t). One way to understand the 
formula above expressing i>{t) as the push- forward of vq under the map Gt 
defined by the mean field evolution is to write 

for some appropriate class of continuous (in some sense to be defined) functions 
T on 'Pi(R.'^). Certainly this class of functions should contain "polynomials" on 
7'i(R'^), i.e. linear combinations of "monomials". A monomial of degree k on 
'P(R'^) is a expression of the form 

Mk,^{p) := / (l>{zi, . . . , Zk)p{dxi) . . .p{dxk) = <j)) 

J(R'i)'= 

where cj) E CbU'R'^)'') — without loss of generality one can assume that (p is 
symmetric in the variables zi, . . . , z^. 
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Specializing the formula above to the case J" = Mj^^ as (j) runs through 
Ci){{'R'^y ) results in the equality 

/ p'^^u{t,dp)= [ {Gtp)^^Mdp). 

JViiR'') JVii'R'') 

Defining 

Fj{t) := I {Gtpr^v^idp), i>l, 

we claim that the sequence {Fj)j>i is a solution of the infinite mean field hierar- 
chy. Indeed, for each initial single particle probability distribution p G 7'i(R''), 
the sequence {{Gtp)^-')j>i is a solution of the infinite mean field hierarchy, 
which is a sequence of linear equations. Therefore (-Fj)j>i is also a solution of 
the infinite mean field hierarchy, being an average under of solutions of this 
hierarchy. 

There is another formulation of this observation. Set :— (R'')'^ , equipped 
with its product topology and the associated Borcl (T-algebra B{Q). For each 
p e Vi(R'^), we denote by the unique Borel probability measure defined 
on ^l by the formula 

=l[p{E,) 

\k>l J k>l 

for each sequence {Ek)k>i of Borel subsets of R** such that E). = R'' for all but 
finitely many fcs. Notice that 

Define 

F{t):= f p^°°u{t,dp)= [ (GtPr'^Mdp); 
then, for each j >1 and each f e R, 

so that {F{t);j)j>i is a solution of the infinite mean field hierarchy. Let us write 
the jth equation of the infinite mean field hierarchy in terms of F{t): 

j 

9tF(t):,- +^div,.(F(i):,+iif(zi,z,+i)):,- =0. 

i=l 

This equality is equivalent to the following weak formulation: for each test 
function e Cl{(Ry) 

j 

i=l 

= /F(f),^i^(z„z,+i)-V,.Vj- 
\ i=i 
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This suggests the following definition. We denote by Vsymi^) the set of 
Borcl probability measures /i on $7 such that 

for each > 1 and each a E &n, where 

Uaizi,Z2, ...):— (2cr(l), ■ • ■ , Zo-(Ar)i ^W+l: Zn+2, ■ ■ ■) ■ 

Definition 1.7.1 (Spohn [89j) Let I C H be an interval. A map 

l3t^P{t) eVsyrni^) 

is a measure-valued solution of the mean field hierarchy if and only if, for each 
j e N* and eacl^ ipj G C^U'R'^y), the map t (P(t),'0j) is of class on 
the interval I and 

9t(P(i),V,) = dt{P{t)..,,^,) = /p(i),^if(z„z,+i) . V,.^A 

for all t £ I. 

With this definition, the map II 3 t i-^ F{t) E Vgymi^) defined above 
is a measure-valued solution of the mean-field hierarchy satisfying the initial 
condition 

F(0)= / p'^°^Mdp)- 

JViiR'') 

This observation raises two natural questions: 

a) are measure-valued solutions of the mean-field hierarchy uniquely deter- 
mined by their initial data? 

b) are all the measure- valued solutions of the mean-field hierarchy defined 
by statistical solutions of the mean field PDE? 

A first useful tool in answering these questions is the following result. 

Theorem 1.7.2 (Hewitt- Savage [52]) For each P G Pgymi^), there exists 
a unique Borel probability measure tt on V{Il.'^) such that 

P = / p^°^TT{dp) . 

^^For each n,k > 1 and each finite dimensional vector space E on R, we denote by 
Cq(R" , E) the set of functions of class C*^ defined on with values in E all of whose 
partial derivatives converge to at infinity. In other words, 

C^(K",E) := {/ e C'=(R",£;) s.t. 8°' fix) -s> as -s> oo for each a e N"} . 

We denote C|5(R") := C|5^{R",R). 



1 . 7. FURTHER RESULTS ON MEAN FIELD LIMITS 



61 



The Hewitt-Savage theorem is often quoted and used as follows. Given 
iPj)j>i a sequence of probability measures such that Pj £ 7'sj^m((R'')^) satis- 
fying the compatibility condition 

(CC) Pk:j = Pj for each fc > j > 1 , 

by the Kolmogorov extension theorem |91) . there exists a unique probability 
measure P G Pgymi^) such that 

P:j = Pj for each j G N* . 

Therefore, by the Hewitt-Savage theorem, there exists a unique Borel probability 
measure tt on V{'Rf') such that 

P, = / p^^nidp) . 

JviR'') 

With this, we can answer the questions a) and b) above. 

First we consider the problem of admissible initial data. In view of the dis- 
cussion above, the initial data for the mean field hierarchy can be equivalently 
either a sequence (-F™)j>i such that Fj" S Vi.symii'^'^y) for j > 1 and satis- 
fying the compatibility conditions (CC), or a unique element F'" € T^i^) such 
that 

F™=i^, j>l, 
or a unique Borel probability measure v^"" on ViiR'^) such that 

f p^'v"^{dp), j->l, 

J-Pi(R<*) 

or equivalently 

F™ = / 

The same reasoning applies for each instant of time t ^ 0. A solution at 
time t of the infinite mean field hierarchy is a sequence {Pj{t))j>i that satisfies 
in particular Pj{t) € 'Pi.saT?i((R-'^)"') for i > 1 together with the compatibility 
condition (CC). Equivalently, this defines a unique probability measure P(t) G 
Vaymi^) such that 'P{t):j = Pj{t) for each j > 1, or a unique Borel probability 
measure 7r(<) on V^Bf-) such that 

P,{t)^ [ p®'ir{t,dp), j>l, 

or equivalently 

V{t)= I p^^7r{t,dp). 

The following important result was proved by H. Spohn |89) . It answers 
questions a) and b) above. 
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Theorem 1.7.3 Under the assumptions (HK1-HK2) on the interaction kernel 
K G C^{Hf' X R'',R''), for each Borel probability measure on Vi^R!^), the 
only measure-valued solution of the infinite mean field hierarchy with initial data 

is 

Y{t)= I {Gtp)^°°v'"{dp), ieR. 

J-Pi(R'') 

In other words, Spohn's theorem proves that the only measure valued solu- 
tion of the infinite mean field hierarchy is the statistical solution of the mean 
field PDE 

where v^" is the initial probability measure on the space Vi (R'') of initial data 
for the mean field PDE such that the initial condition for the mean field hier- 
archy is the element of Vgymi^) given by 

or equivalently the sequence 

= / e Vsymmy) , J>1. 

Notice that, in the case where z/*" = dfin, one has 

pin ^ (jm)0oo ^ equivalently Fj" = for all j > 1 , 

and 

F{t) = (Gtf'")^°° , or equivalently Fj{t) = (Gt/*")®^' for all j>l. 

In other words the statistical solution u{t) of the mean field PDE at time t is 
given by the formula 

Hi) = hit) = ^GtPn = Gt#Sfin . 

That the statistical solution u{t) is a Dirac measure for each t G R, knowing 
that i^(0) = i^*" is a Dirac measure, is equivalent to the propagation of chaos in 
the context of the mean field limit. 

Therefore, Spohn's theorem contains as a particular case the uniqueness 
theorem with factorized initial data which justifies the validity of the mean field 
limit in the approach with the BBGKY hierarchy. 

But more generally, Spohn's theorem shows that solutions of the infinite 
mean field hierarchy coincide with the notion of statistical solutions of the mean 
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field PDE. This important piece of information clarifies the meaning of the 
infinite hierarchy. 

Spohn's uniqueness theorem can also be combined with the following uniform 
stability result on the BBGKY hierarchy to produce a quantitative stability es- 
timate on statistical solutions of the mean field PDE. The following quantitative 
stability estimate uses the following variant of Monge-Kantorovich distance. For 
each P e Vi.symii'R'^)'^^) and Q G 7'i,sy™((R^)^), consider 

DiatMK,i{P,Q) = inf // distMK,i(/^^M ' Mi'jv)p('^-'^M, rf^V) ■ 

pen(P.Q) J JudM y-fidN 

In this formula, P and Q are viewed as Borel probability measures on 7'i(R'^) 
concentrated on the set of M- and TV-particle empirical measures respectively. 

Theorem 1.7.4 Let M,N > 1, and let PI"} e Vi,sy,n{{'R'^V'') and Q^' e 
Vi^symii'R-'^)^)- Assume that the interaction kernel K satisfies the conditions 
(HK1-HK2). Let 1 1— Pnit) and 1 1— >■ QN{t) be respectively the solutions of the 
M -particle and the N -particle Liouville equations with initial data P^ and Q^^ 
respectively. Then 

a) for each t G R, one has 

^l^tMK.l{PM{t),QN{t)) < e2^l*lDistMif,l(^M-Q")- 

b) for each i G R, each m,M,N G N* such that M,N > m, and for each 
bounded and Lipschitz continuous function <j)„i defined on (R'')™, one has 

\{PM:rnit) - QN:m{t), (l)m)\ 

<mfe2^l*lLip(0„)DistMK,i(^'M,Qw) + ("^-l)ll'^mllL ^ 



N 

For a proof of this result, see [44] . 

Notice however that Spohn's uniqueness theorem for the infinite mean field 
hierarchy, even in the particular case of factorized initial data, is more than 
what is needed to justify the mean field limit. It would be enough to prove 
the uniqueness of those solutions of the infinite hierarchy that are limits of 
the sequence of marginals of A^-particle distributions as iV — >■ oo. This weaker 
uniqueness property follows from Theorem ll.7.4| without using Spohn's unique- 
ness theorem. 



1.7.3 Symmetric functions of infinitely many variables 

In various places in this notes — and especially in the last theorem — we 
encountered the idea of viewing elements of Vsymii'R-'^)'^) as Borel probability 
measures on 7'(R'') concentrated on the set of empirical measures. 
In fact, the identification 

1 ^ 

{RT/6n 3 (zi, . . . , z^) ^ liz. '-^N^ ^ P.ymiR') 

k=l 
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can be found in [l^ (see especially p. 330 there). 

This point of view was pushed much further by P.-L. Lions. He constructed a 
complete mathematical framework for handling continuous symmetric functions 
of infinitely many variables that depend weakly on each variable, and used it 
systematically in his 2007-2008 lectures at the College de France on mean field 
games [67]. This remarkable circle of ideas also appears in the recent work 
of Mischler, Mouhot and Wennberg (see [72] and [73]) on the mean field limit 
and on quantitative estimates on the propagation of chaos in classical statistical 
mechanics. 

We introduce, for each X. Y £ , the notation 

dLp{X, Y) inf{e > s.t. #{k = 1, . . . , N \\xk - yk\ > e} < Ne} . 

This quantity is related to the Levy-Prokhorov distance dist^p on Borel prob- 
ability measures on Q in the following manner: 

distLp{^J.x,^J■Y) = inf dLp(X,Y„) , 

where Y^ := (?/cr(i); • ■ • i y<j{N))- We recall the definition of the Levy-Prokhorov 
distance dist^p on V{Q): 

distLp(Pi, P2) := inf i e > I inf / / l|^_j^|>,7r(dxdy) < e I . 

We also recall that the Levy-Prokhorov distance metricizes the weak topology 
on P{Q), so that {V {Q),dist lp) is a compact metric space — see for instance 

Lions' key observation is summarized in the following theorem. 

Theorem 1.7.5 (P.-L. Lions [67J) Let Q be a compact metric space and con- 
sider, for each N > 1, a symmetric function un S C{Q^). Assume that 

sup sup \un{zit ■ ■ ,zm)\ < 00 , 
N>i (zi,...,2jv)eQ" 

and that 

s\xp{\un{X)~un{Y)\ s.t. X,Y eQ^ anddLp{X,Y) < e} ^ 

uniformly in N > 1 as e 0. Then there exists U G C{V{Q)) and a subsequence 
of UN such that 

sup |ujVfc(^) — U{^z) \ as Nk 00 . 
zeQ"k 

This point of view is obviously dual of the Hewitt-Savage theorem quoted 
above. In fact, as noticed by P.-L. Lions, it leads to a very simple proof of the 
Hewitt-Savage theorem. Lions' argument [67] is sketched below. 
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Consider a sequence iPj)j>i of probability measures satisfying the compat- 
ibility condition (CC) above. Lions' idea is to consider 

CiV{Q)) 3U^ L{U) := lim / C/(/iz„)Pjv(d^Jv) e R. 

That the limit on the right hand side of the equality above exists follows from 
considering the case of a monomial. Indeed, when U — M^^^ with the notation 
used in the previous section, i.e. 



U{fi) ^ Mk,4,{fJ.) :^ / 4>{xi,...,Xk)tJ.'^ {dxi...dxk) 
where (j) € C{Q^), a straightforward computation shows that 

Mk,4f^z,,)PN{dZN) = f dp{zi,...,Zk)Pk{dzi...dzk) + 0{l/N). 



Hence 

L{Mk,4,) ^ / (f>{zi, . . . ,Zk)Pk{dzi . . .dzk) 
Jq'' 

for all A; > 1 and all (j) e C{Q''). 

The set of polynomials, i.e. of linear combinations of monomials, is a subal- 
gebra of C{V{Q)) since 

Mk,^{p)Mi^^{^i) = Mk+i, 4,(^4, in) 

with 

0® ^(a^i, ■ • ■,Xk+i) := ^(xi, . . . ,Xk)ipixk+i, ■ ■ ■,Xk+i) ■ 
We also use the convention 

Mo,i = 1 . 

This subalgebra is dense in C{V{Q)) for the topology of uniform convergence 
by the Stone- Weierstrass theorem. 
Since 

/ U{iiz^)PN{dZN) < sup \U{fi)\, 
Jq^ i-iev(Q) 

for eacllll U G C{V{Q)), and since the limit 



L{Mk.4,) := lim / Mk,4,{^JizN)PN{dZM) = I (j^izi, . . . , Zk)Pk{dzi . . . dzk) 

exists for each fc > 1 and each e C{Q^), we conclude that this limit exists for 
each U 9 C{r{Q)). 



^^Since {^{Q), dist^p) is compact, any element U of C(P(Q)) is bounded on 'P{Q), so that 

sup I U (fj.) I < oo . 
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Since obviously L > (being the Umit of Unear functionals defined by Pn 
which is a positive measure), and since 

L(l) = / P^idZN) = 1 for each N > 1 

we conclude that the linear functional L is represented by a unique probability 
measure tt e V{V{Q)), i.e. 



L{U)= [ UifxMdfj,). 

JV{Q) 

Specializing to monomials 

L{Mk,4>)= / ^{zi,...,Zk)Pk{dzi...dzk)= / Mh,4,{iJL)-K{dn) 
Jq>' Jv(q) 



/ (/ 

Jv(Q) \JQ 



(Zl, . . . , Zfc)/!* '{dzi . . . dZk) TT{dfl) 



= / (l){zi,...,Zk) { / n®'^{dzi...dzk)'K{dn) 
and this means that 

^I'^'^nid^i) 



with the unique probability measure tt defined above. This is precisely the 
representation formula in the Hewitt-Savage theorem. 



1.7.4 The case of singular interaction kernels 

The method for proving the mean field limit presented above is based on Do- 
brushin's estimate and, as such, is limited to cases where the interaction kernel 
k is Lipschitz continuous in both its arguments. This is most annoying since 
many interaction kernels of physical interest are singular on the diagonal. All 
the examples presented in the first section (i.e. the Vlasov-Poisson system and 
the vorticity formulation of the Euler equation for incompressible fluids in two 
space dimensions) involve the fundamental solution of the Laplacian, leading to 
interaction kernels that become singular as the distance between the two inter- 
acting particles vanishes. Obtaining rigorous derivations of both models as the 
mean field limit of large particle systems remains an important open problem. 

However, some remarkable results have been obtained in this direction. 

In the case of the vorticity formulation of the two-dimensional Euler equation 
for incompressible fluids, the mean field limit for the dynamics of a large mimbcr 
of vortex centers is analogous to the convergence problem for vortex methods in 
the numerical analysis of the Euler equation. These methods approximate the 
vorticity field lo = (jj{t, x) as follows: 

N 
fe=l 
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where uJk is the (constant) intensity of the vortex centered at Xk{t). (Notice the 
sUght difference with the mean-field hmit discussed above, where each vortex 
would have the same intensity iOk = 1/N.) The motion of the vortex centers 
Xk (t) is governed by the following ODE system: 

N 

i;fc(0 =^w;i^e(a;fe(<) -x,(i)), k^l,...,N, 

1=1 

where K^^ is an approximation of the vortex interaction kernel 

Jx 

(with J standing for the rotation of an angle — §)• 

In one variant of these methods, called "the vortex blob method" , the inter- 
action potential is truncated near the singularity at the origin so as to remain 
smooth as the distance between interacting vortices vanishes. The method based 
on Dobrushin's estimate presented above for proving the mean field limit applies 
without modification to the convergence of the vortex blob method. See chapter 
5 in [71] for more details on this topic. 

In another variant of these methods, called "the vortex point method" , there 
is no regularization of the interaction kernel, i.e. one takes K^ = K m. the system 
of ODEs above governing the motion of vortices. The vorticity field w = uj{t, x) 
is approximated as follows 

where /i > 0, and the vortex centers satisfy 

N 

xl(t)=Y,ulK{x'^{t)~xl(t)), 

l^k 

xl{0) = hk, k£Z\ 
while the vortex intensities are chosen so that 

uj^ = uj{0,hk) , fceZ^. 

The convergence of the vortex point method for initial data in the Schwartz class 
5(R^) of infinitely differentiable functions with rapidly decaying derivatives of 
all orders has been proved in |45| — see also [46t i86t i48j . The key argument is a 
control of the minimal distance between vortex centers in Proposition 1 of |48| . 

One should also mention recent attempts to justify the derivation of the 
Vlasov-Poisson system as the mean field limit of the Liouville equation for a 
large number of identical point particles with (unmollified) Coulomb interaction: 
see [49|, 150] . These papers prove the mean field limit and the propagation of 
chaos for large systems of point particles with singular interaction force field 
of order 0(d^") for a < 1, where d designates the distance between the two 
interacting particles. 
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1.7.5 Prom particle systems to the Vlasov-Maxwell sys- 
tem 

There are several difficulties in adapting the method for proving the mean field 
limit presented above to the case of the Vlasov-Maxwell system; see [43] for a 
detailed discussion of this problem. 

First, the source term in the system of Maxwell's equations is not a prob- 
ability distribution — as in the case of the Vlasov-Poisson system, where the 
electric field E = E{t, x) is given by 

= -Va;^/(i, a:) , with ~ A^4'f(t,x) ^ pf(t,x) . 

(We recall that 

Pfit,x) := / f{t,x,v)dv 

so that pf(t,-) is a probability distribution on if f{t,-,-) is a probability 
distribution on x R^, which can be assumed without loss of generality since 
the total particle number is invariant under the Vlasov-Poisson dynamics.) 

In the case of the Vlasov-Maxwell system, the source term in Maxwell's 
equations is the 4-vector {pf,jf), defined as follows: 

Pfit,x):= [ f{t,x,Od^, jfit,x):= [ viC)fit,x,OdC . 

JR3 JR3 

This difficulty is handled by an idea introduced in earlier works on Vlasov- 
Maxwell type systems [16l[T5]. The idea is to represent the solution of Maxwell's 
equations for the electromagnetic field in terms of a single momentum distribu- 
tion of electromagnetic potential, as follows. Consider the Cauchy problem for 
the wave equation 

where Dt^a; := df — is the d'Alembert operator in the variables t and x. In 
other words, the momentum variable is a simple parameter in the wave equation 
above. The self-consistent electromagnetic field in the Vlasov-Maxwell system 
is represented as 

E{t,x)=-dt[ v{0uf{t,x,0d^-V^ [ Uf{t,x,Odi 

- dtAo{t, x) - Vx0o(i, x) , 
B{t,x)^iotj; / v{^)uj:{t,x,^)d^ + iotj; Ao{t,x) , 

where (j)o = (j)o{t, x) & R and Aq = Ao{t, x) e R'^ are respectively a scalar and 
a vector potential satisfying 



□t,a:(/>o = 0, and Dt^oj^o = , 
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together with appropriate initial conditions so that the formulas above for E 
and B match the prescribed initial conditions in the Cauchy problem for the 
Vlasov-Maxwell system. 

With this representation for the electromagnetic field in the Vlasov-Maxwell 
system, the source term in the field equation is now / itself, a probability 
distribution in the single particle phase space Ri^ x Rj whenever is a 

probability distribution, since the integral 

is an invariant of the dynamics defined by the Vlasov-Maxwell equations. 

A second difficulty is that the solution of the Cauchy problem for the wave 
equation defining u / involves the space-time convolution with Kirchoff 's kernel 
(the forward fundamental solution of the d'Alembert operator): see for instance 
[57] . This formula is equivalent to the retarded potential formula for solving 
the system of Maxwell's equations. At variance with the solution of the Poisson 
equation 

-A^(pf{t,x) ^ Pf{t,x) , 

the formula giving the solution Uf oi the wave equation above in terms of / 
is not local in the time variable t. Physically, this is due to the fact that the 
electromagnetic field is propagated at the speed of light c > in the Vlasov- 
Maxwell system, while the electrostatic field in the Vlasov-Poisson system is 
propagated instantaneously — in other words, the speed of light is considered 
as infinite in the Vlasov-Poisson system. It remains to check that this difference 
in structure between the Vlasov-Poisson and the Vlasov-Maxwell systems does 
not rule out the possibility of an estimate a la Dobrushin. 

This second difficulty was partially handled in an earlier work |31j for the 
simpler Vlasov-d'Alembert system 

{dt + v{0 ■ '^x)I{t, X, - V,cl,f(t, x) ■ v^fit, x,o^o, 

Ot,x,<l)f{t,x) = Pf{t,x) , 

where p/ is defined in terms of / as above, and v{£^) — + as in the 

Vlasov-Maxwell system. 

A third difficulty — albeit a less essential one — is to choose a regularization 
procedure for the interaction potential that does not destroy the delicate invari- 
ance properties of the Vlasov-Maxwell system. Let Xt = Xc (x) be a regularizing 
sequence in R'^, chosen so that the function Xe is even for each e > 0. The 
Vlasov-Maxwell system is regularized by replacing the momentum distribution 
of electromagnetic potential u / with the solution of the Cauchy problem 

Ot,xUf{t, X, ^) = Xe *x Xe *x f{t, X, ^) , 

where -k^ designates the convolution product in the x variable. This regulariza- 
tion procedure, originally due to E. Horst, is such that both the conservation 
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of the total particle number and the conservation of some variant of the total 
energy are satisfied by the mollified system (see [83]). 

The interested reader is referred to [43] for a complete discussion of the 
material presented in this section. 



1.8 The mean field problem in quantum me- 
chanics 



1.8.1 The A^-body problem in quantum mechanics 



At variance with classical mechanics, the A^-body problem in quantum mechan- 
ics is a PDE, and not a system of ODEs. The state at time t of a system of N 
identical point particles is defined by its iV-body wave function 



We assume that the reader is more or less familiar with the formalism of quan- 
tum mechanics, and we shall not attempt to recall more than a few basic facts. 
An excellent introduction to quantum mechanics can be found in |12| . 



*Ar = ■i'N{t,Xi,...,XN) e C, 



xi, . . . ,xm e • 



The meaning of the wave function is such that . . . , xn)\'^ is the (joint) 

probability density of having particle 1 at the position xi, particle 2 at the 
position X2t ■ ■ ,and particle A^ at the position xn at time t. This implies the 
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normalizatiorF^ 

/ \'^N{t,Xi,...,XN)\'^dxi...dxN 
J(R<i)" 

that is satisfied for all i e R. 

We assume that the interaction between a particle at position x and a particle 
at position y is given by a 2-body potential V{x—y) e R. Henceforth, we assume 
that V is even, so that the force exerted by the particle at position y on the 
particle at position x, i.e. —'S/V{x — y), exactly balances the force exerted by the 
particle at position x on the particle at position y, i.e. —\7V{y — x). (Indeed, 
W is odd since V is even.) 

With these data, we can write the Schrodinger equation governing the A^- 
body wave function for a system of N identical particles of mass m with 2-body 
interaction given by the potential V: 

N 

i;i9t*Ar = -|^^A^,*jv+ V{xk-xi)^N, 

fc=l l<k<l<N 



^^This normalization condition is not satisfied by "generalized eigcnfunctions" of an operator 
with continuous spectrum. Consider the two following examples, where f) = L^(R). 

(a) Let H = + ^^'^ (the quantum harmonic oscillator), which has discrete spectrum 
only. The sequence of eigenvalues of H is n+-^ with n £ N. Besides Ker(H— (n+^)I) = Ch„ 
for each n £ N, with 

h„{x) := ^ \ ^, e-^^/^Hr,{x) , where Hr,(x) = (-l)"e-' -^e""" . 

V2"n!7ri/'' dx" 

The function Hn is the nth Hermite polynomial, and the sequence {hn)n>o ^ Hilbert basis 
of fi. In particular one has the orthonormality condition 

hn{x)hm{x)dx = &mn , for all m, n > 

where Smn is the Kronecker symbol (i.e. Smn = if m ^ n and Smn = 1 if m = ri). 

(b) Let P = (the momentum operator), which has continuous spectrum only. The 
spectrum of P is the real line R, and the generalized eigenfunctions of P are the functions 

: X 1—^ ek{x) := e'-^"'"^. For each fc £ R, one has Pe^ = 27r/cefc but ej. ^ Sj. However one 
has the "formula" analogous to the orthonormality condition in case (a): 

/ ek{x)ei{x)dx = So{k - I) 

where So is the Dirac mass at the origin. The integrand on the left hand side of the equality 
above is not an element of L^(R), and the integral is not a Lebesgue integral. The identity 
above should be understood as the Fourier inversion formula on the class iS'(R) of tempered 
distributions on the real line R. 

In the case (a) , if the system is in an eigenstate corresponding with the eigenvalue n + ^ of 
the operator 'H, its wave function is of the form ^ = u}h„ with \u}\ = 1, so that \\ip\\si = 1- 

In the case (b), if the system is in an eigenstate corresponding with the element k of the 
spectrum of the operator P, it cannot be described by any wave function in ii, but only by a 
generalized eigenfunction of P, that does not belong to i^. 

In the discussion below, we shall never consider quantum states described by generalized 
eigenfunctions as in (b), but only quantum states corresponding with normalized wave func- 
tions. 
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where h is the reduced Planck constant. 

The TV-body Schrodinger equation above is a PDE which is the analogue in 
quantum mechanics of Newton's second law 

{mik = , 
N 
ik = -Y.yV{xk-xi), k = l,...,N 

in classical mechanics. 

The question of existence and uniqueness of a solution of the Schrodinger 
equation of A''-particle quantum dynamics is settled by the following remarkable 
result. (Without loss of generality, we assume that h'^/m =1.) 

Theorem 1.8.1 (Kato) Assume that the space dimension is d — 3. If, for 

some R > 0, 

then, for all N > 1, the unbounded operator 

N 

fc=l l<k<l<N 

has a self-adjoint extension on S^n :~ L'^[{Y{?)^) and generates a unitary group 
g-im« onS^N- 

See chapter V, §5.3 in [60] and [59] for a proof of this result. 

This result is the analogue in quantum mechanics of the existence and 
uniqueness of the solution of the Cauchy problem for the motion equations 
in classical mechanics that follow from the Cauchy-Lipschitz theorem. Observe 
that applying the Cauchy-Lipschitz theorem to the system of ODEs resulting 
from Newton's second law applied to each particle would require V to have Lips- 
chitz continuous first order derivatives. The assumptions on the regularity of the 
potential in Kato's result are obviously much less stringent — for instance, the 
Coulomb potential V{z) = C/\z\ obviously satisfies these assumptions. (That 
the Coulomb potential satisfies the assumptions in Kato's theorem is of course 
very satisfying since most computations in atomic physics involve the Coulomb 
interaction between the electrons and the nuclei in atoms and molecules.) 

As in the case of classical mechanics, our goal is to study the behavior 
of '^N{t) ~ e~**^"^5v i'^ the limit as — ?> oo, under appropriate scaling 
assumptions involving the interaction potential V and the particle number N . 

First we need to define the analogue of the mean field scaling used in the 
case of classical mechanics. Our argument for choosing this scaling is based on 
considering the energy of the particle system. 
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In quantum mechanics, the energy of a system of A'^ identical point particles 
of mass m with pairwise interaction described in terms of the real- valued 2-body 
potential V is 



/ ^\Va:^^N{xi,...,XN)fdXi...,dXN 

+ V / V{Xk-Xi)\'^N{Xl,...,XN)\'^dXi...,dXN, 

where the first term is the kinetic energy, while the second term is the potential 
energy. 

Pick a typical length scale L and energy scale E in the system of N parti- 
cles under consideration. We introduce the dimensionless space variables and 
potential as follows: 

Xk ^ , V{xk - Xl) = — ^(Xfc - Xl) . 

Likewise, the iV-particle wave function is scaled as 

*JV(£, ii, . . . , f jv) = L'^^/^Nit, X-i_,...,XN), 

where T is a time scale to be defined later, and where 

. t 

is the dimensionless time variable. 

Observe that 

j \^ N{i-iX\, . . . ,XN)\^dX\ . . .dxN = j \^ N{t,X\, . . . ,XN)'\^dXi . . .dXN = 

with this scaling. 

With these dimensionless quantities, the kinetic and potential energies be- 
come respectively 



Kinetic energy = / 2;^l^«fe*^l^<^*i • • • '^^^ , 
Potential energy = \J / EV{xk — xi)\"^ ^I'^dxi . . . dxN ■ 



l<k<l<N 



Observe that there are A'' terms in the kinetic energy, while the potential 
energy involves ^N{N — 1) terms. Therefore, we scale the system so that the 
kinetic energy and the potential energy are of same order of magnitude, by 
choosing E and L so that 
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Thus 

N 



2n 

k=l l<k<l<N 



fi2 

_ T-dN/2 " 



At this point, we set 



\ k=l l<k<l<N J 



^ „ 1 



N 

k=l l<k<l<N 

and define the time scale T as follows: 

T "^^^ 



H 

With the rescaled time variable t = t/T, the A^-body Schrodinger equation 
in mean-field scaling becomes 

^ 1 

k=l l<k<l<N 

Henceforth we drop all hats on variables, and consider as our starting point 
the scaled equation 

^ 1 

k=l l<k<l<N 

A formal computation based on the fact that V is real-valued shows that 
d 



dt 

so that 



|*Af(t,a;i, . . .,XN)\^dxi . ..dxN = , 

(R<i)" 



/ |*Ar(t,a;i,...,a;Ar)pdxi...da;Ar = 1. 

■/(R<i)" 

for all t eR.. 

The rigorous argument is based on Kato's theorem (Theorem 11.8. ip stated 
above in space dimension d = 3: with the notation 

^ 1 

k=l l<k<l<N 

the one-parameter group e~**^" is unitary on i^jv = L'^ii'R-'^)^), so that 

/ \'^Nit,xi,...,XN)fdxi...dxN = \\e-''""'i';}\\l^ = 1 
for all t eR. 
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1.8.2 The target mean-field equation 

The target quantum mean-field equation is the quantum analogue of the Vlasov- 
Poisson equation. 

Its unknown is the single-particle wave function tJj = 'ijj{t,x) G C, with 

The basic idea in the mean field approximation is the same as in the clas- 
sical case. Since the probability of finding a particle in an infinitesimal volume 
element dx at time t is \ip{t,x)\'^dx, one expects that the action on the fc-th 
particle located at Xk is defined in terms of the potential V by the formula 

1 ^ r 

-yV{xk-xi)^ V{xk~z)\i;{t,z)\''dz. 



N ^ 

; = i 



This suggests that the target mean-field equation for single-particle wave 
function obtained as the limit of the A'^-body Schrodinger equation is 

idt^[t,x) = -\^AAt.x)+^{t,x) [ V{x-y)\i:{t,y)\^dy. 

This equation is known as the Hartree equation — and also sometimes referred 
to as the Schrodinger-Poisson equation when V{z) is the Coulomb potential in 
space dimension 3. Indeed, if 

= zeR^\{0}, 

the equation above is equivalent to the system 

J idt^{t, x) = -^A^iPit, x) + U{t, x)il^{t, x) , 
\ -A,Uit,x) = m,x)\\ 

This system is the quantum analogue of the Vlasov-Poisson system. 

The Hartree equation has been studied extensively by various authors: see 
[ISIES]- The following statement is a consequence of Theorem 3.1 in [30]. 

Proposition 1.8.2 Assume that V is a real-valued, even element o/L°°(R'*). 
For each i/)'" S i/''(R'*), there exists a unique mild solutioi^jl of the Cauchy 

^^A mild solution of the Cauchy problem 

■ u{t) = Au{t) + F[u{t)] , 



lt=o ' 

where A generates a strongly continuous contraction semigroup on a Banach space X and 
F : X — > X is a continuous map, is an X-valued continuous function / 9 t i— > u{t) G X 
defined on / = [0, t] with r £ [0, +oo] that is a solution of the integral equation 



u(t) = e*^n™ + r e'-^-''^^F[u(s)]ds 
Jo 



for each t a I. 
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problem for Hartree 's equation 

( idttp{t,x) = -^A^i;{t,x) + {v*a:\i^\'^)tp{t,x), xeR'^ 

such that ip e Cb(R; H^{Ii''')) n C(R; iJ'' (R^)). Besides, this solution satisfies 
the conservation laws of mass and energy, viz. 

llV'(ir)||i,2(R<i) = Const. 

and 

5l|VxV'(t,-)|li2(R.) + i // V{x-y)\^{t,x)\^\i){t,y)\^dxdy= Const, 
for all t e R. 

1.8.3 The formalism of density matrices 

Before studying the TV-particle Schrodinger equation in the large A'' limit, we 
need another formulation of the quantum iV-body problem. 

Henceforth we denote by (•!•) (or by (•|-)i5„ to avoid ambiguity if needed) 
the inner product of the Hilbert space j^jv := L^((R'^)^), and by | • | or | • 
the associated norm, defined by the formula 

ms,. ■■= imT.- 

To the A^-particle wave function ^'at = '^n{xi, . . . , xn) that is an element of 
Sjn '■= L'^ii^'^)^), one associates the bounded linear operator Dn on Sj defined 
as follows: 

Dn := orthogonal projection on C^jv in ijAr , 

i.e. 

£)jv$Ar = *Ar / W^'i>N{xu...,XN)dxi...dxN = {'Sn\^n)'^n ■ 

This linear operator is called the "density matrix" of the A''-particle system. 
Whenever convenient, we use the notation common in the physics literature 

Dn ■■= |*iv)(*jv| 

which is the equivalent of 

Dn—'^n^L^s,^,, withLvf^ := (*w|-> 
Obviously, Dn is an integral operator, whose integral kernel is 



Dn^xi,..., xn, yi, . . • , 2/jv) = *jv(a;i, . . . , a;jv)*jv(2/i, • • • , 2/iv) • 
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In the sequel, we shall sometimes abuse the notation and use the same letter 
to designate an integral operator and its integral kernel. 

Being an orthogonal projection in Sjn, the operator Dn is self-adjoint and 
nonnegative: 

Dn ^ D% ^ D%>0. 

This property is the quantum analogue of the positivity of the A^-particle dis- 
tribution function in classical statistical mechanics. 

The density matrix Djsj for a system of identical, indistinguishable particles 
satisfies the following symmetry: for all cr G ©at, one has 

S^Dn = Dn , 

where S'o- is the transformation on integral kernels defined by the formula 

SaDN{xi,. . . , xn, 2/1, . . • , Vn) :=£'jv(a;a-i(i), • • • , a;a-i(jv), ya-i(i)' • • • ' Va-^N)) ■ 
This assumption is not to be confused with the condition 

where U^ is the unitary operator defined onS^N by the formula 

U„-^n{xi,. ■ . ,.TAr) := *Ar(x„-i(i), . . ■,X„-i(M)) ■ 

This assumption is an assumption on the statistics of particles, i.e. on whether 
the particles under consideration follow the Bose-Einstein or the Fermi-Dirac 

statistics. Particles following the Bose-Einstein statistics are called bosons, and 
the wave function of any system of N identical bosons satisfies the condition 

Particles following the Fermi-Dirac statistics are called fermions, and the wave 
function of any system of N identical fermions satisfies the condition 

Exercise: Express Ucra' and U* in terms of Ua and U^' for all cr, cr' G 6n, and 
check that Ua- is a unitary operator on Sjn- Express the transformation in 
terms of for all cr e &n- 
In fact, whenever 

Ua^N = u!{ct)'^n with u!{a) e C and |a;(cr)| = 1 

one has 

S^Dn = Dn . 

In particular, the density matrix of a system of N identical bosons (w(cr) = 1) 
or identical fermions (w(cr) = (— l)^'S'i('^)) satisfies 

S^Dn = Dn . 
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Exercise: Assume that, for some & such that \^n\ ~ 1 and for all 

a £ &N, one has 

Ua'i'N = w(cr)vl'7v with uj{a) e C and \u!{<t)\ = 1 . 

Compute uj{aij') in terms of uj{ij) and for all a^a' G &n- Prove that 

aj(cr) — ±1 for all a € &n- 

Assuming that the time-dependent wave function 4'7v(t, •) of a system of N 
identical particles satisfies the A^-body Schrodinger equation with initial data 
i.e. 



^ 1 



N 

fc=l l<k<l<N 



our next task is to find the equation satisfied by :— |5'Ar)(5'7v|. 

Since e~**^" is a unitary operator on provided that the interaction 
potential V satisfies the assumptions of Kato's theorem above, one has 

where 

:=|^P5^)(vI/»|. 

It is instructive to do the analogous computation in finite dimension (i.e. 
with matrices). Assume that A,B £ Af„(C), and set 

M{t) := e*^Be-*^ . 

Then the function R 9 i i-> Af(t) e A/„(C) is obviously of class C°° (and even 
real analytic) and satisfieJ^ 

dt ^ ' dV ' 

= Ae'^Be-'^ - e'^Be-'^A 
= AM{t) - M{t)A = [A, M{t)] . 

since A commutes with e~*"^. 

In view of this elementary computation, we leave it to the reader to verify 
the following statement. 



^^For linear operators A, B on the Hilbert space io, we designate by [A, B\ their commutator 
[A, B\ := AB — BA. There is obviously a difficulty with the domain of [A, B\ viewed as an 
unbounded operator on f) if A and B are unbounded operators on f). This difficulty will be 
deliberately left aside, as we shall mostly consider B i— > \A, B\ as an unbounded operator on 
C{S)), which is different matter. 
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Proposition 1.8.3 Assume that the potential V satisfies the assumptions in 
Kato's theorem (TheoremlMJD, and let *™ g = L'^{{Ti''^)^). 

Let = "^NityXi, . . . ,xn) be the solution of the N-particle Schrddinger 
equation 

N 



k=l l<k<l<N 

t=0 



z.e. 



^N{t,-) - e-**-f^"*i^ 



N ' 



with scaled Hamiltonian 
Hn 



fc=l l<k<l<N 



where designates the Laplacian in the kth variable, while Vki is the operator 
defined as the multiplication by V{xk — xi). 
Then the density matrix 

DNit) |vI/^(i,.))(*A'(t,-)l =e""''"^Ar(0)e"^" 

satisfies the von Neumann equation 

N 



iDN{t)^[HM,DM{t)]=-\Y}^k.DM{t)] + ^ J2 [Vkl,DN]. 
In terms of integral kernels 



N 

k=l l<k<l<N 



idtDN{t,xi, . . . ,XN,yi, ■ ■ ■ ,2/Jv) 

N 

= -^^i^Xk - ^y^)DN{t,Xi,...,XN,yi,---,yN) 
k=l 

T7 E {V{xk~xi)-V{yk-yi))DN{t,xi,...,XN,yi,---,yN)- 



N 

l<k<l<N 

The fundamental properties of the density matrix are propagated under the 
flow associated to the von Neumann equation. 

Proposition 1.8.4 Assume that the potential V satisfies the assumptions in 
Kato's theorem, and let 

Then for each t G R, one has 

D^it) := e-**-f^"L>i^e'*-f^" = DN{t)* = L'jvW^ > 0, 
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with scaled Hamiltonian 

^ 1 

k=l l<k<l<N 

where Ak designates the Laplacian in the kth variable, while Vki is the operator 
defined as the multiplication by V{xk — xi). 

Likewise, if particles are indistinguishable initially, i.e. if 

then 

S,DM{t) = DN{t) 

for all < e R. 
Proof. First 

DN{t)* = (e~**^"L>i^e^*^")* = e~^*^" (i:>i^)*e**^" 

since Hn is self adjoint, so that (e'*^")* = e"'*^". Then 

DN{tf = (e-**-f^"i:»i^e'*^")2 

= e-**^"(L>i^)2e"-f^" = {e-'*"''D';ie'*""f = DN{t) . 

That SaDN{t) = DN{t) for alU e R if S^D]^ = D'jIJ is a straightforward 
consequence of the following lemma. 

Lemma 1.8.5 Let Ua be the unitary operator defined on S^n for each a g ©^r 
by 

Ua-^N{xi,. . .,xn) ■■= *jv(a;CT-i(i), • ■ .,x^-i(n)) . 

Then 

for each < G R and each a e ©at. 

Thus, for each t € R and each cr € Sjv, one has 

S,DN{t) - u^DN{t)u; - ;7,e-**^«i^™e"^"C/; 

provided that 
■ 

Proof of Lemma [mill If $ e C^{{BJ^)^), one has 
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and 



so that 



Therefore 



TJ_pitHM 



.UHn tj 



for each f e R and each tr e &n- ■ 

1.9 Elements of operator theory 

We recall that S) designates a complex, separable Hilbert space io, with inner 
product denoted by (-I-) and norm denoted by | • |. (To avoid ambiguity, we 
also use the notation (-I-)!, and [ ■ [j--, whenever needed.) Wc always assume 
that the inner product is antilincar in its first argument and linear in its second 
argument. The sot of Hilbcirt basis of — i.e. the set of orthonormal and total 
families in — is denoted by HB{SS). 

The set of bounded operators on 9) is denoted by L{Sti), with operator norm 



We recall that endowed with the operator norm || • || is a Banach algebra 

— and even a C*-algebra, since 



A bounded operator A on is said to be compact if A{B(Q, 1)) is relatively 
compact in ij. The set of compact operators on is denoted by /C(^); one 
easily checks that /C(i5) is a closed (for the operator norm) two-sided ideal in 



We recall the following fundamental facts about compact operators: 

• if A is a compact operator on ij, its adjoint A* is also a compact operator; 

• /C(i5) is the operator-norm closure of the set of finite rank operators on ij. 

Next we introduce the trace norm: for each A e C.{St) 



A bounded operator A on is said to be a trace-class operator if and only if 
it has finite trace norm. The set of trace-class operators on is denoted by 
L^{S^)\ in other words. 



= sup \Ax\ . 

|x| = l 



\\A*A\\ = \\Af for each A e £(io) . 



P||i:= sup ^|(Aefe|/fe)| e [0,oo]. 



C\^)) ■.= {A&C{S)) s.t. Il^lli <oo}. 
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One easily checks that C^{^) is a two-sided ideal in C{Sj), that A n- jlj4||i 
defines a norm on C^{^)), and that equipped with the trace norm || • ||i 

is a separable Banach space. Besides, the map 

^:\^)^ A^ A* gc\^) 

is an isometry for the trace norm. The following properties can be checked 
easily: 

(a) C\^) c JC{^j) C C{S)) with \\A\\ < mill for all A e jC^{Sj), and 

(b) for all A € C^i^) and B e C{^) 

\\AB\\i and ||B^||i < ||A||i||B|| , A e C'{Sj) and B e C{S)) . 

There is a natural notion of trace for operators in infinite dimensional spaces, 
that can be defined for each trace-class operator. For each A G C^{Sj), 

trace(^) := ^(^efejefe) , 

fc>0 

for each (efc)fc>o G HB{Sj). We leave it to the reader to check that the expression 
on the right hand side of the previous equality is independent of the choice of 
(efe)fc>o € HB{Sj). 

Perhaps the most important property of the trace is the following identity: 
for all A e jC.^{Sj) and B e jC{Sj), the operators AB and BA are both trace-class 
operators on Sj, and one has 

trace(AS) = trace(B^) . 

Let us recall the notion of polar decomposition of an operator. For each 
A e £{Sj), there exists \A\,U e £{Sj) such that 

UU* = U*U = I, \A\ = \A\* > and A = \A\U . 

In this decomposition 

\A\ = VAA^. 

Then 

A e <^ \A\ e C\9j) and ||A||i = trace(|A|) . 

The polar decomposition is the analogue of polar coordinates in operator theory: 

\A\ is the analogue of jzj for z G C, while U is the analogue of z/\z\ for z G C*. 
Exercise: Let Ag C'^{Sj). Prove that trace (v/AA*) = trace (v^X*~A) = ||^||i. 

The two following facts are important: 

• £^{Sj) is the (topological) dual oi }C{Sj), with duality defined by the trace as 
follows 

C\S)) X IC{S)) 3 {A, K) ^ iT:ace{AK) G C ; 

• C{Sj) is the (topological) dual of >C^(i5), with duality defined by the trace as 
follows 

C{^) X £i(i5) 9 {B,A) h-j> trace(S^) e C. 
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Another important class of bounded operators on Sj is the class of Hilbert- 
Schmidt operators, denoted by £^{Sj), defined as follows: 

A e C^iSj) ^ A* A e C\S)) ^ AA* e C\i:,) . 

The class C^i^)) is a closed two-sided ideal in C{Sj). It is a Banach space for 
the Hilbert-Schmidt norm defined as follows: 

\\Ah:=\\A*A\\Y' ^\\AA*\\\/\ 

One easily checks that 

£i(io)c/:2(^)c/C(io). 

An important particular case is 55 = i^(A'), where A" is a measurable subset 
of R'', equipped with the Lebesgue measure. In that case, a bounded operator 
A on L'^{X) is a Hilbert-Schmidt operator if and only if it is an integral operator 
of the form 

Af{x) = / a{x, y)f{y)dv for aU / e L\X) , 



with 

In that case 



a e L^{X X X). 



\ 1/2 

\a{x,y)\^dxdy\ = \\a\\L'^(x><x) 



References for the material presented so far in this section are chapter VI of 
[22] and chapter XIX, section 1 of [54]. 

According to the discussion above, each A e £-^(X'^(R")) is an integral 
operator with integral kernel a G L^(R" x R"). The integral kernel of a trace- 
class operator on i^(R") has an additional interesting property, recalled below. 

Lemma 1.9.1 If A E £^(L^(R")), its integral kernel a = a{x,y) is such that 
the maitP^ 

z a{x, X + z) belongs to CbiR"; L^(R")) . 

Exercise. Here is an outline of the proof of Lemma [1.9. II Let A £ with 
^1 :== L^{W). 

1) Prove that the integral kernel a of A can be represented as 

aix, y) = X! ^rien{x)fn{y) 



^^For each topological space X and each vector space i? on R equipped with the norm || ■ \\e, 
we designate by C(X, E) the set of continuous maps from X to E, and we denote 

Ci,(X, E) := {/ e C{X, E) s.t. sup < oo} . 

xex 

In other words, Cb{X, E) is the set of bounded continuous functions defined on X with values 
in E. 
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with e„, /„ € Sj such that 



e^b = i/nb = 1 and with 



< 00 . 



2) Prove that, for each <p G Sj, one has 




as \h\ -> 0. 
3) Prove that 



I \a{x,x + z) — a{x,x + z')\dx ^ as|2; — 2;'|^0 

R" 



and conclude. 

In particular, if A e £^(L^(R")), then its integral kernel a = a{x, y) is such 
that the map x ^ a{x,x) is well-defined as an element of L^(R"), and 



The converse of this last statement is obviously false: if A G C{S)) is an 
integral operator with integral kernel a{x,y), that the map x i— >■ a{x,x) belongs 
to L^R") does not imply in general that A e jC^{L'^{W)). 



Exercise: Construct an example of operator that belongs to £^(L^(R"')) but 
not to /:1(L2(r«)). 



In the sequel, we shall need the notion of "partial trace" for a trace-class 
operator defined on the tensor product of two Hilbert spaces. Let Sji and Sj2 
be two complex, separable Hilbert spaces. 

Definition 1.9.2 For each A e jC^{Sji®Sj2), the first partial trace of A, denoted 

Ai =tracei(A) e£^(i5i), 
is the unique element of C^{S)i) such that 

tracei5j(A:iS) = traccij.^j^^ » /^J) 
holds for each B € C{S)i), where I^^ is the identity on S)2- 

The second partial trace of A is the element of £^(^2) defined analogously. 

Let A be an integral operator with integral kernel a = a{xi, .T2, yi, 2/2) in the 
case where S)j = L^(R'"^) for j = 1,2. In that case, iji (g) i32 = L^(W'^+''^), 
and one easily checks that A-^i is the integral operator with integral kernel 




a{x, x)dx . 
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There is an easy generalization of the notion of partial trace to the case 
of trace-class operators defined on A^-fold tensor products of separable Hilbert 
spaces. Given N complex, separable Hilbert spaces Sji, . . . , S^n and a trace-class 
operator A e C^{Sji . . . (g) i^Ar), for each 1 < fc < A^, we define 

A,k = tracefe(A) 

to be the first partial trace tracci (A) in the decomposition 

iOl (g) . . . <g)SjN — i^l fg • • ■ ® ^k) (g fg • • • (gl^Ar) ■ 

If one views A e C^iSjii^ Sjn) such that A ^ A* > and tracG(A) = 1 
as the operator analogue of an TV-particle probability measure, then A.j^ is the 
analogue of its (first) fc-particle marginal considered above in our description of 
the BBGKY hierarchy for the mean field limit in classical statistical mechanics. 

More details on partial traces can be found in [71 [S] . 

Remark: All the quantum states of the A^-particle system considered so far are 
defined by an A^-particle wave function \E'Ar(i, •) that belongs to i^((R'')^) and 
satisfies the normalization condition ||\['Ar(t, ■)I1l2((R'*)") = 1 for each t e R. An 
A^-particle density matrix D^it) is associated to this iV-particle wave function 
by the formula DM{t) :— |^Ar(t, •))(^jv(^, OIj ^^id we have seen that this A^- 
particle density matrix D^it) satisfies 

DN{t) = DN{t)* = DN{tf , and tracep^vW) = 1 • 

Conversely, if 0^(1) E £^(L^((R'*)^)) satisfies the conditions above, it is an or- 
thogonal projection with rank 1, and therefore there exists $Ar(t, •) G L^((R'^)^) 
such that ||$Ar(t, •)||i2((Rd)jv) = 1 and Dpfit) := |$Ar(t, •))($Ar(t, •)!. Besides, 
the function ^(t, •) is defined uniquely up to multiplication by a complex number 
of modulus 1. 

Such quantum states are referred to as "pure states" . 

However, there also exist more general quantum states, that are described 
by an A^-particle density matrix DN{t) e £^(L^((R'')^)) satisfying 

DN{t) = Djv (0* > , and trace{DN{t)) = 1 . 

The condition 0^(1) ~ DN{t)* > is equivalent to the fact that 

{<i>N\DN{mN) LHiR^r) > for ah $jv e L^HKX) ■ 

Equivalently, the spectrum of 0^(1) consists of nonnegative eigenvalues (since 
DN{t) G £^(L^((R'')^)), we already know that D]\f{t) is compact, so that its 
spectrum consists of eigenvalues only). 

With the normalization condition trace(L)jv(i)) = 1, one sees that -Djv(0 is 
of the form 

^iv(t)-5I^"WI^"(^'-))<^"(*'-)l 

ri>l 
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where {ipn{t, ■))n>i is a complete orthonormal system in L^(R'^)^) and 

A„(i)>0, ^A„(t) = l. 

n>l 

In other words, 0^(1) is a (possibly infinite) convex combination of the density 
matrices \ipn(t, •)){'ipn(t, ■)\ corresponding with pure states. 

While the Schrodinger equation governs the evolution of quantum states 
described by wave functions, the von Neumann equation governs the evolution 
of mixed as well as pure states in quantum mechanics. If -Djv(i) satisfies the 
von Neumann equation 

where 

N 1 

fe=l l<k<l<N 

and V is such that iHpf is the generator of a unitary group on i^(R'*)^) (as in 
Kato's theorem for instance), then the eigenvalues of D^it) satisfy 

A„ (t) = A,i (0) , for each n>l. 

The eigenfunctions of DN{t) are transformed as 

i/,„(i, .) = e"'*^"?/'«(0, •) , for each n > 1 . 

Throughout the following sections, we shall be dealing with pure states, 
although most of our arguments could easily be extended to mixed states. 

The notion of partial trace discussed above leads to a very natural way of 
constructing mixed states, discussed in the exercise below. 

Exercise. Let '^'m ^ L^((R'')*^) be an M-particle wave function such that 
ll*Af||L2((R<i)Af) = 1, and let iV > 1 be such that N < M. Set Dm ■= |1'A/)(*Af| 
and consider the operator ZJjv Dm-.n- 

1) Check that D^ ~ D^j > 0, and that trace(-DAr) = 1. 

2) Assume that Dn = \ip){-<P\ with ip G L'^{(R'^)'^) such that |lV'||L2((R,d)iv) = 1 
and compute trace(I?Ar|V') (V'l) in terms of 4' a/. 

3) Under which condition on vj/j^j does one have Dn = Dj^l (Hint: use the 
result of question 2 and the equality case in the Cauchy-Schwarz inequality to 
find that \1/m is of the form '^Mixi, . . . , xm) = ■ ■ ■ , XN)4'ixN+i, ■ ■ ■ , xm)-) 

This exercise is inspired from the second paragraph in §14 of [64], where it is 
explained that mixed states arise for instance as quantum states for subsystems 
of a larger system in some pure state. The result obtained in the exercise above 
suggests that the state of such a subsystem (particles 1 to N) is pure if and 
only if the subsystem considered (the N first particles are independent from the 
M-N remaining particles.) 
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1.10 The BBGKY hierarchy in quantum me- 
chanics 

In this section, we outline the BBGKY method presented above in the case of 
classical statistical mechanics, and explain how it leads to a rigorous derivation 
of the mean-field Hartree equation from the ./V-particle Schrodinger equation. 

First we introduce some elements of notation. 

For all iV > 1, we set Xjv := (xi, . . . , Xn) and Ijv := (j/i, • ■ • , Vn)', likewise 
for all A'' > 1 and 1 < k < N, we adopt the notation = {zk+i, ■ ■ ■ , zn)- 

Given an integral operator Dm acting on L^{(Rf')^) with integral kernel 
Dm{Xisi,Yisi) we denote by D^-.k the operator with integral kernel 

DN:k{Xk,Yk):= [ DN{Xk,Z%-Yk,Z%)dZ%. 
In terms of operators, 

-Djv:fe = tracefeCfAr) . 
Here again, wc use the same letter to designate an integral operator and its 
integral kernel. This obvious abuse of notation is deliberately chosen to avoid 
the unnecessary multiplication of mathematical symbols. 

1.10.1 The quantum BBGKY hierarchy 

First wc explain how to deduce the BBGKY hierarchy from the A'"-particle 

Schrodinger equation. 

Theorem 1.10.1 Assume that V is a real-valued, even Junction belonging to 
L°°(R''). Let D'j!} ^ [D^SY > fee a trace-class operator with trace(Di^') ^ 1 on 
= io®^ = L'^{{'R'^)^) with Sj = L^(R''), satisfying the indistinguishability 
symmetry, i.e. 

S„D^ = Dj^ for each a e Sjv . 

Let 

DN{t) := e-^*^"i:»™e'*^" , t G R , 

where 

HN = -^,^Ak + - 

fe=l l<k<l<N 
is the N -particle Schrodinger operator. Here designates the Laplacian acting 
on the k-th variable Xk, while Vki designates the multiplication by V{xk — X[). 

Then the sequence of partial traces of Dm satisfies the BBGKY hierarchy of 
differential equations 

^ N k ^ 
iDM:k = — \ "Yy^j^DM-.k] H — y^[^',fe+l? -PjV:fc+l]:fc 

+ ^ [Vlrn,DM:k] 
l<l<m<k 
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for all k > 1, with the convention 

Dn:n = Dn , and D^-.k = whenever k > N . 
In terms of integral kernels, the BBGKY hierarchy reads 

k 

idtDN:k{t,Xk,Yk) = -\Y,i^^i - ^Vi)^N:k{t,Xk,Yu) 

H i^y] / {V{xj - z) -V{yj - z))DN:k+i{t,Xk,z,Yk,z)dz 

m i^i^i-^m)-V{yi-ym))DN:k{t,Xk,Yk) 

l<l<m<k 

for all k> 1. 

Proof. Consider the free A''-body Schrodinger group 

WNit) := exp {-it^{Ai + ... + An)) , 

and set 

fjv(t) := WN{t)DN{t)WN{-t) = WN{t)DN{t)WN{tr . 
Starting from the von Neumann equation, one deduces that 

lFN{t) = ^ [WNmimWN{tr,FN{t)]. 
l<l<m<N 

Since WN{t) is a unitary operator on L^((R'^)^), the operator Fjv(t) is trace- 
class on _L^((R'')^). Since V € L°°(R''), all the operators Vim are bounded on 
L2((R'i)^), so that [WN{t)Vi„^WN{t)*,FN{t)] is trace-class on ^^((Rd^w^ rp^, 
king partial traces of both sides of the equality above leads to 

iFN:k{t) = j^ J2 [WN{t)VimWN{tr,FN{t)]..k 
l<l<m<N 

for all A; > 1. 

Then we must distinguish three cases, 
a) If A; < Z < m < iV, then 

[WN{t)VlmWN{tr,FN{t)]..k = [Vlm{t),FN:m{t)]:k = 0, 

with the notation 

VUt) := WMmmW^ity = e-'4(^'+^'")v;„e^*5(^'+^'") . 

Indeed, let A e C{L^((Ii'^)'')), and denote by Ip the identity in C{L^i(R'^)P)). 
Then 

[A(8l^_k,Vlm{t)]=0, 
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so that 

tvacesjk ( [Vim {t) , FN-.m {t)]:kA)= tracejj^ ( [VJ,„ (t) , Fjv™ (i)] A O 7^_fe) 

= tracej, J[A (g) 7„_fe, V;m(t)]Fjv:m(i)) = , 

with the notation S^k ■= L'^H'R^)''). 
h)Iil<l<m<k<N, then 

[WN{t)VlmWN{tr,FN{t)]..k = [Wk{t)VlmWk{tr,FN:k{t)] • 

Indeed, with the same notation as above, for each A e jC{L^ {{'R'^)'')) , 

tvacesj^{[WN{t)VimWN{t)* , FN{t)]A (g, In^u) 
= tracei5„([yz™(t), FN{t)]A In-u) 
= trace^^ {Vim{t)FN{t)A 7jv-fc) 

- traccjj^ {FN{t)Vim{t)A 7;v-fc) 
= tracej5^(7jv(t)^ lN-kVim{t)) 

- trace^^ {FN{t)Vim{t)A iN-k) 
= tracej5^(Fiv(t)(AVim(i)) iN-k) 

- tvaces^^{FN{t){Vim{t)A) 7jv-fc) 
= tracej5,(i^Ar:fc(f)[A, t//™(t)]) 

= tracei5,([Vi™(t),FAr:fc(i)]A) 

= tracei5,([VFfc(t)V(„H^fc(t)*,i^;v:fcW]^) • 

c) It remains to treat the case l<l<k<m<N. Denote hy a E &n the 
transposition exchanging A: + 1 and m. At this point we use the symmetry of 
both operators and Wn, i-e. the fact that 

S^DNit) = DN{t) and S„WN{t) = Wn^I) 

for each i € R. First, one can check by direct inspection that 

S„[Vi,k+l{t),FN{t)\ = [Vi,m{t),S„FN{t)] = [Vi,m{t),FN{t)] . 

Then, for each A e £(L^((R'^)'^)) and with the same notation as above, 
tTace^,{[Vi,m{t),FN{t)]:kA) = trace^„([V^,„(i),7^jv(t)]^ 

= tvacesj^{[Vl^m{t),FN{t)]Sa{A(g,Im-k)) 
= tTacesj^{Sa[Vl,m{t), FN{t)]Sl{A(g) Im-k)) 
= tTacesj^{Sa[Vl^m{t),FN{t)]{A Im-k)) 

= tTacesj^{[Vi,k+i{t),FN{t)]{A Im-k)) 
= tracers, ([t^,fe+i(t), Fjv(i)]:fe^) 
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(The second equality above follows from the fact that Sa{A(^Ijn-k) = Ac^Im-k 
since a acts as the identity on {1, . . . , k}; the third equality follows from the iden- 
tities trace(5'^T) = trace(C/„TC/*) = trace(r) and S^{BA) = U„BU*U„AUa = 
Sc{B)Sa{A), while the fourth equality is based on the relation = Id.) 
Thus 

[Vl,„.{t).FN{t)]:k = [Vl,k+l{t).FN{t)]..k = [VLk+l{t),FN:k+l{t)]:k 

whenever l<l<k<m<N — the last equality being obvious. 
Therefore 

^ k N ^ 

lFN:k{t) = j:^Yl E [Vl,m{t)^PN{t)]:k + J^ I^- W , ^^JV:fe W] 

(=1 ■m=k+l l<l<m<k 



N-k ^ 



Y,[hk+l{t).FN:k+l{t)]:k + ^ [yim.FN:k{t)]. 



;=i i<;<m<fc 
Undoing the conjugation with VFfe(t), we next arrive at the equality 

k 



N 

3 



It remains to check that 



l<l<m<k 

,=1 V J-.k 



/ -.k 

For each A e £(L^((R'^)'^)) and with the same notation as above, one has 
trace (^(^e-"5^'=+4y,-fc+i,i?Ar:fe+i]e'*^'^^+^^ A^ 

= trace (^e~''^^''+' [Vj.k+i , DN-.k+iV^"^"^'^ A®h 
= trace (^[yj-fc+i,I?W:fe+i]e^*^^''-+^(A(» /i)e-*4^'=+i 

= trace([V,-,fe+i,L'7V:fc+i]^«)/i) 
= trace([V,-,fe+i,i:'jV:fe+i]:feA) 



smce 



and this concludes the proof. ■ 

References for this section are [Ml [Tl [TOl 
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1.10.2 The infinite quantum mean field hierarchy 

From the quantum BBGKY hierarchy, we formaUy deduce an infinite hierarchy 
that is the quantum analogue of the infinite mean field hierarchy in classical 
statistical mechanics. The formal argument is as follows: assume that 

DN:k{t) ^ Dk{t) 

for each fixed A; > 1 as A'' ^ oo in some sense to be made precise, so that 
Then we pass to the limit as iV — >■ oo in both terms 

l<l<m<k 

and 

/ }Vj,k+l,DN:k+l\:k ■ 



0, 



The first term vanishes as N ^ oo since 

^ E [Vlm,D...k]=0('^ 
l<l<m<k ^ 

while the second term satisfies 
N — k ^ 

y^[^-,fc+l;-PjV:fc+l]:fc Ef^.fc+l, -Pfc+l]:fc 

for each fc > 1 as TV oo. 

Therefore, we expect that the sequence of operators {Dk{t))k>i should sat- 
isfy the infinite hierarchy of diflFerential equations 

k k 

Wk = -^^[Aj,Dk] + Y^[Vj,k+i,Dk+i]:k , A: > 1 . 

3=1 3=1 

In terms of integral kernels 

DN:k{t,Xh,Yh) Dh{t,Xh,Yh) , 

for each fc > 1 as A'' — >■ oo, where Dk{t,Xk,Yk) is the integral kernel of the 
operator Dk{t). Then, the sequence of functions {Dk{t, Xk,Yk))k>i satisfies the 
infinite hierarchy of integro-differential equations indexed by fc > 1: 



idtDk{t,Xk,Yk) = -iE(^-. - Ay^)Dk{t,Xk,Yu) 

3=1 

+ [ {V{xj - z) -Vivj - z))Dk+i{t,Xk,z,Yk,z)dz . 
3=1 
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Of course, the argument above is purely formal and remains to be justified 
rigorously. This is achieved by the following theorem. 

Theorem 1.10.2 Assume that the potential V is an even, real-valued function 
such that 

V e C(R'^ \ {0}) n LL(R-'') , V ^0 at oo, andV> -V^ 
for some Vm > 0. Assume that the initial density matrix 

= with [ \^';}{XN)\^dXN^i 

satisfies the indistinguishability symmetry, i.e. 

S^D';} = D'S for each cr G Sat , 

and the energy hound 

N 

^nHH \V^^^^?i{XM)\'dXM 

+ 4 E / V{xi~x^WS{XN)\''dXN^O{N). 



l<;<m<A' 



Assume that, for each k > \ and in the limit N ^ oo, 

D^N-.k C\l\{R'^)'') weak-*. 

Then 

(a) the sequence {DN:k)k>i indexed by N is relatively compact in the infi- 
nite product space ns.>ii°°(R;>C^(i^((R'*)''))), endowed with the product of 
weak-* topologies — L°°(R; >C^(L^((R'')'^))) being the dual of the Banach space 
Li(R;/C(i2((R'^)fc))) equipped with the norm 



\K\\\ := / \\K{t)\\^dt:, 



(h) each limit point {Dk)k>i of this sequence as N oo solves the infinite 
quantum mean field hierarchy written above in the sense of distributions, with 
initial data 

Dk{Q)^Dl\ k>l. 

Although the proof of this result relies on relatively soft functional analytic 
techniques, we shall not give it in detail. We just sketch the key arguments and 
refer the interested reader to [7] for the complete proof. 

The general strategy is to use the uniform bounds on the operator D^lt) 
deduced from the condition trace(£'jv(^)) ~ 1 and the conservation of energy to 
obtain the compactness property in statement (a) of the theorem above. 
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Next, one writes each equation in the BBGKY hierarchjo in the sense of 
distributions and passes to the Hmit as iV — >■ oo for each fc > 1 fixed. This 
strategy uses in particular various observations hsted below. 

A first important ingredient in the proof is the Cauchy-Schwarz inequality 

\DN.,k{t,Xk,Yk:)\^ < DN.,kit,Xk,Xk)DN.,k{t,Yk,Yk) 
since D]\j(t) is a nonnegative operator. It implies in particular the bound 

\DN.,k{t,Xk,Yk)\^dXkdYk < 1 



that is uniform in t € R and N,k > 1. We recall the convention D^-k = 
whenever k > N and the equality 



smce 



/ DN;kit, Zk, Zk)dZk = 1 , 

J(R<i)'= 

/ DN{t, Zn, Zi^)dZ]^ = 1 . 



A second crucial ingredient in the proof is the H^ bound 



1 / \V^,^M{t,XM)\^dXM<^ + ^—V^, 



2 



for each I — 1,...,N. This estimate is a consequence of the initial energy 
bound postulated on together with the conservation of energy under the 
Schrodinger group, implying that 



T-j y. f \V{xk-Xi)\\<fN{t,XN)\''dXN 



l<fe<i<A' •^(^'') 

<£N + UN~mm [ \-9N{t,XN)\^dXN = £n + ~ l)Vm 

■/(R<i)« 

and of the assumption of indistinguishable particles, so that one has, for each 
pair fc, I such that 1 < k I < N, 

f \V,,-^N{t,XN)\^dXN= j |V,,*Ar(t,Xjv)pdXjv. 

J(R'i)" ■/(R'i)" 

A third key step in the proof is based on the following result bearing on 
trace-class operators. 



^"The BBGKY hierarchy has been derived from the A^-particle Schrodinger equation in the 
case where the potential V belongs to L°°(R'*) — see Theorem ll.10.il However its vahdity 
in the sense of distributions under the assumptions of Theorem 11.10.21 results from the same 
arguments as those used in the derivation of the infinite hierarchy and presented below. 
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Lemma 1.10.3 Consider a sequence of integral operators An G £^(L^(R'')) 
with integral kernels a„ = an{x,y) satisfying 

An -^0 m C^{L^(R'^)) weak-* 

as n oo, and 

sup / \an{z, z + h) — an{z, z)\dz — >■ as — >■ . 

n>l jRd 

Then, for each x e Cc{BJ^) 

an{z, z)x{z)dz ^0 as n oo . 



See section 2 of [7] for the proof of this lemma. 

This result is applied in the following context. For each 9 G Cc((R'')'^), 
consider 



ow,fc(t,z,u;) // Q{Xk)Q{Yk)DN:k+i{t,Xn,z,Yk,w)dXkdYk 

J J(R<i)fex(R'i)'= 

Then by the Sobolev estimate above, one has 



\aN,k{t,z,z + h) - aN,k{t,z,z)\dz < |/i|||6|||oo ( — ^ H 



Thus 



as iV — >■ oo, where 



aN.kiz,z)x{z)dz ^ / ak{z,z)x{z)dz 

Rt* JR<i 



J J(R<')'=x(R<')'' 

and where x S C°°(R) is such that < x < 1, with x(^) = 1 for 1/i? < z < i? 
and supp(x) C [1/2R, i? + 1]. This is the key step in obtaining the limit of the 
interaction term 

{V{xk - z) - V{yk - z))DN:k+i{t, Xk, z,Yk, z)dz . 



Some remarks are in order before going further. 

Observe first that the assumptions on the potential V include the (repulsive) 
Coulomb potential in space dimension d — 3 
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between identical charged particles. This case is obviously of fundamental im- 
portance in view of the (many) applications to atomic physics and quantum 
chemistry, in which the interacting particles are electrons. 

Our second observation is of a more mathematical nature. Under the as- 
sumptions on the potential used in the theorem, the interaction integrand 

{V{xj - z) -V{yj - z))Dk+i{t,Xk,z,Yk,z) 

may fail to belong to L^(Rf ) for each (or almost every) Xk, Yk. Yet our analysis 
defines the interaction integral 

/ (Vixj - z) - V{yj - z))Dk+iit,Xk,z,Yk,z)dz 

as a Radon measure in (X„, y„) — instead of a function of these same variables. 
A general reference for the material presented in this section is [7] . 

1.11 The mean field limit in quantum mechanics 
and and Hartree's equation 

1.11.1 Mathematical statement of the mean field limit 

The quantum analogue of the mean field limit obtained in classical statistical 
mechanics is the following result. 

Theorem 1.11.1 Let V be a real-valued, even, bounded measurable function 
defined on R'^, such that V at infinity. Let the initial single-particle wave 
function G H^{'R,'^) satisfy = 1, and let ip = i^itjx) be the solution 

of the Cauchy problem for the Hartree equation 

f [idt + iA,)^(i, x) + ^{t, x) j /{x- z)\''dz = , 

U(0,x) = V"(a;). 

Let ^-i^ = (^*")®^ and let *Ar(t) = e^**^"*]^ be the solution of N -particle 
Schrodinger equation in the mean field scaling, i.e. with 

j=l l<k<l<N 

Then, for each k > 1, one has 

DN:k{t) := \^N{t)){^N{t)\:k Dk{t) 

in i°°(R;£i(L2((R'i)fc))) weak-* as N ^ oo, where 

Dk{t) = mr^){m^^\\ 
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— in other words 

k 

Dk{t,Xk,Yk) ^Y[tP{t,x,)ij{t,yj). 

As already observed in the case of the mean field limit in classical mechanics, 
using the hierarchy approach requires knowing that the mean field PDE is well- 
posed: see Proposition lf .8.21 above. 

Here again, we only sketch the key arguments in the proof, and refer the 
interested reader to [5S] or [7] for a complete account. 

By Theorem ll.il. 11 we know that the sequence {DN;k)k>i indexed by N is 
relatively compact in 

J|l°°(R;/:1(l2((r^)'=))) 

k>l 

equipped with the product topology, where L°°(R; £^(L^((R'')'^))) is endowed 
with the weak-* topology for each fc > 1. 

Besides, Theorem 11.11.11 implies that any limit point {Dk)k>i of that se- 
quence is a solution of the infinite quantum mean field hierarchy with initial 
data 

:= K^/;™ (i/,™)®'^ I , fc>l. 

Proposition 1.11.2 LetV he a real-valued, even, bounded measurable Junction 
defined on R'', such that V ^ Q. Let the initial single-particle wave function 
V-™ G iJi{R'^) satisfy = 1, and let ^ € C(R; i?i(R'^)) nCi(R; ^^(Rrf)) 

be the solution of the Cauchy problem for the Hartree equation 

f {idt + iA^)V'(t, x) + iP{t, x) J ^V{x~ z)\iP{t, z)fdz = , 
U(0,x)=V"(x). 
Then the sequence 

Pk{t) = mt,r'')m,-r''\ 

is a particular solution of the Cauchy problem for the infinite quantum mean 
field hierarchy with initial data 

fc>l. 

Proof. Since?/; e C(R;i/i(R''))nCi(R;i2(Rd)) satisfies the Hartree equation, 
one has 

idti^it,x)i;{t,y)) - -i(A, - Ay)^P{t,x)^{t,y) 

+ ^{t, x)iJ{t~^ f (Vix -z)- Viy - z))|V(t, z)\^dz 
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in the sense of distributions on x x R^. Next, an elementary computation 
shows that, for each k > 1, 

k k 

idtPk{t, Xk,Yk) = i^Yi xi)'4'{t, yi)dt{il>{t, Xj)ip{t, yj)) 
1=1 1=1 

k k 

= -5 I] n - Ay.)i;{t, Xj)i;{t,yj) 

1=1 1=1 

+ y2pk{t,Xk,Yk) I {V{xj-z)-V{yj-z))m,z)\^dz 

k k 
3=1 1=1 

+ y,l {V{xj - z) - V{yj - z))pk{t,Xk,Yk)m,z)\^dz 
i=i 

k 

= -\Y,{^.,-^v,)Pkit^Xk,Yk) 
i=i 

+ X] / (Vi^j - - ^iVj - ^))Pk+i{t,Xk,z,Yk,z)dz . 

Therefore, the sequence {pk{t, Xk,Yk))k>i is a solution of the infinite quantum 
mean field hierarcliy. ■ 

It remains to prove that all limit points (-Dfc)fc>i of the sequence (-Djv:fc)/c>i 
as TV — > oo fall in some uniqueness class for the infinite quantum mean field 
hierarchy. If this is the case, each limit point {Dk)k>i, being a solution of the 
infinite quantum mean field hierarchy with the same initial data as {pk)k>i, 
must coincide with it for all t € R. By compactness and uniqueness of the limit 
point, one concludes that 

DN:k{t) ^ Pk{t) := |V(t,-)®')(V'(i,-)®'=| , for each k > 1 

in the weak-* topology of L°°(R; £i(L2((Rd)fe))) as TV ^ oo. 

Since the integral kernels of the operators pk are of factorized form, i.e. 

fe 

Pk{t,Xk,Yk) = ]J i/'(t,a;j)V'(t,yj) , k>l, 

this last statement is the quantum analogue of the propagation of chaos studied 
above in the context of the mean field limit in classical statistical mechanics. 

In order to carry out this program, two things need to be checked: 
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a) one needs an estimate for the interaction terms 



for J ■ = 1, . . . , fc, and 

b) one needs an estimate for the growth of 



as fc — )■ oo . 



In the next section, we present an abstract framework explaining the role of 
each estimates. 

1.11.2 A tool for studying infinite hierarchies 

We first review a seemingly unrelated result, namely the Nirenberg-Ovcyannikov 
abstract analogue of the Cauchy-Kovalevska theorem. More precisely, we recall 
below a variant of this theorem due to Nishida. 

The Nirenberg-Ovcyannikov abstract Cauchy-Kovalevska theorem 

The Nirenberg-Ovcyannikov abstract Cauchy-Kovalevska theory bears on the 
Cauchy problem for differential equations of the form 



in situations where one can apply neither the Cauchy-Lipschitz theorem nor the 
Peano existence theorem (see chapter 1 in [56|). Such situations arise typically 
when the solution u(t) belongs to function spaces that are not invariant under 
the mapping v i— >■ F(t, v). We first introduce the functional setting and present 
the main result obtained independently by Nirenberg T6 and Ovcyannikov [78] . 

The setting involves a family (a scale) of Banach spaces (i?r)r >o indexed by 
the positive real numbers. The norm in Br is denoted by || • \\r- This scale of 
Banach spaces is nonincreasing, in the sense that, for each r,r' > 0, one has 

(HBr) r' < r ^ Br C Br' with ||w||r' < \\v\\r for each v £ Br ■ 

Next, we specify how the time-dependent "vector field" F behaves with 
respect to the scale of Banach spaces (-Br)r>o- 

Our first assumption bears on the joint dependence of F in both of its 
arguments: 

(a) there exists ro > s.t. 



u{t) ^ F{t,u{t)) , 



{HFa) 



0<r'<r<ro^Fe C(R x Br, Br') 

and Br 3 V ^ F{t, v) E Br' is linear for alH € R . 



A second assumption bears on the dependence of F in its second argument, 
uniformly in the time variable t: 
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(b) there exists C > s.t. 

(HFb) 0<r' <r <ro^\\F{t,v)\\r' <C-^^. 

r — r 

We state only the uniqueness part of the Nirenberg-Ovcyannikov theorem, 
and refer the interested reader to [77 for a the complete statement — that is 
an existence theorem as well, and moreover bears on possibly nonlinear "vector 
fields" F. 

Theorem 1.11.3 For all r e (0,ro) and all a > 0, the only solution u of the 
Cauchy problem 

(u{t)^ Fit,u{t)), 
\u{0) = 0, 

such that u G ((— a, a), Br) is u = 0. 

In order to fully appreciate the meaning of assumption (HFf,), we briefly 
explain how this theorem is related to the classical Cauchy-Kovalevska theorem. 
Consider the PDE with unknown u = u{t, z) G C, 

dtu{t, z) = P{t, z, dz)u{t, z) , 

where t G R and z G C. It is assumed that P is a 1st order linear differential 
operator of the form 

P[t, z, d^)u{t, z) := A{t, z)d^u{t, z) + B{t, z)u{t, z) , 

with coefficients A and B bounded and analytic on (— T, T) x fi^^, where fJ^ C C 
is the open strip defined as follows: 

Ur := {z G C s.t. |Im(z)| < r} . 

Denote by Br the space of bounded holomorphic functions on fi^ with 

\\v\\r := sup \v{z)\ . 

We recall that, for each function / that is holomorphic and bounded in the 
strip Qri one has 

< sup gnj/(z)| z G r!,, and fc > 1 , 

(r — r Y 

because of Cauchy's estimates for holomorphic functions (see for instance chap- 
ter 4, section 2.3 in 3 ). 
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Therefore, one has 



\\P{t,-,d,)u{t,-)\\r' < sup \A{t,z)\ sup \dMt,z)\ 



+ sup \B{t^z)\ sup z)| 




+ sup \B{t,z)\ sup z)| 



zGflrrx zGO,r 



sup^eo, \u{t,z)\ 



r — r 



with 



C := sup \A{t,z)\+ro sup \B{t,z)\. 




In other words, the 1st order differential operator P{t, z, dz) satisfies assumption 
(HF;,) as a consequence of the Cauchy inequahties for bounded holomorphic 
functions in the strip fir- 

Observe that differential operators Q{t, z, dz) of order higher than 1 do not 
satisfy assumption (HFb), but an inequality of the form 



This is a weaker inequality that the one in (HFb), which follows again from 
Cauchy's estimates, this time on higher order derivatives of holomorphic func- 
tions. 

In fact, this is not surprising: the Cauchy-Kovalevska theorem guarantees 
the local existence and uniqueness of a local analytic solution of the Cauchy 
problem for the PDE 



defined for \t\ < T and | Im(z)| < r* for some T > and some r* > 0. This is 
obviously impossible in general if Q is of order higher than 1 . Consider indeed 
the simplest possible example where Q{t,z,dz) = d^. It is well known that 
the Cauchy problem for the backward heat equation is not well-posed in any 
one of the class B^. Indeed, for each t > 0, the semigroup e*^^ maps the set 
of distributions with compact support on R into the set of functions that are 
entire holomorphic (i.e. holomorphic on C). 

Application of the abstract Cauchy-Kovalevska theorem to infinite 
hieretrchies 

In this section, we give an abstract theorem for proving the uniqueness of the 
solution of infinite hierarchies based on the abstract Cauchy-Kovalevska theorem 
presented above. 



\\Q{t,;dz)u{t,-)\\r' <C 



sup^so, \u{t,z)\ 
(r — r'Y 



with k = order of Q . 
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Our setting involves a sequence of Banach spaces (-B„)„>i. For each n > 1, 
the norm in En is denoted by | • |„. 

Along with the Banach spaces En, our setting also involves 

a) a strongly continuous group of isometries J7„(t) defined on En for each n > 1, 

and 

b) a sequence of bounded linear operators Ln,n+i '■ -E'n+i — ^ En defined for each 
n > 1. 

It will be assumed that the sequence of operators Ln,n+i satisfies the follow- 
ing bound: there exists C > such that 

(HL) II 

\\c{E„+i,En) < C'" > for each n > 1 . 

In the setting so defined, we consider the infinite hierarchy of differential 
equations 

Un{t) = Un{t)Ln^n+lUn+l{-t)Un+l{t) , n>l. 

Theorem 1.11.4 Let t* > and let Un G C^{[0,t*], En) for each n > 1 be a 
solution of 

( Un{t) = Un{t)Ln,n+lUn+l{-t)Un+l{t) , n > 1 , 

\w„(0) =0. 

Assume there exists R > such that 
{Hu) sup ||u„(t)||„ < i?" for alln>l. 

0<t<t' 

Then 

Un{t) = for all t e [0, t*] and each n> 1 . 

Proof. We apply the abstract Cauchy-Kovalevska theorem in the following 
setting. First we choose the scale of Banach spaces: for each r > 0, we set 

Br :=lv = {Vn)n>l € J| -En s.t. ||i;||,. := ^r"|'u„|„ < oo > . 

[ n>l n>l J 

Obviously, if < r' < r, one has 
and 

n>l n>l 

for each v E Br- 

The mapping F is defined on sequences v = (un)n>i as follows: 

F{t,v) := {Un{t)Ln,n+lUn+l{-t)Vn+l)n>l ■ 
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With this definition, observing that 

„n+l _ „rn+l 

nr"" < r'" + r'^'^r + . . . + rV^-^ + r" = 



assumption (HL) implies that 

in+l 



\\F{t,v)\\r> <cY,nr'>n\n<cY, 



r 

—, \Vn\n 



n>l n>l 



n>\ 

for each r, r' > such that r' < r, where C is the constant that appears in the 
control of the norm of Ln.n+i in condition (HL). 

This imphes that F satisfies assumption (HF^). As for assumption (HFa), 
it follows from assumption (HFft) and the fact that Un{t) is a linear group of 
isometrics on En for each n > I. 

Thus, if a sequence u{t) = {u{t)„)n>i with u{t)n G En for each n > 1 and 
each t e [0,t*] satisfies the growth condition (Hu) and the infinite hierarchy, 
then u{t) € Br for all t e [0,t*] and for each r e (0, 1/R). 

We deduce from the differential equation that u € C^{[0,t*], Br) for each 
r G {Q,l/R). By the abstract Cauchy-Kovalevska theorem, we conclude that 
u{t) = for each t G [0, t*], which means that 

Unit) = for each t G [0, t*] and each n > 1 . 

This is precisely the expected uniqueness property. ■ 

Notice that assumption (Hu) is essential for the uniqueness property above. 
Here is an easy counterexample, showing that this assumption cannot be dis- 
pensed with. 

Consider the infinite hierarchy of ODEs 

j yk{t) ^ kyk+i{t) , k>l, 

This hierarchy has obvious factorized solutions: if x is the solution of the Riccati 
equation 

r x{t) = x{tf , 

1 x{Q) = x'^ , 



I.e. 



then 



xit) 



1 - tx'"- 
yk{t) = x(t)\ k>l 
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is a solution of the infinite hierarcliy with initial data 

yl^-.^ixn", k>l. 

The infinite hierarchy with — for all fc > has a unique solution 
satisfying the growth condition (Hu), which is yk{t) = for alH e R and k > 1. 
Consider now the function E defined by 



m = 



e-i/* ift>0 
if t < 



that is of class C°° on the real line and satisfies i?(")(0) = for all n > 0. Then 
the formula 

yk{t) = j^^—T:E^''-^'>{t), teR&ndk>l, 
(fc - 1)! 

also defines a solution of the infinite hierarchy. 

That this solution does not grows exponentially as fc — ^ oo corresponds with 
the fact that E is not an analytic function on the line, since it vanishes on the 
half-line without being identically 0, which contradicts the principle of isolated 
zeros. 

References for this section are [SSI E] ■ 

1.11.3 Application to the Hartree limit in the bounded 
potential case 



Finally we explain how to prove Theorem 11.11.11 with the formalism discussed 
above. 

Proof of Theorem 11.11.11 

Step 1: First we write the infinite quantum mean field hierarchy 



in the form 



where 



while 



^jit) = Uj{t)Lj^j+iUj+i{~t)uj+i{t) , j>l 

3 



fc=i 



As before, designates the Laplacian acting on the fcth variable in L^((R'*)-'), 
for each fc = 1, . . . , j, while 

Uj{t) Uj{t)Dj{t) . 
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Step 2: Next we choose appropriate norms and function spaces so as to apply 
Theorem lim Define Ej := (L"^ {(R'^y )) with trace norm denoted by || • ||i 
(or II • 111 J- in cases where the previous notation would lead to ambiguities). 

The operator Uj{t) defined in Step 1 is an isometry of Ej for all j > 1 and 
all t € R, since it consists of conjugating elements of Ej with the unitary group 
associated to the j-particle free Schrodinger equation. Indeed, we recall the 
following elementary result. 

Lemma 1.11.5 Let be a separable Hilhert space. For each A G L^{f)) and 
each unitary operator U on Sj, then UAU* £ 'C^(jo) and one has 

\\UAU*h = Pill- 

Exercise: prove the lemma above. (Hint: use the definition of the trace-norm.) 

Now we need to control the interaction term. 
First, since V G L°°, one has 

\\[Vk.j+i,D,+i]\\i < \\Vk,j+iDj+i\\i + \\D,+iVk,j+i\\i 
< 2\\Vk,j+i\\\\Dj+i\\i < 2\\V\\l^\\Dj+i\\i 

for each {Dj)j>i with Dj e £^ {L^ {{R'^y )) for each j > 1. 
Next, we use the following auxiliary result. 

Lemma 1.11.6 Set Sj = L2(R™+") and i^i = L^{R"'). Let K e C^{S^) be an 
integral operator of the form 

K(f>{x,z)= / k{x,z,y,w)(j){y,w)dydw. 

Let Kz^w be the operator defined on S]i by 

Kz^ni^ix) ^ k{x,z,y,w)il;{y)dy 

for a.e. 2,wGR". Then 

(1) the map h ^ [z ^ K,^,+h] belongs to C{Rl, L^R"^ , C\^^i))); 

(2) for a.e. z e R" the operator K^^z is trace-class on L-^(R™) and 



trace 


/ Kz,zdz 











I trace|if,,^|dz< ||if||i 

^R" 



This result is an amplification of Lemma Fl.O.ll — see Lemma 2.1 in [7]. 
We apply the lemma above to 

K=[Vu,,+,,D,+,], 
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which is an integral operator with integral kernel 

k{Xj,Xj+i,Yj,yj+i) 
= {V{xk ~ Xj+i) ~ - ?/j+i))Dj+i(Xj,Xj+i, Yj-,yj+i) . 



This shows that 



||[yfcj+i,£'j+i]:j||i = trace / K^^^^^^^^^dz^+y 

< ||i^||i<2||T/|U~||i?,+i|li . 

Therefore we conclude that 

\\Lj^J + l\\ciE,^,,E,)<mV\\L^, 

which means that the infinite quantum mean field hierarchy with bounded po- 
tential satisfies assumption (HL) in Theorem 1 1.1 1.41 

Step 3: Finally we control the growth of Dj as j ^ oo, where Dj is a limit point 
oiDN::j in L°°(R+; £i(L2((Rrf)J))) weak-* as TV cx). First we recall that, for 
each > 1 and each t > 0, one has 

DN{t) =DN{t)* > 0, and ||L»Ar(t)||i trace at (t) = 1. 

Applying the lemma above shows that 

\\DN:j{t)\\i < 1 for aU t > 0, N,j > 1. 

(We recall the convention Dn-j = whenever j > N.) Since 

DN:j^Dj{t) in i°°(R+;/:i(i2((R'')J))) weak-*, 

we conclude that 

\\Dj{t)\\i < 1 for alH > 0, j > 1 . 

Step 4: One can apply Theorem II . 1 1 .41 since the infinite hierarchy under con- 
sideration satisfies the bound (HL) with C = 2\\V\\l°°, while each limit point 
Dj of the sequence Dpf.^j in _L°°(R+; £^(L^((R'^)^))) as N oo satisfies the 
growth estimate (Hu) with R = 1. 
Indeed, setting 

v,{t) ■.^u,it)iD,it)-m,r')m,r'\) 

we know that the sequence {vj)j>i is a solution of the infinite mean field hier- 
archy, which satisfies 

\\vM\i,: < \\DM\i,, + \\ m,-r')m,-r mil,, <i + \mr)\\%^2. 
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Since Vj{0) = for each j > 1, we conclude from Theorem 1 1 . 1 1 .41 that Vj{t) = 
for each j > 1 and each t > 0. In other words, 

for each j > 1, which is precisely the desired result. ■ 

Notice that the verification of the growth condition (Hu) on the sequence 
{Dj)j>i is trivial in this problem. This is not always the case, and we shall 
see some examples below, where checking this crucial condition can be a serious 
difficulty. 

References for this section are [551 IZ] • 



1.12 Other mean field limits in quantum me- 
chanics 

1.12.1 Derivation of the Schrodinger-Poisson equation 

The quantum analogue of the Vlasov-Poisson system is the following variant of 
the Hartree equation, where the interaction potential is the repulsive Coulomb 
potential (between particles of equal charges): 

Viz) ^ 



\-k\z\ 

in space dimension d = 3. It takes the form 

' idti>{t, x) = -iAj;V(^> 2;) + x) , a; € , 

Jr3 V - v\ 

or equivalently 

J idt^(t, x) ^ -^Axi){t, x) + U{t, x)ip{t, , x e , 
j -A,Uit,x) - m,x)\\ 

This last form of the mean field equation explains why, as mentioned above, 
the Hartree equation in this case is also known as the "Schrodinger-Poisson" 
equation. 

Although the Coulomb potential is unbounded near the origin, the Schrodin- 
ger-Poisson equation can be derived following the strategy described above: see 
[35] for this very interesting result. In that case, the right choice of spaces En 
is as follows. Define 



Sr-= ^JI~^.,, j-l,...,iV 
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With the same notation as above, the sequence of Banach spaces in Theorem 
II. 11.41 is defined as foUows: 

E^ := {D e C^L^iK'r)) \Si... S^DSi . . . S„ e C\L^{{R^r))} , 

for all n > 1, with norm 

\\D\\E„=tTa.ce\Si...SnDSi...Sn\. 

With this choice of spaces, the inequality in assumption (HL) in Theorem 
1 1.1 1.41 follows from Hardy's inequality 

\2 



I ^dx<C [ \Vfix)\'dx 



(see Lemma 7.1 in |32l). 

The estimate (Hu) is a follows from the conservation laws 

^ trace(i/]^'L»Ar(t)) =0 for all m > , 

where 

N 



N 

fc=l l<k<l<N 

(with Vki denoting the multiplication by V{xk — xi)) and 

DN{t) = |e**-f^"(?/'"')®^)(e**-^"(i/'™)^^| . 

A crucial step in the proof of the growth estimate (Hu) is to compare powers of 
the TV-particle Hamiltonian Hn with powers of the TV-particle free Schrodinger 
operator 



N 



k=l 

— see Propositions 4.1 and 5.1 in [35]. Verifying the growth condition (Hu) is 
perhaps the most technical parts of [35]. 

A more direct proof, by a somewhat different argument, was proposed later 
in [SU] — see section 11.12.41 below. 



1.12.2 Derivation of the nonlinear Schrodinger equation 

One can consider the mean field limit with even more singular interactions, viz. 
V{z) = So{z). Equivalently, one can prove the mean field limit with a potential 
whose range shrinks to zero as the number of particles tends to infinity. Specif- 
ically, one can consider the TV-particle Schrodinger operator in space dimension 
1: 

N N 

Hn := -5 E + 1^ E - xi)) 
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where U >0 belongs to the Schwartz class 5(R), and 7 S (0, 1) . 

The corresponding mean-field equation is the cubic nonlinear Schrodinger 
equation 

idtipit, x) = —^dl'ilj{t, x) + a\^p{t, x)\'^ip{t, x) , a; G R , 

where 

a = U {z)dz . 
Jr 

In this case again, the strategy outlined in Theorem [TTTL4] applies with some 
appropriate choice of norms. The right choice of Banach spaces in this case is 
as follows: 

En {D G C'^iL^R^')) \Si... S^DSi . . . 5„ G CHL^CRn)} , 
with norm 

\\D\\e^ = (trace(|5i...5„7?5i...5„|2))i/2, 
where, as in the previous example, 

S, := ^Jl-d^, j = l,...N. 

Controling the interaction operator in the n-th equation of the infinite mean 
field hierarchy associated to the Schrodinger operator Hn is essentially equiv- 
alent to controling n integral operators on L^(R") with integral kernels of the 
form 

Dn+i{xi, . . . ,Xn,Xi,yi, . . . ,yn,xi) 

in terms of ||-D„+i||,i. 

If Dn is the density matrix associated to a factorized wave function, i.e. if 

Dn{xi,...,Xn,yi,---, Vn) = '/-(s^l) • ■ • '4){Xn)4>{yi) ■ ■ ■ ^j{yn) , 

then Dn e En if and only if G -ff^(R) and ||-Dn||n — HV'll/J'i- In this case, 
£'n+i(a;i, ■ • ■ ,Xn,xi,yi, . . . ,yn,xi) 

= ■lp{xi)\ll'{xi)\'^1p{x2) . ■ . tl){Xn)ll^{yi) ■ ■ ■ IpiVn) 

and the estimate in assumption (HL) reduces to the fact that the Sobolev space 
-ff^(R) is an algebra. (Indeed, we recall that 

u,v€ H\R) =>uve H^R) , 

with 

ll^llffi < C'llwllffi ll'i'llffi for aU u,v e H'^{R) 

for some positive constant C. This elementary fact is most easily checked in 
terms of the Fourier transforms of u and v.) 
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Proving the bound in assumption (HL) for general density matrices (and 
not only in the factorized case) is a rather straightforward generalization of the 
(classical) proof that -ff^(R) is an algebra. Proving the bound in assumption 
(Hu) is somewhat more technical and will not be discussed here. 

We refer the interested reader to [2] — see also [1] — for a detailed statement 
of the result and a complete proof thereof. 

The analogous result in space dimension higher than 1 is considerably more 
involved — notice that iJ^(R'^) is not an algebra whenever d > 2, so that the 
argument above for the bound in assumption (HL) is no longer valid. The proof 
of the mean field limit for interaction potentials shrinking to a Dirac measure 
can be found in [30l [33l [Ml EI]- An interesting approach to the uniqueness 
problem in the corresponding infinite hierarchy can be found in j62l . 

A general reference for this and the previous section is [38]. However Gerard's 
survey in |38| predates Pickl's contributions to the subject, which are described 
in section [1 . 1 2 .41 below. 

1.12.3 The time-dependent Hartree-Fock equations 

In the quantum mean field theories considered so far, the asymptotic A^-particle 
wave function always approached a factorized state of the form 

N 

^Nit,Xi,..., Xn) Y\. ^(*' ^fc) 
fc=l 

in the limit as oo. In fact, we have assumed that the initial state is already 
of this form, i.e. 

N 

*Ar(0,Xi,...,a;Ar) = V'"(2;fe) 

k=l 

in all the situations considered so far. In particular, ^I'jv is symmetric in its 
space variables, i.e. 

*jv(i, a;<T(l), ■ • -jXaiN}) = '^Nit,Xi, . . .,Xn) 

for all (T G &N and all xi,...,xn G R'^- In other words, all the mean field 
theories considered above apply to the case of bosons; the case of fermions 
requires a separate study. 

In the case of fermions, one should replace the initial factorized state with 

*Ar(0,a;i,...,a;jv) = det(V'r(a^O)i<fc,i<iv 

where (V'l, . . . , iPn) is an orthonormal system in L^(R''), i.e. 

/ iJf^)^rix)dx^6ki, k,l = l,...,N. 
Jr.-' 
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This special form of A''-particle wave function is referred to as a Slater determi- 
nant. 

Elementary computations show that Slater determinants satisfy 

/ I det{iljl"'{xi))i<k,i<Nfdxi ...dxN = N\ 

so that *jv(0, . . .) is normahzed in L2((R,rf)JV). 

Moreover, one easily checks that the density matrix associated to a Slater 
determinant has its first marginal given by 



1 ^ 

Dn:i{0) := |*iv(0, . . .))(*iv(0, . . = ^ E IV'f )(V'ri • 



k=l 

In other words, NDn;i{0) is the orthogonal projector on span({^i, . . . ,iPn}), 
so that 

D N:i{0)*= Dn. , i{0) = NDn. ,i{Of>0, and tracc(L>jv,i(0)) = 1 . 

Higher order marginals of the density matrix associated with a Slater determi- 
nant are given in terms of the first marginal by the following expressions: 

DN:n{0) = ^"^^P^' ^iVaW^" E {-^f'''^''^U, 

for all n = 1,. . . ,N, where we recall that, for each permutation ct G 6„, the 
operator Ua is defined by 

Ua^nixi, ...,Xn)= $„(a;^-i(i), . . . , a;^-i(„)) 

for all e L^iiR'')")- 

Consider the A'^-particle Schrodinger operator in the mean field scaling, i.e. 



N 1 ^ 

fe=l l<k<l<N 



for all A'' > 1, where we recall that Vki designates the multiplication by V{xk—xi ) 
acting on ^^((Rd)^). 

The expected mean field equation in the case of fermions is the time-depen- 
dent Hartree-Fock (TDHF) equation, written below in the language of operators 

zF{t) = [-^A,F(t)] + [V,2,F{tr\l - U,2)U . 

Assuming that F{t) is an integral operator on L^(R'') with integral kernel 
F{t, x, y), the equation above can be recast in terms of integral kernels as follows: 

idtF{t, X, y) = -i(A, - Ay)F{t, x, y) 

+ [ {V{x -z)- V{y - z)){F{t, X, y)F{t, z, z) - F{t, z, y)F{t, x, z))dz , 
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or equivalently 

idtF{t,x,y) ^- ^{A,- Ay)F{t,x,y) 



+ F{t,x,y) / {V{x - z) -V{y - z))F{t,z,z)dz 
{V{x - z) - V{y - z))F{t, X, z)F{t, z, y)dz . 

The existence and uniqueness of the solution of the Cauchy problem for the 
TDHF equation is extensively studied in [TOirroilTT] . 

The mean field limit in the case of N fermions with 2-body bounded inter- 
action is summarized in the following statement. 

Theorem 1.12.1 Assume that V is an even function belonging to L°°{IV^), 
and let . . . , ip^^) be an orthonormal system in L^(R''). Let 



where 



N 



and 



^ 1 

k=l l<k<l<N 

*5^'(a;i, . . . , xn) = -^^^ det{'ipk{xi))i<k,i<N ■ 

is the Slater determinant built on the orthonormal system ("01, . . . ,iPn)- 
Then the density operator DN{t) :— {'i' N{t)){^ N(t)\ satisfies 

DN-.nit) - ^"(^- ^ ^ 

m as N ^oo for each n>l, and 

DN:l{t) - FN{t) ^ 

in {L'^ (EV^)) as N oo, where F]\f{t) is the solution of the Cauchy problem 
iFN{t) = i-^A, FN{t)] + [Vi2,FN{t)^^{I~Uu)\:i , 



1 ^ 



k=l 

See [U [H] for proofs of this result and of more general theorems in the same 
direction. 

There is another, more familiar formulation of the TDHF equation. Seek 
the operator F{t) in the form 

1 ^ 
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or equivalently its integral kernel in the form 

F{t,x,y) = — ^■4'k{t,x)'ilJk{t,v), 



N ■ 



where {ilji{t, ■),... ,iljpf{t, ■)) is an orthonormal system in L'^{'R'^). Then the 
functions ipk{t, •) should satisfy the system of PDEs 



idtMt,x)=-lA^Ht,x)+Mt,x) [ V{x- z)^y2\Mt,z)\^dz 

— / V{x- z)^i{t,z)^Pk{t,z)dz, l<k<N. 



N 
k=l 

In this formulation, the functions ipk indexed by fc = 1, . . . , are referred to as 
"molecular orbitals" in the context of quantum chemistry. 

These equations are to be compared with the Hartree equation discussed 
above, i.e. 

idt'il}{t,x) = -^A^ipit^x) +i:{t,x) / V{x- z)\'ip{t,z)fdz. 
Obviously the term 

N 



Mt,x) f V{x~z)^Y.\Mt,z)\'dz 



1=1 

is the analogue in the Hartree- Fock case of the term 

ip{t,x) / V{x ~ z)\i:{t,z)\'^dz 
in Hartree's equation. The additional term in the Hartree-Fock equation, i.e. 

N 

N 

fe=i 

is called the "exchange interaction integral" and is special to the case of fermions: 
see for instance [64] (Problem 1 on p. 233). 

Approaching the A'^-particle wave function of a system of fermions by a 
Slater determinant is obviously the most natural idea. However, other choices 
are also possible — and often made in practice, especially for the purpose of 
numerical computations in the context of quantum chemistry. One such method 
is the theory of multi-configuration time-dependent Hartree-Fock equations [ll], 
where the A^-particle functions is approximated by a linear combination of Slater 
determinants. For the first detailed analysis of the multi-configuration ansatz. 



I N r 

^7 V^i (2;) / V{x- z)ipi (t, z)ipk {t, z)dz 
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which has a very rich mathematical structure, we refer the interested reader to 

Hartree-Fock equations are used primarily in the context of quantum chem- 
istry. Therefore, one has to take into account the fohowing interactions 

(a) electron-electron, 

(b) electron-nuclei, 

(c) nuclei-nuclei. 

In many situations, the nuclei are considered either as fixed (as in [66]), or 
as macroscopic objects governed by a system of ODEs (as in [21]), while the 
quantum description of the electron system involves the Hartree-Fock equations. 
In that case, the only part of the Hamiltonian that is approximated by a non- 
linear, self-consistent mean-field interaction is the repulsive force (a) between 
electrons. 

All these interactions involve the Coulomb potential; the derivation of the 
time dependent Hartree-Fock theory from the A^-particle Schrodinger equation 
with Coulomb interactions remains an open problem at the time of this writing. 

1.12.4 Pickl's approach to quantum mean field fimits 

There are other approaches to the mean field limit in quantum mechanics — 
see for instance the work of Rodnianski and Schlein 85 , which proposes an 
error estimate for the mean field limit, that uses the formalism of Fock spaces 
in second quantization, as well as the more recent references [35( 136] . 

In the present section, we give a brief presentation of yet another approach 
of the mean field limit that avoids the technicalities of BBGKY hierarchies as 
well as the formalism of Fock spaces. The discussion below is based on [80] . 

Pickl's idea is to consider some appropriate quantity that measures the dis- 
tance between the first marginal D]^-i{t) of the iV-particle density matrix Diy{t) 
and the single particle density matrix \tp{t, ■)){ip{t, •)| built on the solution ip of 
the Hartree equation. 

Let "i/N be the solution of the iV-particle Schrodinger equation 

idt^N{t,xi,...,XN) =-ffAr*Ar(i,a;i,...,a;Ar), a;i,...,a;Ar e R'^, 
with 

^ 1 

fc=l l<k<l<N 

Let -ip be the solution of Hartree's equation 

tdt^p{t,x) = -^A,^{t,x) + {V *,\^\^)^lj{t,x) , xeR^ 

Consider the quantity 



£N{t) tmce{DN:i{t){I ~ m, •))(V'(i, ■)!)) • 
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Pickl's derivation of the mean field limit is based on proving that 

Swit) ^ as N —i' oo . 

This estimate implies that DM:i{t) converges to •)) •)! in operator 
norm, in view of the lemma below. 

Henceforth we denote hy S^n the iV-particle Hilbert space S)n i^((R'')^) 
for each iV > 1, and we set 

p{t) 

Lemma 1.12.2 Let N > I and 5" at be a symmetric element of such that 
= 1- -^e^ |*Af)(*Af| and let ip G satisfy = 1- Then 

trace(L'Ar:i(/ - as N ^ oo 

if and only if 

Dn:1 ^ operator norm as N ~> oo . 

This elementary result is statement (a) in Lemma 2.3 of [80j. In some sense, 
it can be considered as a variant of the exercise at the end of section [TTni Notice 
that ^/^(t,-) is a symmetric function of xi, . . . ^ xj^ for all t e R by Lemma 
11.8.51 and because of the choice of the initial data ^'atI^ ^ = )®^, which is 
itself a symmetric function of a: i, ... , xjv- 

Next we seek to control the evolution of £jv(i)- We recall that, for k,l = 
1, . . . ,N, the notation Vki designates the operator on S)n defined as 

(Vfci$jv)(a;i, . . .,xn) V{xk - xi)^n{xi, ■ ■ ■,xn) ■ 

Since trace(DjV:i(0) = 1j oi^^ has 

SNit) - -4trace(Dw:iW/oW) 
dt 

— ~ tra,ce{DN:i{t)p{t)) — tTace{DN:i{t)p{t)) 

= - tTSice{i[^A, DN:i{t)]p{t)) ~ tTace{iDN:i{t)[lA, p{t)]) 

+ tmce{i[Vi2,DN:2{t)]:ip{t}) + tTace{iDN.,i{t)[Vi2, p{t)^%i) . 

First 

[^A,DN:litMt)+DN:lit)[^A,pit)] - [^A,DN:lit)pit)], 

SO that 

trace([iA,Dwa(i)]p(i))+trace(Djv:i(i)[5A,p(t)]) 

= tracc([iA, DN:i{t)p{t)]) = . 
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Next 

trace([yi2,Dw:2(t)]:i/5(i)) = tTeice{[Vi2, DN:2it)]{p{t) (g, I)) , 
tTaceiDN:i{t)[Vi2,p{tr%i) - trace(i^jv:2 W([T42, ^ I)) ■ 

Therefore 

^N{t) = ^trace(z[yi2,l?jv:2(i)]/5(i) ® /) +trace(ii^jV:2W[l^i2,pW®']:i) 
= -i trace(DjV:2 W[n^Vi2 - (T^ * IV-l*, ■)!') ® I, P{t) ® I]) ■ 

The resuh of this computation is summarized in the fonowing lemma. 
Lemma 1.12.3 With the same notation as above, one has 

£N{t) = TN{t) 

where 

TN{t):=2lnitTa.ce{DN:2it){p{t)(E)I)A{{I-p{t)) (g)I)) 

and 

A — i 

N 



A:=^V,2-iV*m,-)\')®I. 



The core of Pickl's argument is the following estimate for J-N{t) in terms of 
£N{t) and N. 

Proposition 1.12.4 Assuming that r > 1 and denoting by r' — the dual 
Holder exponent of r, one has 

i^ivwi < mv\\L^Am.-)\\L^A£N{t) + j^) . 

The proof of Proposition 1 1.1 2 .41 is sketched below. The interested reader is 
referred to the original paper [SD] for more details. 
For the sake of notational simplicity, we set 

p{t) ■.= I-p{t). 

Elementary computations show that 

TN{t) ^ 2lniiia.ce{D N. , 2{t) p{t)'^'^ A{p{t) ® p{t))) 
+ 2 lniiYa.ce{DN:2{t)p{t)®'^Ap{t)®'^) 
+ 2 Imtrace(£)jV:2(i)(pW ® p{t))Ap{t)®'^) . 

The term 

Imtrace(DAr:2(0(p(i) ®p[t))A{p{t) ® p{t))) 

vanishes identically by symmetry of ^Ar(i, •). 

The first term on the right hand side is mastered as follows: observe that 

(/ ® p{t))A{i ® p{t)) = -^{v^m, on ® p{t) 
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— where V * \ip{t, •)p designates the operator on Sji = i^(R'^) defined by 



I tva.ce{DN:2it)p{t)^^A{p{t) p(i)))| 
trace(£)jv:2(i)(p(*) ^ ® ^ ^ 

< \Mt) /)(/ p{t))A{I ^)lll|£»iV:2||l 

< ||(p(t)0 7)(7 0p(i))A(7 0p(t))||. 



Then, observing that A is self-adjoint since ^ is a real-valued, even function, 
one has 



\Mt)^I){I®p{t))AiI^p{t))f 

= \Mt) /)((/ p{t))A{i pm\pit) ® 

= ^\mt,-)\iv*m,-rmtr)) Pit)® Pirn 



since the solution of Hartree's equation satisfies \\tp{t, •)||l2 = 1 for all t gR. 
The third term satisfies 



I tva.ce{DN:2{t){pit) p{t))Ap{tn 
< II (Pit) I)A{p{t) 7) II II (7 p{t))DN..2it){1 pm\ 1 
= ||(p(i) I)A{p{t) 7)11 trace((7 0p(t))DAr:2(f)(/ 0p(t))) 

< \\{p{t) 7)^11 tracc(L'Ar:ip(t)) 
= \\{pit)(^I)A\\£M{t). 



||(p(t)0 7)Af = ||(p(t)0 7)A2(p(t) 7)11 

< (2|iF2|u.+2|i(y*iv(i,-)n'iko iiwt,-)nL.' 

= (2||y||i.. + 2||y*IV(i,-)nii.O IIV'(t,-)lli..' 

<4||F||i..||V(t,-)lli..', 



again because the solution tp of Hartree's equation satisfies ||V'(i, OIU^ = 1 for 
all t e R. 



ci>^{Vi.m,-)\^)cp. 



Therefore 




<j;^\\iv*mt,-)m\Lrmt,-fhr' 
= ^\\v*mt,-r\\hAmr)\\l..' 
<^\\v\\hAmt,-)\\L^r^. 



Indeed, 



Now 
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Controlling the second term in TN{t) is slightly more complicated. 

We shall need the following elements of notation. First we denote by p{t) 
the projection p{t) = /i — p{t) = /i — ■)){tp{t, ■)]. 

For each j ^ 1, . . . , N, we set Ij := /p,,, and for each T e C{Sji), we define 
Tj, jv e C{SjN) for j = 1, . . . , iV as 

Tj^N := Ij^i®T ® In-3 ■ 

For each a G {0, Ij'f^' ' '^^ and each vr e >C(iOi) such that tt* = tt = tt^, we 
define 

TV 

Lemma 1.12.5 Lei 

1 ^ 

MA,(t) = ^-nfc,A,[p(i)] 

w/iere 

ae{o,i}{i. ■■-."> 

a(l) + ... + a(N) = fc 

Besides 

^3,nW)\* = n^- , and Hj- ^[p(t)]nfe,jv[p(t)] - (^.-fcHj- 
/or each j,k — 1, . . . , N , and 

N 
k=0 

In other words, the relation 

^ k 

MAr(t) = ^-nfc,A,[p(i)] 

fe=0 

is the spectral decomposition of the self-adjoint element M^it) of C{S^n)- These 
elementary computations are summarized in formula (6) of |80) . 
Define, for each a > 0, 

N / i.\a 

fc=i ^ ^ 
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and 

N 



fe=i 



that is the pseudo-inverse of M^{t) extended by on Ker(Mjv(f)). In other 
words 

M^(OM^"(t) = In- no,jv[p(<)] 

In particular 

for each j = 1, . . . , N. 

After this observation, the term that remains to be estimated, i.e. 

trace(DiV:2{t)p(i)®'^(/ - 

is recast as 

tTace{DN{t)pi,N{t)p2,N{t)AM^J^{t)M^^^^{t)pi,Nit)p2,Nit)). 

Then 

I trace(£);v {t)pi,N {t)p2,N {t)AM]^'' {t)M-^/^ {t)pi,N {t)P2,N {t))\ 

< ||M-'/'(i)pi,;v(i)P2,Jv(i)*Jv(i, •))lhjv 

X \\Ml/^{t)Api^Nit)p2,N{t)^N{t, -M 

The second factor on the right hand side is mastered by observing that 

\\M'J\t)Ap^,N{t)p2,NmN{t,-))\\L 
= (*Ar(t, ■)\pi^N{t)p2,N{t)AMN{t)Api^N{t)p2^N{t)-^N{t, •))) 

= ■)\piMit)p2.N{t)Ap3,N{t)Api,N{t)p2,Nit)'^N{t, •))> 

+ jf{^N{t, ■)\pl,N{t)p2,N{t)Api^N{t)Api^N{t)p2,N{t)'^N{t, •))> 

since ^'jv(t, •) is symmetric. The first summand is bounded by 

^\\pi,Nm^piMt)\\m.)\\ps,NmNmL 

< \\piMt)^^PhN{t)\\ci^^)£N{t) 

<^v\\l.4^{t,■)\\l.r'SN{t), 

while the second is bounded by 
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Finally, we use again the fact that ^^{t,-) is symmetric to conclude that 

N{N - l)\\M^'/\t)p,,N{t)p2,NmN{t, •))||L 

= 2 {'i'N{t,-),M^\t)p,^N{t)pk^NmN{t,-)) 
l<j<k<N 

< N^^Nit, •), M^Hm^^mNit, •)) 

= N^{^Nit,-),MNmN{t,-)) 
^N^£N{t). 

Putting all these estimates together establishes the bound in Proposition 1 1 . 1 2 .41 
Applying Gronwall's estimate, we deduce from Proposition 1 1 . 1 2 .41 the follow- 
ing result, stated in space dimension 3 for the sake of simplicity. The assump- 
tions on the potential V have been adapted to match those in Kato's Theorem 

[mi 

Theorem 1.12.6 Let V be a real-valued even function such that 

for some r > 1 and some R > 0. 

Let ^™ e L'^CRP) be such that \\ip™\\]^2 = 1, and assume that the Cauchy 
problem for Hartree's equation with initial data tp"^ has a unique solution tp g 
C(R; L^'' (R"^)), where r' = is the dual Holder exponent of r. Let ^ n be 
the solution of the Cauchy problem for the N -particle Schrddinger equation with 
initial data (V'*")®^, and Dm {t) := l^-jvC^, •)>(*Ar(i, •)!• 

Then, for each < G R, one has 

DN:i{t) — > \ip(t, ■)){'>p{t, ■)\ in operator norm as N ^ +oo , 
with the estimate 



tra.ce{DN:i{t){I - m,-))m,-)\)) 



< ^ ( exp 



10\\V\\L2rmS,-)h2,-'ds 



- 1 



Some remarks are in order after this rather remarkable result. 

The first obvious advantage of Pickl's argument is its simplicity. The proof 
of Proposition 11.1 2T4l and the resulting Theorem 11.12.61 provide a much simpler 
derivation of the mean field limit for potentials with a Coulomb singularity at 
the origin than the Erdos-Yau proof [32] based on the BBGKY hierarchy. Ac- 
cording to Theorem 3.1 in '39|, if t/j™ e iJ'°(R'^), then the unique mild solution 
ip of Hartree's equation belongs to C(R+; i/'^(R'^)) under the same assumption 
on V as in Theorem [LlMl Therefore V G C(R; L^''' (R^)) by Sobolev's em- 
bedding theorem provided that k is chosen large enough. In fact, even stronger 
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singularities at the origin than Coulomb in the interaction potential can be 
handled with this method. 

A second advantage of Pickl's method is its versatility. For instance, it can 
also be applied to the derivation of the Gross-Pitaevski: see [5T] . In this latter 
reference, Pickl uses appropriate nonlinear functions of the operator MN{t), 
which are defined in terms of the spectral decomposition provided by Lemma 
11.12.51 The same type of method also applies to study the convergence rate in 
the mean field limit j63j in a more general and systematic way than in Theorem 
11.12.61 above. 

On the other hand, unlike methods based on the BBGKY hierarchy, Pickl's 
argument |80] seems to be limited to pure states and to bosonic particles so 
far. Therefore, in spite of its appealing simplicity, Pickl's method does not 
really supersedes all previous derivations of the mean field limits in quantum 
mechanics. 

Besides, even if BBGKY hierarchies are a rather cumbersome mathematical 
object, they are ubiquitous in the field of nonequilibrium statistical physics, and 
appear in the treatment of quantum as well as classical mechanical models. This 
accounts for the importance of these hierarchies in these lectures. 

1.13 Afterword 

These notes are an updated and (much) expanded version of my earlier (2003) 
lectures |35]. I express my gratitude to my colleagues A. Muntean, J. Radema- 
cher and A. Zagaris for their kind invitation to the NDNS+ 2012 Summer School 
at Universiteit Twente in Enschede. The first part of these notes (sections 1-6) 
were also presented in a joint course with L. Desvillettes at the 37th Summer 
School on Mathematical Physics in Ravello (2012). I am also very grateful to 
Profs. T. Ruggeri and S. Rionero, and to the participants for their positive 
feedback, as well as for the most enjoyable hospitality of the Ravello school. 

As in [42] , the core of the material presented here is the mean field limit of 
A^-particle systems governed by the equations of classical mechanics in the case 
where the iV-particle empirical measure is an exact weak solution of the mean 
field PDE, and the mean field limit of the A^-particle Schrodinger equation. 

Unlike in f42] , I have chosen to introduce both the approach based on empir- 
ical measures and the method based on the BBGKY hierarchy on the mean field 
limit of the A^-particle systems in classical mechanics with Lipschitz continuous 
interaction kernel K. The discussion of the propagation of chaos in section [1.61 
is much more detailed here than in |42j . The connection between the empirical 
measure and the hierarchy approaches is hopefully clarified by the discussion in 
sections [TOl [1711 and [1731 

The most recent developments — in particular results obtained since the 
publication of [42] — are surveyed in sections 11.71 and 11.121 However, I have 
chosen to give less many details on the derivation of the time-dependent Hartree- 
Fock equations from the A^-particle Schrodinger equation. The discussion of 
the time-dependent Hartree-Fock equations in the present notes is limited to a 
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presentation of the equations and to the statement of the main result obtained 
on this problem. My earlier notes [42] contain a sketch of the proof of the mean 
field limit leading to the time-dependent Hartree-Fock equations; the interested 
reader is referred to [42] for more details in this direction. 

I have limited myself to a discussion of evolution problems. Yet there are 
many important results bearing on steady problems: see for instance the work 
of Caglioti-Lions-Marchioro-Pulvirenti [23] and Kiefiling [^ on the Onsager 
statistical theory of vortices in incompressible fluids. I have also neglected to 
discuss the similarities and the differences between the mean field limit of large 
particles systems and the convergence of particle methods for mean field equa- 
tions. The interested reader will find information on these issues in chapter 
5 of [7T] in the context of the vortex formulation of the Euler equations for 
two dimensional incompressible fluids. See also, in the context of Vlasov type 
equations, [TS] [28l and the references therein. 

I have also deliberately chosen to avoid discussing results that would involve 
both the mean field limit (as the particle number N tends to infinity) and 
the classical limit of quantum mechanics (assuming that the typical action of an 
individual particle is large compared to Planck's constant h) . One of the reasons 
for this choice is that this part of the theory is technically more involved and 
perhaps not as complete as what has been presented above. Yet there are many 
important contributions of in this directions, and the interested reader should 
have a look at [ST] [71], to quote only a few references exploring that topic. 

Another conspicuous omission in the second part of these notes is the formal- 
ism of Fock spaces and quantum field theory. This mathematical setting allows 
considering data where the particle number can be infinite — more precisely, 
the particle number N is an unbounded operator in this formalism. The mean 
field limit can be formulated in this setting: see for instance [9] [35l [36] , together 
with the papers [511140] already mentioned above, involving the classical limit in 
addition to the case of infinitely many particles. Interestingly, there is a striking 
analogy between the classical limit of quantum mechanics in formulated in terms 
of Egorov's theorem (Theorem 25.3.5 in |55j ) and the point of view introduced 
in [35l|36] for the quantum mean field limit: see [35]on p. 1024. 

I am indebted in various ways to several colleagues for the material presented 
in these notes. First I wish to thank my collaborators R. Adami, C. Bardos, 
B. Ducomet, L. Erdos, A. Gottlieb, N. Mauser, V. Ricci, A. Teta and H.- 
T. Yau, with whom I have worked on mean field limits on various occasions. 
Besides, the material presented here is in (many) places far from original, and 
the present notes owe much to earlier presentations of mean field limits — such 
as [22] or unpublished notes by C. Villani, as well as lectures by P.-L. Lions 
at the College de France. I also learned a lot from a series of lectures by C. 
Mouhot presenting his joint work with S. Mischler and B. Wennberg [72l[73] on 
quantitative estimates for the propagation of chaos. Most of what I know on 
the derivation of Vlasov type equations from classical particle dynamics comes 
from discussions with J. Batt, H. Neunzert, J. Wick on various occasions in the 
late 80s. I learned of Dobrushin's beautiful estimate [29] from a lecture given 
by M. Pulvirenti at Ecole Normale Superieure in 1997. 
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Another key idea in these lectures is the use of the Nirenberg-Ovcyannikov 
abstract version of the Cauchy-Kovalevska theorem to handle BBGKY type 
hierarchies. I got acquainted with this idea in a seminar given by S. Ukai at 
the Sone seminar in 1998 in Kyoto. Ukai used this technique in a shghtly 
different context (the derivation of the Boltzmann equation from the dynamics 
of a large number of hard spheres in the Boltzmann-Grad limit), and for a 
slightly different purpose (he wanted to simplify the uniform stability estimates 
in Lanford's argument). With his usual modesty, Ukai presented his result as a 
failed attempt at simplifying Lanford's notoriously intricate proof, and did not 
immediately published it — I came across his paper |93j only very recently. In 
the case of the infinite mean field hierarchies for quantum iV-particle dynamics, I 
find the argument based on the abstract variant the Cauchy-Kovalevska theorem 
much more elegant and illuminating than the conventional argument based on 
the Duhamel series as in Lanford's proof, and thus I have chosen to present the 
former in these notes, as in [42]. The argument based on the Duhamel series 
can be found in many places — see for instance the proof of Theorem 5.7 on p. 
610 in ii. 

S. Ukai was at the origin of several fundamental results on the mathematical 
analysis of kinetic models — such as the existence and uniqueness of global solu- 
tions of the Boltzmann equation for all initial data near Maxwellian equilibrium, 
the existence and uniqueness of global solutions of the Vlasov-Poisson system 
[94] in space dimension 2, the hydrodynamic limit of the Boltzmann equation 
leading to the incompressible Navier-Stokes equations in the regime of "small" 
initial data. In view of his infiuence on some of the topics discussed above, these 
notes are dedicated to his memory. 
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